Geant4-11
Data Structures | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | Static Protected Attributes | Private Types
G4GDMLWriteSolids Class Referenceabstract

#include <G4GDMLWriteSolids.hh>

Inheritance diagram for G4GDMLWriteSolids:
G4GDMLWriteMaterials G4GDMLWriteDefine G4GDMLWrite G4GDMLWriteSetup G4GDMLWriteParamvol G4GDMLWriteStructure

Data Structures

class  G4ThreeVectorCompare
 

Public Member Functions

void AddAuxiliary (G4GDMLAuxStructType myaux)
 
void AddElement (const G4Element *const)
 
virtual void AddExtension (xercesc::DOMElement *, const G4LogicalVolume *const)
 
void AddIsotope (const G4Isotope *const)
 
void AddMaterial (const G4Material *const)
 
void AddModule (const G4int depth)
 
void AddModule (const G4VPhysicalVolume *const topVol)
 
void AddPosition (const G4String &name, const G4ThreeVector &pos)
 
virtual void AddSolid (const G4VSolid *const)
 
virtual void DefineWrite (xercesc::DOMElement *)
 
virtual void ExtensionWrite (xercesc::DOMElement *)
 
void FirstpositionWrite (xercesc::DOMElement *element, const G4String &name, const G4ThreeVector &pos)
 
void FirstrotationWrite (xercesc::DOMElement *element, const G4String &name, const G4ThreeVector &rot)
 
G4String GenerateName (const G4String &, const void *const)
 
G4ThreeVector GetAngles (const G4RotationMatrix &)
 
virtual void MaterialsWrite (xercesc::DOMElement *)
 
void PositionWrite (xercesc::DOMElement *element, const G4String &name, const G4ThreeVector &pos)
 
void RotationWrite (xercesc::DOMElement *element, const G4String &name, const G4ThreeVector &rot)
 
void ScaleWrite (xercesc::DOMElement *element, const G4String &name, const G4ThreeVector &scl)
 
void SetOutputFileOverwrite (G4bool flag)
 
virtual void SetupWrite (xercesc::DOMElement *, const G4LogicalVolume *const)=0
 
virtual void SolidsWrite (xercesc::DOMElement *)
 
virtual void StructureWrite (xercesc::DOMElement *)=0
 
virtual void SurfacesWrite ()=0
 
virtual G4Transform3D TraverseVolumeTree (const G4LogicalVolume *const, const G4int)=0
 
virtual void UserinfoWrite (xercesc::DOMElement *)
 
G4Transform3D Write (const G4String &filename, const G4LogicalVolume *const topLog, const G4String &schemaPath, const G4int depth, G4bool storeReferences=true)
 

Static Public Member Functions

static void SetAddPointerToName (G4bool)
 

Protected Member Functions

void AddAuxInfo (G4GDMLAuxListType *auxInfoList, xercesc::DOMElement *element)
 
void AtomWrite (xercesc::DOMElement *, const G4double &)
 
void BooleanWrite (xercesc::DOMElement *, const G4BooleanSolid *const)
 
void BoxWrite (xercesc::DOMElement *, const G4Box *const)
 
void ConeWrite (xercesc::DOMElement *, const G4Cons *const)
 
void CutTubeWrite (xercesc::DOMElement *, const G4CutTubs *const)
 
DepthMapTypeDepthMap ()
 
void DWrite (xercesc::DOMElement *, const G4double &)
 
void ElconeWrite (xercesc::DOMElement *, const G4EllipticalCone *const)
 
void ElementWrite (const G4Element *const)
 
void EllipsoidWrite (xercesc::DOMElement *, const G4Ellipsoid *const)
 
void EltubeWrite (xercesc::DOMElement *, const G4EllipticalTube *const)
 
G4bool FileExists (const G4String &) const
 
 G4GDMLWriteSolids ()
 
void GenericPolyconeWrite (xercesc::DOMElement *, const G4GenericPolycone *const)
 
void GenTrapWrite (xercesc::DOMElement *, const G4GenericTrap *const)
 
void HypeWrite (xercesc::DOMElement *, const G4Hype *const)
 
void IsotopeWrite (const G4Isotope *const)
 
void MaterialWrite (const G4Material *const)
 
void MEEWrite (xercesc::DOMElement *, const G4double &)
 
G4String Modularize (const G4VPhysicalVolume *const topvol, const G4int depth)
 
void MultiUnionWrite (xercesc::DOMElement *solElement, const G4MultiUnion *const)
 
xercesc::DOMAttr * NewAttribute (const G4String &, const G4double &)
 
xercesc::DOMAttr * NewAttribute (const G4String &, const G4String &)
 
xercesc::DOMElement * NewElement (const G4String &)
 
void OpticalSurfaceWrite (xercesc::DOMElement *, const G4OpticalSurface *const)
 
void OrbWrite (xercesc::DOMElement *, const G4Orb *const)
 
void ParaboloidWrite (xercesc::DOMElement *, const G4Paraboloid *const)
 
void ParaWrite (xercesc::DOMElement *, const G4Para *const)
 
void PolyconeWrite (xercesc::DOMElement *, const G4Polycone *const)
 
void PolyhedraWrite (xercesc::DOMElement *, const G4Polyhedra *const)
 
void Position_vectorWrite (xercesc::DOMElement *, const G4String &, const G4String &, const G4ThreeVector &)
 
void PropertyConstWrite (const G4String &, const G4double, const G4MaterialPropertiesTable *)
 
void PropertyVectorWrite (const G4String &, const G4PhysicsFreeVector *const)
 
void PropertyWrite (xercesc::DOMElement *, const G4Material *const)
 
void PropertyWrite (xercesc::DOMElement *, const G4OpticalSurface *const)
 
PhysVolumeMapTypePvolumeMap ()
 
void PWrite (xercesc::DOMElement *, const G4double &)
 
void Rotation_vectorWrite (xercesc::DOMElement *, const G4String &, const G4String &, const G4ThreeVector &)
 
void RZPointWrite (xercesc::DOMElement *, const G4double &, const G4double &)
 
void Scale_vectorWrite (xercesc::DOMElement *, const G4String &, const G4String &, const G4ThreeVector &)
 
void ScaledWrite (xercesc::DOMElement *, const G4ScaledSolid *const)
 
void SphereWrite (xercesc::DOMElement *, const G4Sphere *const)
 
void TessellatedWrite (xercesc::DOMElement *, const G4TessellatedSolid *const)
 
void TetWrite (xercesc::DOMElement *, const G4Tet *const)
 
void TorusWrite (xercesc::DOMElement *, const G4Torus *const)
 
void TrapWrite (xercesc::DOMElement *, const G4Trap *const)
 
void TrdWrite (xercesc::DOMElement *, const G4Trd *const)
 
void TubeWrite (xercesc::DOMElement *, const G4Tubs *const)
 
void TwistedboxWrite (xercesc::DOMElement *, const G4TwistedBox *const)
 
void TwistedtrapWrite (xercesc::DOMElement *, const G4TwistedTrap *const)
 
void TwistedtrdWrite (xercesc::DOMElement *, const G4TwistedTrd *const)
 
void TwistedtubsWrite (xercesc::DOMElement *, const G4TwistedTubs *const)
 
void TWrite (xercesc::DOMElement *, const G4double &)
 
VolumeMapTypeVolumeMap ()
 
void XtruWrite (xercesc::DOMElement *, const G4ExtrudedSolid *const)
 
void ZplaneWrite (xercesc::DOMElement *, const G4double &, const G4double &, const G4double &)
 
virtual ~G4GDMLWriteSolids ()
 

Protected Attributes

G4GDMLAuxListType auxList
 
xercesc::DOMElement * defineElement = nullptr
 
xercesc::DOMDocument * doc = nullptr
 
std::vector< const G4Element * > elementList
 
xercesc::DOMElement * extElement = nullptr
 
std::vector< const G4Isotope * > isotopeList
 
std::vector< const G4Material * > materialList
 
xercesc::DOMElement * materialsElement = nullptr
 
G4bool overwriteOutputFile = false
 
std::vector< const G4PhysicsFreeVector * > propertyList
 
G4String SchemaLocation
 
std::vector< const G4VSolid * > solidList
 
xercesc::DOMElement * solidsElement = nullptr
 
xercesc::DOMElement * userinfoElement = nullptr
 

Static Protected Attributes

static G4bool addPointerToName = true
 
static const G4double kAngularPrecision = DBL_EPSILON
 
static const G4double kLinearPrecision = DBL_EPSILON
 
static const G4double kRelativePrecision = DBL_EPSILON
 
static const G4int maxTransforms = 8
 

Private Types

using DepthMapType = std::map< G4int, G4int >
 
using PhysVolumeMapType = std::map< const G4VPhysicalVolume *, G4String >
 
using VolumeMapType = std::map< const G4LogicalVolume *, G4Transform3D >
 

Detailed Description

Definition at line 73 of file G4GDMLWriteSolids.hh.

Member Typedef Documentation

◆ DepthMapType

using G4GDMLWrite::DepthMapType = std::map<G4int, G4int>
privateinherited

Definition at line 55 of file G4GDMLWrite.hh.

◆ PhysVolumeMapType

using G4GDMLWrite::PhysVolumeMapType = std::map<const G4VPhysicalVolume*, G4String>
privateinherited

Definition at line 54 of file G4GDMLWrite.hh.

◆ VolumeMapType

using G4GDMLWrite::VolumeMapType = std::map<const G4LogicalVolume*, G4Transform3D>
privateinherited

Definition at line 53 of file G4GDMLWrite.hh.

Constructor & Destructor Documentation

◆ G4GDMLWriteSolids()

G4GDMLWriteSolids::G4GDMLWriteSolids ( )
protected

Definition at line 71 of file G4GDMLWriteSolids.cc.

◆ ~G4GDMLWriteSolids()

G4GDMLWriteSolids::~G4GDMLWriteSolids ( )
protectedvirtual

Definition at line 77 of file G4GDMLWriteSolids.cc.

78{
79}

Member Function Documentation

◆ AddAuxiliary()

void G4GDMLWrite::AddAuxiliary ( G4GDMLAuxStructType  myaux)
inherited

Definition at line 398 of file G4GDMLWrite.cc.

399{
400 auxList.push_back(myaux);
401}
G4GDMLAuxListType auxList
Definition: G4GDMLWrite.hh:136

References G4GDMLWrite::auxList.

◆ AddAuxInfo()

void G4GDMLWrite::AddAuxInfo ( G4GDMLAuxListType auxInfoList,
xercesc::DOMElement *  element 
)
protectedinherited

Definition at line 94 of file G4GDMLWrite.cc.

96{
97 for(auto iaux = auxInfoList->cbegin(); iaux != auxInfoList->cend(); ++iaux)
98 {
99 xercesc::DOMElement* auxiliaryElement = NewElement("auxiliary");
100 element->appendChild(auxiliaryElement);
101
102 auxiliaryElement->setAttributeNode(NewAttribute("auxtype", (*iaux).type));
103 auxiliaryElement->setAttributeNode(NewAttribute("auxvalue", (*iaux).value));
104 if(((*iaux).unit) != "")
105 {
106 auxiliaryElement->setAttributeNode(NewAttribute("auxunit", (*iaux).unit));
107 }
108
109 if(iaux->auxList)
110 {
111 AddAuxInfo(iaux->auxList, auxiliaryElement);
112 }
113 }
114 return;
115}
xercesc::DOMElement * NewElement(const G4String &)
Definition: G4GDMLWrite.cc:192
xercesc::DOMAttr * NewAttribute(const G4String &, const G4String &)
Definition: G4GDMLWrite.cc:155
void AddAuxInfo(G4GDMLAuxListType *auxInfoList, xercesc::DOMElement *element)
Definition: G4GDMLWrite.cc:94

References G4GDMLWrite::AddAuxInfo(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by G4GDMLWrite::AddAuxInfo(), G4GDMLWriteStructure::TraverseVolumeTree(), and G4GDMLWrite::UserinfoWrite().

◆ AddElement()

void G4GDMLWriteMaterials::AddElement ( const G4Element * const  elementPtr)
inherited

Definition at line 350 of file G4GDMLWriteMaterials.cc.

351{
352 for(std::size_t i = 0; i < elementList.size(); ++i) // Check if element is
353 { // already in the list!
354 if(elementList[i] == elementPtr)
355 {
356 return;
357 }
358 }
359 elementList.push_back(elementPtr);
360 ElementWrite(elementPtr);
361}
std::vector< const G4Element * > elementList
void ElementWrite(const G4Element *const)

References G4GDMLWriteMaterials::elementList, and G4GDMLWriteMaterials::ElementWrite().

Referenced by G4GDMLWriteMaterials::MaterialWrite().

◆ AddExtension()

void G4GDMLWrite::AddExtension ( xercesc::DOMElement *  ,
const G4LogicalVolume * const   
)
virtualinherited

Definition at line 80 of file G4GDMLWrite.cc.

82{
83 // Empty implementation. To be overwritten by user for specific extensions
84 // related to attributes associated to volumes
85}

Referenced by G4GDMLWriteStructure::TraverseVolumeTree().

◆ AddIsotope()

void G4GDMLWriteMaterials::AddIsotope ( const G4Isotope * const  isotopePtr)
inherited

Definition at line 336 of file G4GDMLWriteMaterials.cc.

337{
338 for(std::size_t i = 0; i < isotopeList.size(); ++i) // Check if isotope is
339 { // already in the list!
340 if(isotopeList[i] == isotopePtr)
341 {
342 return;
343 }
344 }
345 isotopeList.push_back(isotopePtr);
346 IsotopeWrite(isotopePtr);
347}
void IsotopeWrite(const G4Isotope *const)
std::vector< const G4Isotope * > isotopeList

References G4GDMLWriteMaterials::isotopeList, and G4GDMLWriteMaterials::IsotopeWrite().

Referenced by G4GDMLWriteMaterials::ElementWrite().

◆ AddMaterial()

void G4GDMLWriteMaterials::AddMaterial ( const G4Material * const  materialPtr)
inherited

Definition at line 364 of file G4GDMLWriteMaterials.cc.

365{
366 for(std::size_t i = 0; i < materialList.size(); ++i) // Check if material is
367 { // already in the list!
368 if(materialList[i] == materialPtr)
369 {
370 return;
371 }
372 }
373 materialList.push_back(materialPtr);
374 MaterialWrite(materialPtr);
375}
void MaterialWrite(const G4Material *const)
std::vector< const G4Material * > materialList

References G4GDMLWriteMaterials::materialList, and G4GDMLWriteMaterials::MaterialWrite().

Referenced by G4GDMLWriteStructure::TraverseVolumeTree().

◆ AddModule() [1/2]

void G4GDMLWrite::AddModule ( const G4int  depth)
inherited

Definition at line 361 of file G4GDMLWrite.cc.

362{
363 if(depth < 0)
364 {
365 G4Exception("G4GDMLWrite::AddModule()", "InvalidSetup", FatalException,
366 "Depth must be a positive number!");
367 }
368 if(DepthMap().find(depth) != DepthMap().end())
369 {
370 G4Exception("G4GDMLWrite::AddModule()", "InvalidSetup", FatalException,
371 "Adding module(s) at this depth is already requested!");
372 }
373 DepthMap()[depth] = 0;
374}
@ FatalException
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
DepthMapType & DepthMap()
Definition: G4GDMLWrite.cc:73

References G4GDMLWrite::DepthMap(), FatalException, and G4Exception().

◆ AddModule() [2/2]

void G4GDMLWrite::AddModule ( const G4VPhysicalVolume *const  topVol)
inherited

Definition at line 327 of file G4GDMLWrite.cc.

328{
329 G4String fname = GenerateName(physvol->GetName(), physvol);
330 G4cout << "G4GDML: Adding module '" << fname << "'..." << G4endl;
331
332 if(physvol == nullptr)
333 {
334 G4Exception("G4GDMLWrite::AddModule()", "InvalidSetup", FatalException,
335 "Invalid NULL pointer is specified for modularization!");
336 return;
337 }
338 if(dynamic_cast<const G4PVDivision*>(physvol))
339 {
340 G4Exception("G4GDMLWrite::AddModule()", "InvalidSetup", FatalException,
341 "It is not possible to modularize by divisionvol!");
342 return;
343 }
344 if(physvol->IsParameterised())
345 {
346 G4Exception("G4GDMLWrite::AddModule()", "InvalidSetup", FatalException,
347 "It is not possible to modularize by parameterised volume!");
348 return;
349 }
350 if(physvol->IsReplicated())
351 {
352 G4Exception("G4GDMLWrite::AddModule()", "InvalidSetup", FatalException,
353 "It is not possible to modularize by replicated volume!");
354 return;
355 }
356
357 PvolumeMap()[physvol] = fname;
358}
#define G4endl
Definition: G4ios.hh:57
G4GLOB_DLL std::ostream G4cout
G4String GenerateName(const G4String &, const void *const)
Definition: G4GDMLWrite.cc:132
PhysVolumeMapType & PvolumeMap()
Definition: G4GDMLWrite.cc:66
string fname
Definition: test.py:308

References FatalException, test::fname, G4cout, G4endl, G4Exception(), G4GDMLWrite::GenerateName(), G4VPhysicalVolume::GetName(), G4VPhysicalVolume::IsParameterised(), G4VPhysicalVolume::IsReplicated(), and G4GDMLWrite::PvolumeMap().

◆ AddPosition()

void G4GDMLWriteDefine::AddPosition ( const G4String name,
const G4ThreeVector pos 
)
inlineinherited

Definition at line 73 of file G4GDMLWriteDefine.hh.

74 {
76 }
static const G4double pos
xercesc::DOMElement * defineElement
void Position_vectorWrite(xercesc::DOMElement *, const G4String &, const G4String &, const G4ThreeVector &)
const char * name(G4int ptype)

References G4GDMLWriteDefine::defineElement, G4InuclParticleNames::name(), pos, and G4GDMLWriteDefine::Position_vectorWrite().

Referenced by TessellatedWrite(), and TetWrite().

◆ AddSolid()

void G4GDMLWriteSolids::AddSolid ( const G4VSolid * const  solidPtr)
virtual

Definition at line 1137 of file G4GDMLWriteSolids.cc.

1138{
1139 for(std::size_t i = 0; i < solidList.size(); ++i) // Check if solid is
1140 { // already in the list!
1141 if(solidList[i] == solidPtr)
1142 {
1143 return;
1144 }
1145 }
1146
1147 solidList.push_back(solidPtr);
1148
1149 if(const G4BooleanSolid* const booleanPtr =
1150 dynamic_cast<const G4BooleanSolid*>(solidPtr))
1151 {
1152 BooleanWrite(solidsElement, booleanPtr);
1153 }
1154 else if(const G4ScaledSolid* const scaledPtr =
1155 dynamic_cast<const G4ScaledSolid*>(solidPtr))
1156 {
1157 ScaledWrite(solidsElement, scaledPtr);
1158 }
1159 else if(solidPtr->GetEntityType() == "G4MultiUnion")
1160 {
1161 const G4MultiUnion* const munionPtr =
1162 static_cast<const G4MultiUnion*>(solidPtr);
1163 MultiUnionWrite(solidsElement, munionPtr);
1164 }
1165 else if(solidPtr->GetEntityType() == "G4Box")
1166 {
1167 const G4Box* const boxPtr = static_cast<const G4Box*>(solidPtr);
1168 BoxWrite(solidsElement, boxPtr);
1169 }
1170 else if(solidPtr->GetEntityType() == "G4Cons")
1171 {
1172 const G4Cons* const conePtr = static_cast<const G4Cons*>(solidPtr);
1173 ConeWrite(solidsElement, conePtr);
1174 }
1175 else if(solidPtr->GetEntityType() == "G4EllipticalCone")
1176 {
1177 const G4EllipticalCone* const elconePtr =
1178 static_cast<const G4EllipticalCone*>(solidPtr);
1179 ElconeWrite(solidsElement, elconePtr);
1180 }
1181 else if(solidPtr->GetEntityType() == "G4Ellipsoid")
1182 {
1183 const G4Ellipsoid* const ellipsoidPtr =
1184 static_cast<const G4Ellipsoid*>(solidPtr);
1185 EllipsoidWrite(solidsElement, ellipsoidPtr);
1186 }
1187 else if(solidPtr->GetEntityType() == "G4EllipticalTube")
1188 {
1189 const G4EllipticalTube* const eltubePtr =
1190 static_cast<const G4EllipticalTube*>(solidPtr);
1191 EltubeWrite(solidsElement, eltubePtr);
1192 }
1193 else if(solidPtr->GetEntityType() == "G4ExtrudedSolid")
1194 {
1195 const G4ExtrudedSolid* const xtruPtr =
1196 static_cast<const G4ExtrudedSolid*>(solidPtr);
1197 XtruWrite(solidsElement, xtruPtr);
1198 }
1199 else if(solidPtr->GetEntityType() == "G4Hype")
1200 {
1201 const G4Hype* const hypePtr = static_cast<const G4Hype*>(solidPtr);
1202 HypeWrite(solidsElement, hypePtr);
1203 }
1204 else if(solidPtr->GetEntityType() == "G4Orb")
1205 {
1206 const G4Orb* const orbPtr = static_cast<const G4Orb*>(solidPtr);
1207 OrbWrite(solidsElement, orbPtr);
1208 }
1209 else if(solidPtr->GetEntityType() == "G4Para")
1210 {
1211 const G4Para* const paraPtr = static_cast<const G4Para*>(solidPtr);
1212 ParaWrite(solidsElement, paraPtr);
1213 }
1214 else if(solidPtr->GetEntityType() == "G4Paraboloid")
1215 {
1216 const G4Paraboloid* const paraboloidPtr =
1217 static_cast<const G4Paraboloid*>(solidPtr);
1218 ParaboloidWrite(solidsElement, paraboloidPtr);
1219 }
1220 else if(solidPtr->GetEntityType() == "G4Polycone")
1221 {
1222 const G4Polycone* const polyconePtr =
1223 static_cast<const G4Polycone*>(solidPtr);
1224 PolyconeWrite(solidsElement, polyconePtr);
1225 }
1226 else if(solidPtr->GetEntityType() == "G4GenericPolycone")
1227 {
1228 const G4GenericPolycone* const genpolyconePtr =
1229 static_cast<const G4GenericPolycone*>(solidPtr);
1230 GenericPolyconeWrite(solidsElement, genpolyconePtr);
1231 }
1232 else if(solidPtr->GetEntityType() == "G4Polyhedra")
1233 {
1234 const G4Polyhedra* const polyhedraPtr =
1235 static_cast<const G4Polyhedra*>(solidPtr);
1236 PolyhedraWrite(solidsElement, polyhedraPtr);
1237 }
1238 else if(solidPtr->GetEntityType() == "G4Sphere")
1239 {
1240 const G4Sphere* const spherePtr = static_cast<const G4Sphere*>(solidPtr);
1241 SphereWrite(solidsElement, spherePtr);
1242 }
1243 else if(solidPtr->GetEntityType() == "G4TessellatedSolid")
1244 {
1245 const G4TessellatedSolid* const tessellatedPtr =
1246 static_cast<const G4TessellatedSolid*>(solidPtr);
1247 TessellatedWrite(solidsElement, tessellatedPtr);
1248 }
1249 else if(solidPtr->GetEntityType() == "G4Tet")
1250 {
1251 const G4Tet* const tetPtr = static_cast<const G4Tet*>(solidPtr);
1252 TetWrite(solidsElement, tetPtr);
1253 }
1254 else if(solidPtr->GetEntityType() == "G4Torus")
1255 {
1256 const G4Torus* const torusPtr = static_cast<const G4Torus*>(solidPtr);
1257 TorusWrite(solidsElement, torusPtr);
1258 }
1259 else if(solidPtr->GetEntityType() == "G4GenericTrap")
1260 {
1261 const G4GenericTrap* const gtrapPtr =
1262 static_cast<const G4GenericTrap*>(solidPtr);
1263 GenTrapWrite(solidsElement, gtrapPtr);
1264 }
1265 else if(solidPtr->GetEntityType() == "G4Trap")
1266 {
1267 const G4Trap* const trapPtr = static_cast<const G4Trap*>(solidPtr);
1268 TrapWrite(solidsElement, trapPtr);
1269 }
1270 else if(solidPtr->GetEntityType() == "G4Trd")
1271 {
1272 const G4Trd* const trdPtr = static_cast<const G4Trd*>(solidPtr);
1273 TrdWrite(solidsElement, trdPtr);
1274 }
1275 else if(solidPtr->GetEntityType() == "G4Tubs")
1276 {
1277 const G4Tubs* const tubePtr = static_cast<const G4Tubs*>(solidPtr);
1278 TubeWrite(solidsElement, tubePtr);
1279 }
1280 else if(solidPtr->GetEntityType() == "G4CutTubs")
1281 {
1282 const G4CutTubs* const cuttubePtr = static_cast<const G4CutTubs*>(solidPtr);
1283 CutTubeWrite(solidsElement, cuttubePtr);
1284 }
1285 else if(solidPtr->GetEntityType() == "G4TwistedBox")
1286 {
1287 const G4TwistedBox* const twistedboxPtr =
1288 static_cast<const G4TwistedBox*>(solidPtr);
1289 TwistedboxWrite(solidsElement, twistedboxPtr);
1290 }
1291 else if(solidPtr->GetEntityType() == "G4TwistedTrap")
1292 {
1293 const G4TwistedTrap* const twistedtrapPtr =
1294 static_cast<const G4TwistedTrap*>(solidPtr);
1295 TwistedtrapWrite(solidsElement, twistedtrapPtr);
1296 }
1297 else if(solidPtr->GetEntityType() == "G4TwistedTrd")
1298 {
1299 const G4TwistedTrd* const twistedtrdPtr =
1300 static_cast<const G4TwistedTrd*>(solidPtr);
1301 TwistedtrdWrite(solidsElement, twistedtrdPtr);
1302 }
1303 else if(solidPtr->GetEntityType() == "G4TwistedTubs")
1304 {
1305 const G4TwistedTubs* const twistedtubsPtr =
1306 static_cast<const G4TwistedTubs*>(solidPtr);
1307 TwistedtubsWrite(solidsElement, twistedtubsPtr);
1308 }
1309 else
1310 {
1311 G4String error_msg = "Unknown solid: " + solidPtr->GetName() +
1312 "; Type: " + solidPtr->GetEntityType();
1313 G4Exception("G4GDMLWriteSolids::AddSolid()", "WriteError", FatalException,
1314 error_msg);
1315 }
1316}
Definition: G4Box.hh:56
Definition: G4Cons.hh:78
void TwistedtrdWrite(xercesc::DOMElement *, const G4TwistedTrd *const)
void TorusWrite(xercesc::DOMElement *, const G4Torus *const)
void TetWrite(xercesc::DOMElement *, const G4Tet *const)
void TrapWrite(xercesc::DOMElement *, const G4Trap *const)
void HypeWrite(xercesc::DOMElement *, const G4Hype *const)
void TwistedtrapWrite(xercesc::DOMElement *, const G4TwistedTrap *const)
void ScaledWrite(xercesc::DOMElement *, const G4ScaledSolid *const)
void MultiUnionWrite(xercesc::DOMElement *solElement, const G4MultiUnion *const)
void PolyhedraWrite(xercesc::DOMElement *, const G4Polyhedra *const)
void TessellatedWrite(xercesc::DOMElement *, const G4TessellatedSolid *const)
void GenTrapWrite(xercesc::DOMElement *, const G4GenericTrap *const)
void TwistedboxWrite(xercesc::DOMElement *, const G4TwistedBox *const)
void TubeWrite(xercesc::DOMElement *, const G4Tubs *const)
void ParaboloidWrite(xercesc::DOMElement *, const G4Paraboloid *const)
std::vector< const G4VSolid * > solidList
void SphereWrite(xercesc::DOMElement *, const G4Sphere *const)
void BoxWrite(xercesc::DOMElement *, const G4Box *const)
void OrbWrite(xercesc::DOMElement *, const G4Orb *const)
void TwistedtubsWrite(xercesc::DOMElement *, const G4TwistedTubs *const)
void EltubeWrite(xercesc::DOMElement *, const G4EllipticalTube *const)
void BooleanWrite(xercesc::DOMElement *, const G4BooleanSolid *const)
void ElconeWrite(xercesc::DOMElement *, const G4EllipticalCone *const)
void ConeWrite(xercesc::DOMElement *, const G4Cons *const)
void GenericPolyconeWrite(xercesc::DOMElement *, const G4GenericPolycone *const)
xercesc::DOMElement * solidsElement
void CutTubeWrite(xercesc::DOMElement *, const G4CutTubs *const)
void EllipsoidWrite(xercesc::DOMElement *, const G4Ellipsoid *const)
void PolyconeWrite(xercesc::DOMElement *, const G4Polycone *const)
void ParaWrite(xercesc::DOMElement *, const G4Para *const)
void TrdWrite(xercesc::DOMElement *, const G4Trd *const)
void XtruWrite(xercesc::DOMElement *, const G4ExtrudedSolid *const)
Definition: G4Hype.hh:69
Definition: G4Orb.hh:56
Definition: G4Para.hh:79
Definition: G4Tet.hh:56
Definition: G4Trd.hh:63
Definition: G4Tubs.hh:75
G4String GetName() const
virtual G4GeometryType GetEntityType() const =0

References BooleanWrite(), BoxWrite(), ConeWrite(), CutTubeWrite(), ElconeWrite(), EllipsoidWrite(), EltubeWrite(), FatalException, G4Exception(), GenericPolyconeWrite(), GenTrapWrite(), G4VSolid::GetEntityType(), G4VSolid::GetName(), HypeWrite(), MultiUnionWrite(), OrbWrite(), ParaboloidWrite(), ParaWrite(), PolyconeWrite(), PolyhedraWrite(), ScaledWrite(), solidList, solidsElement, SphereWrite(), TessellatedWrite(), TetWrite(), TorusWrite(), TrapWrite(), TrdWrite(), TubeWrite(), TwistedboxWrite(), TwistedtrapWrite(), TwistedtrdWrite(), TwistedtubsWrite(), and XtruWrite().

Referenced by BooleanWrite(), MultiUnionWrite(), ScaledWrite(), and G4GDMLWriteStructure::TraverseVolumeTree().

◆ AtomWrite()

void G4GDMLWriteMaterials::AtomWrite ( xercesc::DOMElement *  element,
const G4double a 
)
protectedinherited

Definition at line 53 of file G4GDMLWriteMaterials.cc.

55{
56 xercesc::DOMElement* atomElement = NewElement("atom");
57 atomElement->setAttributeNode(NewAttribute("unit", "g/mole"));
58 atomElement->setAttributeNode(NewAttribute("value", a * mole / g));
59 element->appendChild(atomElement);
60}
static constexpr double mole
Definition: G4SIunits.hh:279
static constexpr double g
Definition: G4SIunits.hh:168

References g, mole, G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by G4GDMLWriteMaterials::ElementWrite(), G4GDMLWriteMaterials::IsotopeWrite(), and G4GDMLWriteMaterials::MaterialWrite().

◆ BooleanWrite()

void G4GDMLWriteSolids::BooleanWrite ( xercesc::DOMElement *  solElement,
const G4BooleanSolid * const  boolean 
)
protected

Definition at line 141 of file G4GDMLWriteSolids.cc.

143{
144 G4int displaced = 0;
145
146 G4String tag("undefined");
147 if(dynamic_cast<const G4IntersectionSolid*>(boolean))
148 {
149 tag = "intersection";
150 }
151 else if(dynamic_cast<const G4SubtractionSolid*>(boolean))
152 {
153 tag = "subtraction";
154 }
155 else if(dynamic_cast<const G4UnionSolid*>(boolean))
156 {
157 tag = "union";
158 }
159
160 G4VSolid* firstPtr = const_cast<G4VSolid*>(boolean->GetConstituentSolid(0));
161 G4VSolid* secondPtr = const_cast<G4VSolid*>(boolean->GetConstituentSolid(1));
162
163 G4ThreeVector firstpos, firstrot, pos, rot;
164
165 // Solve possible displacement of referenced solids!
166 //
167 while(true)
168 {
169 if(displaced > 8)
170 {
171 G4String ErrorMessage = "The referenced solid '" + firstPtr->GetName() +
172 +"in the Boolean shape '" + +boolean->GetName() +
173 +"' was displaced too many times!";
174 G4Exception("G4GDMLWriteSolids::BooleanWrite()", "InvalidSetup",
175 FatalException, ErrorMessage);
176 }
177
178 if(G4DisplacedSolid* disp = dynamic_cast<G4DisplacedSolid*>(firstPtr))
179 {
180 firstpos += disp->GetObjectTranslation();
181 firstrot += GetAngles(disp->GetObjectRotation());
182 firstPtr = disp->GetConstituentMovedSolid();
183 ++displaced;
184 continue;
185 }
186 break;
187 }
188 displaced = 0;
189 while(true)
190 {
191 if(displaced > maxTransforms)
192 {
193 G4String ErrorMessage = "The referenced solid '" + secondPtr->GetName() +
194 +"in the Boolean shape '" + +boolean->GetName() +
195 +"' was displaced too many times!";
196 G4Exception("G4GDMLWriteSolids::BooleanWrite()", "InvalidSetup",
197 FatalException, ErrorMessage);
198 }
199
200 if(G4DisplacedSolid* disp = dynamic_cast<G4DisplacedSolid*>(secondPtr))
201 {
202 pos += disp->GetObjectTranslation();
203 rot += GetAngles(disp->GetObjectRotation());
204 secondPtr = disp->GetConstituentMovedSolid();
205 ++displaced;
206 continue;
207 }
208 break;
209 }
210
211 AddSolid(firstPtr); // At first add the constituent solids!
212 AddSolid(secondPtr);
213
214 const G4String& name = GenerateName(boolean->GetName(), boolean);
215 const G4String& firstref = GenerateName(firstPtr->GetName(), firstPtr);
216 const G4String& secondref = GenerateName(secondPtr->GetName(), secondPtr);
217
218 xercesc::DOMElement* booleanElement = NewElement(tag);
219 booleanElement->setAttributeNode(NewAttribute("name", name));
220 xercesc::DOMElement* firstElement = NewElement("first");
221 firstElement->setAttributeNode(NewAttribute("ref", firstref));
222 booleanElement->appendChild(firstElement);
223 xercesc::DOMElement* secondElement = NewElement("second");
224 secondElement->setAttributeNode(NewAttribute("ref", secondref));
225 booleanElement->appendChild(secondElement);
226 solElement->appendChild(booleanElement);
227 // Add the boolean solid AFTER the constituent solids!
228
229 if((std::fabs(pos.x()) > kLinearPrecision) ||
230 (std::fabs(pos.y()) > kLinearPrecision) ||
231 (std::fabs(pos.z()) > kLinearPrecision))
232 {
233 PositionWrite(booleanElement, name + "_pos", pos);
234 }
235
236 if((std::fabs(rot.x()) > kAngularPrecision) ||
237 (std::fabs(rot.y()) > kAngularPrecision) ||
238 (std::fabs(rot.z()) > kAngularPrecision))
239 {
240 RotationWrite(booleanElement, name + "_rot", rot);
241 }
242
243 if((std::fabs(firstpos.x()) > kLinearPrecision) ||
244 (std::fabs(firstpos.y()) > kLinearPrecision) ||
245 (std::fabs(firstpos.z()) > kLinearPrecision))
246 {
247 FirstpositionWrite(booleanElement, name + "_fpos", firstpos);
248 }
249
250 if((std::fabs(firstrot.x()) > kAngularPrecision) ||
251 (std::fabs(firstrot.y()) > kAngularPrecision) ||
252 (std::fabs(firstrot.z()) > kAngularPrecision))
253 {
254 FirstrotationWrite(booleanElement, name + "_frot", firstrot);
255 }
256}
int G4int
Definition: G4Types.hh:85
double z() const
double x() const
double y() const
void FirstrotationWrite(xercesc::DOMElement *element, const G4String &name, const G4ThreeVector &rot)
static const G4double kLinearPrecision
void RotationWrite(xercesc::DOMElement *element, const G4String &name, const G4ThreeVector &rot)
static const G4double kAngularPrecision
void FirstpositionWrite(xercesc::DOMElement *element, const G4String &name, const G4ThreeVector &pos)
void PositionWrite(xercesc::DOMElement *element, const G4String &name, const G4ThreeVector &pos)
G4ThreeVector GetAngles(const G4RotationMatrix &)
virtual void AddSolid(const G4VSolid *const)
static const G4int maxTransforms
Definition: xmlparse.cc:187

References AddSolid(), FatalException, G4GDMLWriteDefine::FirstpositionWrite(), G4GDMLWriteDefine::FirstrotationWrite(), G4Exception(), G4GDMLWrite::GenerateName(), G4GDMLWriteDefine::GetAngles(), G4VSolid::GetName(), G4GDMLWriteDefine::kAngularPrecision, G4GDMLWriteDefine::kLinearPrecision, maxTransforms, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), pos, G4GDMLWriteDefine::PositionWrite(), G4GDMLWriteDefine::RotationWrite(), CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by AddSolid().

◆ BoxWrite()

void G4GDMLWriteSolids::BoxWrite ( xercesc::DOMElement *  solElement,
const G4Box * const  box 
)
protected

Definition at line 292 of file G4GDMLWriteSolids.cc.

294{
295 const G4String& name = GenerateName(box->GetName(), box);
296
297 xercesc::DOMElement* boxElement = NewElement("box");
298 boxElement->setAttributeNode(NewAttribute("name", name));
299 boxElement->setAttributeNode(
300 NewAttribute("x", 2.0 * box->GetXHalfLength() / mm));
301 boxElement->setAttributeNode(
302 NewAttribute("y", 2.0 * box->GetYHalfLength() / mm));
303 boxElement->setAttributeNode(
304 NewAttribute("z", 2.0 * box->GetZHalfLength() / mm));
305 boxElement->setAttributeNode(NewAttribute("lunit", "mm"));
306 solElement->appendChild(boxElement);
307}
static constexpr double mm
Definition: G4SIunits.hh:95
G4double GetYHalfLength() const
G4double GetZHalfLength() const
G4double GetXHalfLength() const

References G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4Box::GetXHalfLength(), G4Box::GetYHalfLength(), G4Box::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ ConeWrite()

void G4GDMLWriteSolids::ConeWrite ( xercesc::DOMElement *  solElement,
const G4Cons * const  cone 
)
protected

Definition at line 310 of file G4GDMLWriteSolids.cc.

312{
313 const G4String& name = GenerateName(cone->GetName(), cone);
314
315 xercesc::DOMElement* coneElement = NewElement("cone");
316 coneElement->setAttributeNode(NewAttribute("name", name));
317 coneElement->setAttributeNode(
318 NewAttribute("rmin1", cone->GetInnerRadiusMinusZ() / mm));
319 coneElement->setAttributeNode(
320 NewAttribute("rmax1", cone->GetOuterRadiusMinusZ() / mm));
321 coneElement->setAttributeNode(
322 NewAttribute("rmin2", cone->GetInnerRadiusPlusZ() / mm));
323 coneElement->setAttributeNode(
324 NewAttribute("rmax2", cone->GetOuterRadiusPlusZ() / mm));
325 coneElement->setAttributeNode(
326 NewAttribute("z", 2.0 * cone->GetZHalfLength() / mm));
327 coneElement->setAttributeNode(
328 NewAttribute("startphi", cone->GetStartPhiAngle() / degree));
329 coneElement->setAttributeNode(
330 NewAttribute("deltaphi", cone->GetDeltaPhiAngle() / degree));
331 coneElement->setAttributeNode(NewAttribute("aunit", "deg"));
332 coneElement->setAttributeNode(NewAttribute("lunit", "mm"));
333 solElement->appendChild(coneElement);
334}
static constexpr double degree
Definition: G4SIunits.hh:124
G4double GetOuterRadiusPlusZ() const
G4double GetStartPhiAngle() const
G4double GetDeltaPhiAngle() const
G4double GetInnerRadiusMinusZ() const
G4double GetInnerRadiusPlusZ() const
G4double GetOuterRadiusMinusZ() const
G4double GetZHalfLength() const

References degree, G4GDMLWrite::GenerateName(), G4Cons::GetDeltaPhiAngle(), G4Cons::GetInnerRadiusMinusZ(), G4Cons::GetInnerRadiusPlusZ(), G4VSolid::GetName(), G4Cons::GetOuterRadiusMinusZ(), G4Cons::GetOuterRadiusPlusZ(), G4Cons::GetStartPhiAngle(), G4Cons::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ CutTubeWrite()

void G4GDMLWriteSolids::CutTubeWrite ( xercesc::DOMElement *  solElement,
const G4CutTubs * const  cuttube 
)
protected

Definition at line 894 of file G4GDMLWriteSolids.cc.

896{
897 const G4String& name = GenerateName(cuttube->GetName(), cuttube);
898
899 xercesc::DOMElement* cuttubeElement = NewElement("cutTube");
900 cuttubeElement->setAttributeNode(NewAttribute("name", name));
901 cuttubeElement->setAttributeNode(
902 NewAttribute("rmin", cuttube->GetInnerRadius() / mm));
903 cuttubeElement->setAttributeNode(
904 NewAttribute("rmax", cuttube->GetOuterRadius() / mm));
905 cuttubeElement->setAttributeNode(
906 NewAttribute("z", 2.0 * cuttube->GetZHalfLength() / mm));
907 cuttubeElement->setAttributeNode(
908 NewAttribute("startphi", cuttube->GetStartPhiAngle() / degree));
909 cuttubeElement->setAttributeNode(
910 NewAttribute("deltaphi", cuttube->GetDeltaPhiAngle() / degree));
911 cuttubeElement->setAttributeNode(
912 NewAttribute("lowX", cuttube->GetLowNorm().getX() / mm));
913 cuttubeElement->setAttributeNode(
914 NewAttribute("lowY", cuttube->GetLowNorm().getY() / mm));
915 cuttubeElement->setAttributeNode(
916 NewAttribute("lowZ", cuttube->GetLowNorm().getZ() / mm));
917 cuttubeElement->setAttributeNode(
918 NewAttribute("highX", cuttube->GetHighNorm().getX() / mm));
919 cuttubeElement->setAttributeNode(
920 NewAttribute("highY", cuttube->GetHighNorm().getY() / mm));
921 cuttubeElement->setAttributeNode(
922 NewAttribute("highZ", cuttube->GetHighNorm().getZ() / mm));
923 cuttubeElement->setAttributeNode(NewAttribute("aunit", "deg"));
924 cuttubeElement->setAttributeNode(NewAttribute("lunit", "mm"));
925 solElement->appendChild(cuttubeElement);
926}
double getZ() const
double getX() const
double getY() const
G4ThreeVector GetHighNorm() const
G4double GetStartPhiAngle() const
G4double GetZHalfLength() const
G4double GetInnerRadius() const
G4ThreeVector GetLowNorm() const
G4double GetDeltaPhiAngle() const
G4double GetOuterRadius() const

References degree, G4GDMLWrite::GenerateName(), G4CutTubs::GetDeltaPhiAngle(), G4CutTubs::GetHighNorm(), G4CutTubs::GetInnerRadius(), G4CutTubs::GetLowNorm(), G4VSolid::GetName(), G4CutTubs::GetOuterRadius(), G4CutTubs::GetStartPhiAngle(), CLHEP::Hep3Vector::getX(), CLHEP::Hep3Vector::getY(), CLHEP::Hep3Vector::getZ(), G4CutTubs::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ DefineWrite()

void G4GDMLWriteDefine::DefineWrite ( xercesc::DOMElement *  element)
virtualinherited

Implements G4GDMLWrite.

Definition at line 137 of file G4GDMLWriteDefine.cc.

138{
139#ifdef G4VERBOSE
140 G4cout << "G4GDML: Writing definitions..." << G4endl;
141#endif
142 defineElement = NewElement("define");
143 element->appendChild(defineElement);
144}

References G4GDMLWriteDefine::defineElement, G4cout, G4endl, and G4GDMLWrite::NewElement().

◆ DepthMap()

G4GDMLWrite::DepthMapType & G4GDMLWrite::DepthMap ( )
protectedinherited

Definition at line 73 of file G4GDMLWrite.cc.

74{
75 static DepthMapType instance;
76 return instance;
77}
std::map< G4int, G4int > DepthMapType
Definition: G4GDMLWrite.hh:55

Referenced by G4GDMLWrite::AddModule(), and G4GDMLWrite::Modularize().

◆ DWrite()

void G4GDMLWriteMaterials::DWrite ( xercesc::DOMElement *  element,
const G4double d 
)
protectedinherited

Definition at line 63 of file G4GDMLWriteMaterials.cc.

65{
66 xercesc::DOMElement* DElement = NewElement("D");
67 DElement->setAttributeNode(NewAttribute("unit", "g/cm3"));
68 DElement->setAttributeNode(NewAttribute("value", d * cm3 / g));
69 element->appendChild(DElement);
70}
static constexpr double cm3
Definition: G4SIunits.hh:101

References cm3, g, G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by G4GDMLWriteMaterials::MaterialWrite().

◆ ElconeWrite()

void G4GDMLWriteSolids::ElconeWrite ( xercesc::DOMElement *  solElement,
const G4EllipticalCone * const  elcone 
)
protected

Definition at line 337 of file G4GDMLWriteSolids.cc.

339{
340 const G4String& name = GenerateName(elcone->GetName(), elcone);
341
342 xercesc::DOMElement* elconeElement = NewElement("elcone");
343 elconeElement->setAttributeNode(NewAttribute("name", name));
344 elconeElement->setAttributeNode(
345 NewAttribute("dx", elcone->GetSemiAxisX() / mm));
346 elconeElement->setAttributeNode(
347 NewAttribute("dy", elcone->GetSemiAxisY() / mm));
348 elconeElement->setAttributeNode(NewAttribute("zmax", elcone->GetZMax() / mm));
349 elconeElement->setAttributeNode(
350 NewAttribute("zcut", elcone->GetZTopCut() / mm));
351 elconeElement->setAttributeNode(NewAttribute("lunit", "mm"));
352 solElement->appendChild(elconeElement);
353}
G4double GetSemiAxisX() const
G4double GetSemiAxisY() const
G4double GetZMax() const
G4double GetZTopCut() const

References G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4EllipticalCone::GetSemiAxisX(), G4EllipticalCone::GetSemiAxisY(), G4EllipticalCone::GetZMax(), G4EllipticalCone::GetZTopCut(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ ElementWrite()

void G4GDMLWriteMaterials::ElementWrite ( const G4Element * const  elementPtr)
protectedinherited

Definition at line 116 of file G4GDMLWriteMaterials.cc.

117{
118 const G4String name = GenerateName(elementPtr->GetName(), elementPtr);
119
120 xercesc::DOMElement* elementElement = NewElement("element");
121 elementElement->setAttributeNode(NewAttribute("name", name));
122
123 const std::size_t NumberOfIsotopes = elementPtr->GetNumberOfIsotopes();
124
125 if(NumberOfIsotopes > 0)
126 {
127 const G4double* RelativeAbundanceVector
128 = elementPtr->GetRelativeAbundanceVector();
129 for(std::size_t i = 0; i < NumberOfIsotopes; ++i)
130 {
131 G4String fractionref = GenerateName(elementPtr->GetIsotope(i)->GetName(),
132 elementPtr->GetIsotope(i));
133 xercesc::DOMElement* fractionElement = NewElement("fraction");
134 fractionElement->setAttributeNode(
135 NewAttribute("n", RelativeAbundanceVector[i]));
136 fractionElement->setAttributeNode(NewAttribute("ref", fractionref));
137 elementElement->appendChild(fractionElement);
138 AddIsotope(elementPtr->GetIsotope(i));
139 }
140 }
141 else
142 {
143 elementElement->setAttributeNode(NewAttribute("Z", elementPtr->GetZ()));
144 AtomWrite(elementElement, elementPtr->GetA());
145 }
146
147 materialsElement->appendChild(elementElement);
148 // Append the element AFTER all the possible components are appended!
149}
double G4double
Definition: G4Types.hh:83
G4double * GetRelativeAbundanceVector() const
Definition: G4Element.hh:167
G4double GetZ() const
Definition: G4Element.hh:131
const G4Isotope * GetIsotope(G4int iso) const
Definition: G4Element.hh:170
G4double GetA() const
Definition: G4Element.hh:139
size_t GetNumberOfIsotopes() const
Definition: G4Element.hh:159
const G4String & GetName() const
Definition: G4Element.hh:127
void AtomWrite(xercesc::DOMElement *, const G4double &)
void AddIsotope(const G4Isotope *const)
xercesc::DOMElement * materialsElement
const G4String & GetName() const
Definition: G4Isotope.hh:87

References G4GDMLWriteMaterials::AddIsotope(), G4GDMLWriteMaterials::AtomWrite(), G4GDMLWrite::GenerateName(), G4Element::GetA(), G4Element::GetIsotope(), G4Element::GetName(), G4Isotope::GetName(), G4Element::GetNumberOfIsotopes(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZ(), G4GDMLWriteMaterials::materialsElement, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by G4GDMLWriteMaterials::AddElement().

◆ EllipsoidWrite()

void G4GDMLWriteSolids::EllipsoidWrite ( xercesc::DOMElement *  solElement,
const G4Ellipsoid * const  ellipsoid 
)
protected

Definition at line 356 of file G4GDMLWriteSolids.cc.

358{
359 const G4String& name = GenerateName(ellipsoid->GetName(), ellipsoid);
360
361 xercesc::DOMElement* ellipsoidElement = NewElement("ellipsoid");
362 ellipsoidElement->setAttributeNode(NewAttribute("name", name));
363 ellipsoidElement->setAttributeNode(
364 NewAttribute("ax", ellipsoid->GetSemiAxisMax(0) / mm));
365 ellipsoidElement->setAttributeNode(
366 NewAttribute("by", ellipsoid->GetSemiAxisMax(1) / mm));
367 ellipsoidElement->setAttributeNode(
368 NewAttribute("cz", ellipsoid->GetSemiAxisMax(2) / mm));
369 ellipsoidElement->setAttributeNode(
370 NewAttribute("zcut1", ellipsoid->GetZBottomCut() / mm));
371 ellipsoidElement->setAttributeNode(
372 NewAttribute("zcut2", ellipsoid->GetZTopCut() / mm));
373 ellipsoidElement->setAttributeNode(NewAttribute("lunit", "mm"));
374 solElement->appendChild(ellipsoidElement);
375}
G4double GetSemiAxisMax(G4int i) const
G4double GetZTopCut() const
G4double GetZBottomCut() const

References G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4Ellipsoid::GetSemiAxisMax(), G4Ellipsoid::GetZBottomCut(), G4Ellipsoid::GetZTopCut(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ EltubeWrite()

void G4GDMLWriteSolids::EltubeWrite ( xercesc::DOMElement *  solElement,
const G4EllipticalTube * const  eltube 
)
protected

Definition at line 378 of file G4GDMLWriteSolids.cc.

380{
381 const G4String& name = GenerateName(eltube->GetName(), eltube);
382
383 xercesc::DOMElement* eltubeElement = NewElement("eltube");
384 eltubeElement->setAttributeNode(NewAttribute("name", name));
385 eltubeElement->setAttributeNode(NewAttribute("dx", eltube->GetDx() / mm));
386 eltubeElement->setAttributeNode(NewAttribute("dy", eltube->GetDy() / mm));
387 eltubeElement->setAttributeNode(NewAttribute("dz", eltube->GetDz() / mm));
388 eltubeElement->setAttributeNode(NewAttribute("lunit", "mm"));
389 solElement->appendChild(eltubeElement);
390}
G4double GetDy() const
G4double GetDx() const
G4double GetDz() const

References G4GDMLWrite::GenerateName(), G4EllipticalTube::GetDx(), G4EllipticalTube::GetDy(), G4EllipticalTube::GetDz(), G4VSolid::GetName(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ ExtensionWrite()

void G4GDMLWrite::ExtensionWrite ( xercesc::DOMElement *  )
virtualinherited

Definition at line 88 of file G4GDMLWrite.cc.

89{
90 // Empty implementation. To be overwritten by user for specific extensions
91}

Referenced by G4GDMLWrite::Write().

◆ FileExists()

G4bool G4GDMLWrite::FileExists ( const G4String fname) const
protectedinherited

Definition at line 53 of file G4GDMLWrite.cc.

54{
55 struct stat FileInfo;
56 return (stat(fname.c_str(), &FileInfo) == 0);
57}

References test::fname.

Referenced by G4GDMLWrite::Write().

◆ FirstpositionWrite()

void G4GDMLWriteDefine::FirstpositionWrite ( xercesc::DOMElement *  element,
const G4String name,
const G4ThreeVector pos 
)
inlineinherited

Definition at line 68 of file G4GDMLWriteDefine.hh.

70 {
71 Position_vectorWrite(element, "firstposition", name, pos);
72 }

References G4InuclParticleNames::name(), pos, and G4GDMLWriteDefine::Position_vectorWrite().

Referenced by BooleanWrite().

◆ FirstrotationWrite()

void G4GDMLWriteDefine::FirstrotationWrite ( xercesc::DOMElement *  element,
const G4String name,
const G4ThreeVector rot 
)
inlineinherited

Definition at line 63 of file G4GDMLWriteDefine.hh.

65 {
66 Rotation_vectorWrite(element, "firstrotation", name, rot);
67 }
void Rotation_vectorWrite(xercesc::DOMElement *, const G4String &, const G4String &, const G4ThreeVector &)

References G4InuclParticleNames::name(), and G4GDMLWriteDefine::Rotation_vectorWrite().

Referenced by BooleanWrite().

◆ GenerateName()

G4String G4GDMLWrite::GenerateName ( const G4String name,
const void * const  ptr 
)
inherited

Definition at line 132 of file G4GDMLWrite.cc.

133{
134 G4String nameOut;
135 std::stringstream stream;
136 stream << name;
138 {
139 stream << ptr;
140 };
141
142 nameOut = G4String(stream.str());
143 std::vector<char> toremove = { ' ', '/', ':', '#', '+' };
144 for(auto c : toremove)
145 {
146 if(G4StrUtil::contains(nameOut, c))
147 {
148 std::replace(nameOut.begin(), nameOut.end(), c, '_');
149 }
150 }
151 return nameOut;
152}
static G4bool addPointerToName
Definition: G4GDMLWrite.hh:131
G4bool contains(const G4String &str, std::string_view ss)
Check if a string contains a given substring.

References G4GDMLWrite::addPointerToName, G4StrUtil::contains(), and G4InuclParticleNames::name().

Referenced by G4GDMLWrite::AddModule(), G4GDMLWriteStructure::AssemblyWrite(), BooleanWrite(), G4GDMLWriteStructure::BorderSurfaceCache(), BoxWrite(), ConeWrite(), CutTubeWrite(), G4GDMLWriteStructure::DivisionvolWrite(), ElconeWrite(), G4GDMLWriteMaterials::ElementWrite(), EllipsoidWrite(), EltubeWrite(), G4GDMLParser::ExportRegions(), GenericPolyconeWrite(), GenTrapWrite(), HypeWrite(), G4GDMLWriteMaterials::IsotopeWrite(), G4GDMLWriteMaterials::MaterialWrite(), MultiUnionWrite(), OpticalSurfaceWrite(), OrbWrite(), ParaboloidWrite(), G4GDMLWriteParamvol::ParametersWrite(), G4GDMLWriteParamvol::ParamvolAlgorithmWrite(), G4GDMLWriteParamvol::ParamvolWrite(), ParaWrite(), G4GDMLWriteStructure::PhysvolWrite(), PolyconeWrite(), PolyhedraWrite(), G4GDMLWriteMaterials::PropertyConstWrite(), G4GDMLWriteMaterials::PropertyVectorWrite(), G4GDMLWriteMaterials::PropertyWrite(), PropertyWrite(), G4GDMLWriteStructure::ReplicavolWrite(), ScaledWrite(), G4GDMLWriteSetup::SetupWrite(), G4GDMLWriteStructure::SkinSurfaceCache(), SphereWrite(), TessellatedWrite(), TetWrite(), TorusWrite(), TrapWrite(), G4GDMLWriteStructure::TraverseVolumeTree(), TrdWrite(), TubeWrite(), TwistedboxWrite(), TwistedtrapWrite(), TwistedtrdWrite(), TwistedtubsWrite(), and XtruWrite().

◆ GenericPolyconeWrite()

void G4GDMLWriteSolids::GenericPolyconeWrite ( xercesc::DOMElement *  solElement,
const G4GenericPolycone * const  polycone 
)
protected

Definition at line 547 of file G4GDMLWriteSolids.cc.

549{
550 const G4String& name = GenerateName(polycone->GetName(), polycone);
551 xercesc::DOMElement* polyconeElement = NewElement("genericPolycone");
552 const G4double startPhi = polycone->GetStartPhi();
553 polyconeElement->setAttributeNode(NewAttribute("name", name));
554 polyconeElement->setAttributeNode(
555 NewAttribute("startphi", startPhi / degree));
556 polyconeElement->setAttributeNode(
557 NewAttribute("deltaphi", (polycone->GetEndPhi() - startPhi) / degree));
558 polyconeElement->setAttributeNode(NewAttribute("aunit", "deg"));
559 polyconeElement->setAttributeNode(NewAttribute("lunit", "mm"));
560 solElement->appendChild(polyconeElement);
561
562 const std::size_t num_rzpoints = polycone->GetNumRZCorner();
563 for(std::size_t i = 0; i < num_rzpoints; ++i)
564 {
565 const G4double r_point = polycone->GetCorner(i).r;
566 const G4double z_point = polycone->GetCorner(i).z;
567 RZPointWrite(polyconeElement, r_point, z_point);
568 }
569}
void RZPointWrite(xercesc::DOMElement *, const G4double &, const G4double &)
G4double GetStartPhi() const
G4double GetEndPhi() const
G4int GetNumRZCorner() const
G4PolyconeSideRZ GetCorner(G4int index) const

References degree, G4GDMLWrite::GenerateName(), G4GenericPolycone::GetCorner(), G4GenericPolycone::GetEndPhi(), G4VSolid::GetName(), G4GenericPolycone::GetNumRZCorner(), G4GenericPolycone::GetStartPhi(), G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), G4PolyconeSideRZ::r, RZPointWrite(), and G4PolyconeSideRZ::z.

Referenced by AddSolid().

◆ GenTrapWrite()

void G4GDMLWriteSolids::GenTrapWrite ( xercesc::DOMElement *  solElement,
const G4GenericTrap * const  gtrap 
)
protected

Definition at line 780 of file G4GDMLWriteSolids.cc.

782{
783 const G4String& name = GenerateName(gtrap->GetName(), gtrap);
784
785 std::vector<G4TwoVector> vertices = gtrap->GetVertices();
786
787 xercesc::DOMElement* gtrapElement = NewElement("arb8");
788 gtrapElement->setAttributeNode(NewAttribute("name", name));
789 gtrapElement->setAttributeNode(
790 NewAttribute("dz", gtrap->GetZHalfLength() / mm));
791 gtrapElement->setAttributeNode(NewAttribute("v1x", vertices[0].x()));
792 gtrapElement->setAttributeNode(NewAttribute("v1y", vertices[0].y()));
793 gtrapElement->setAttributeNode(NewAttribute("v2x", vertices[1].x()));
794 gtrapElement->setAttributeNode(NewAttribute("v2y", vertices[1].y()));
795 gtrapElement->setAttributeNode(NewAttribute("v3x", vertices[2].x()));
796 gtrapElement->setAttributeNode(NewAttribute("v3y", vertices[2].y()));
797 gtrapElement->setAttributeNode(NewAttribute("v4x", vertices[3].x()));
798 gtrapElement->setAttributeNode(NewAttribute("v4y", vertices[3].y()));
799 gtrapElement->setAttributeNode(NewAttribute("v5x", vertices[4].x()));
800 gtrapElement->setAttributeNode(NewAttribute("v5y", vertices[4].y()));
801 gtrapElement->setAttributeNode(NewAttribute("v6x", vertices[5].x()));
802 gtrapElement->setAttributeNode(NewAttribute("v6y", vertices[5].y()));
803 gtrapElement->setAttributeNode(NewAttribute("v7x", vertices[6].x()));
804 gtrapElement->setAttributeNode(NewAttribute("v7y", vertices[6].y()));
805 gtrapElement->setAttributeNode(NewAttribute("v8x", vertices[7].x()));
806 gtrapElement->setAttributeNode(NewAttribute("v8y", vertices[7].y()));
807 gtrapElement->setAttributeNode(NewAttribute("lunit", "mm"));
808 solElement->appendChild(gtrapElement);
809}
G4double GetZHalfLength() const
const std::vector< G4TwoVector > & GetVertices() const

References G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4GenericTrap::GetVertices(), G4GenericTrap::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ GetAngles()

G4ThreeVector G4GDMLWriteDefine::GetAngles ( const G4RotationMatrix mtx)
inherited

Definition at line 50 of file G4GDMLWriteDefine.cc.

51{
52 G4double x, y, z;
53 G4RotationMatrix mat = mtx;
54 mat.rectify(); // Rectify matrix from possible roundoff errors
55
56 // Direction of rotation given by left-hand rule; clockwise rotation
57
58 static const G4double kMatrixPrecision = 10E-10;
59 const G4double cosb = std::sqrt(mtx.xx() * mtx.xx() + mtx.yx() * mtx.yx());
60
61 if(cosb > kMatrixPrecision)
62 {
63 x = std::atan2(mtx.zy(), mtx.zz());
64 y = std::atan2(-mtx.zx(), cosb);
65 z = std::atan2(mtx.yx(), mtx.xx());
66 }
67 else
68 {
69 x = std::atan2(-mtx.yz(), mtx.yy());
70 y = std::atan2(-mtx.zx(), cosb);
71 z = 0.0;
72 }
73
74 return G4ThreeVector(x, y, z);
75}
CLHEP::Hep3Vector G4ThreeVector
double zz() const
double yz() const
double zx() const
double yx() const
double zy() const
double xx() const
double yy() const

References CLHEP::HepRotation::rectify(), CLHEP::HepRotation::xx(), CLHEP::HepRotation::yx(), CLHEP::HepRotation::yy(), CLHEP::HepRotation::yz(), CLHEP::HepRotation::zx(), CLHEP::HepRotation::zy(), and CLHEP::HepRotation::zz().

Referenced by G4GDMLWriteStructure::AssemblyWrite(), BooleanWrite(), MultiUnionWrite(), G4GDMLWriteParamvol::ParametersWrite(), G4GDMLWriteStructure::PhysvolWrite(), and G4GDMLWriteStructure::TraverseVolumeTree().

◆ HypeWrite()

void G4GDMLWriteSolids::HypeWrite ( xercesc::DOMElement *  solElement,
const G4Hype * const  hype 
)
protected

Definition at line 438 of file G4GDMLWriteSolids.cc.

440{
441 const G4String& name = GenerateName(hype->GetName(), hype);
442
443 xercesc::DOMElement* hypeElement = NewElement("hype");
444 hypeElement->setAttributeNode(NewAttribute("name", name));
445 hypeElement->setAttributeNode(
446 NewAttribute("rmin", hype->GetInnerRadius() / mm));
447 hypeElement->setAttributeNode(
448 NewAttribute("rmax", hype->GetOuterRadius() / mm));
449 hypeElement->setAttributeNode(
450 NewAttribute("inst", hype->GetInnerStereo() / degree));
451 hypeElement->setAttributeNode(
452 NewAttribute("outst", hype->GetOuterStereo() / degree));
453 hypeElement->setAttributeNode(
454 NewAttribute("z", 2.0 * hype->GetZHalfLength() / mm));
455 hypeElement->setAttributeNode(NewAttribute("aunit", "deg"));
456 hypeElement->setAttributeNode(NewAttribute("lunit", "mm"));
457 solElement->appendChild(hypeElement);
458}
G4double GetInnerStereo() const
G4double GetZHalfLength() const
G4double GetOuterStereo() const
G4double GetOuterRadius() const
G4double GetInnerRadius() const

References degree, G4GDMLWrite::GenerateName(), G4Hype::GetInnerRadius(), G4Hype::GetInnerStereo(), G4VSolid::GetName(), G4Hype::GetOuterRadius(), G4Hype::GetOuterStereo(), G4Hype::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ IsotopeWrite()

void G4GDMLWriteMaterials::IsotopeWrite ( const G4Isotope * const  isotopePtr)
protectedinherited

Definition at line 103 of file G4GDMLWriteMaterials.cc.

104{
105 const G4String name = GenerateName(isotopePtr->GetName(), isotopePtr);
106
107 xercesc::DOMElement* isotopeElement = NewElement("isotope");
108 isotopeElement->setAttributeNode(NewAttribute("name", name));
109 isotopeElement->setAttributeNode(NewAttribute("N", isotopePtr->GetN()));
110 isotopeElement->setAttributeNode(NewAttribute("Z", isotopePtr->GetZ()));
111 materialsElement->appendChild(isotopeElement);
112 AtomWrite(isotopeElement, isotopePtr->GetA());
113}
G4int GetZ() const
Definition: G4Isotope.hh:90
G4int GetN() const
Definition: G4Isotope.hh:93
G4double GetA() const
Definition: G4Isotope.hh:96

References G4GDMLWriteMaterials::AtomWrite(), G4GDMLWrite::GenerateName(), G4Isotope::GetA(), G4Isotope::GetN(), G4Isotope::GetName(), G4Isotope::GetZ(), G4GDMLWriteMaterials::materialsElement, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by G4GDMLWriteMaterials::AddIsotope().

◆ MaterialsWrite()

void G4GDMLWriteMaterials::MaterialsWrite ( xercesc::DOMElement *  element)
virtualinherited

Implements G4GDMLWrite.

Definition at line 321 of file G4GDMLWriteMaterials.cc.

322{
323#ifdef G4VERBOSE
324 G4cout << "G4GDML: Writing materials..." << G4endl;
325#endif
326 materialsElement = NewElement("materials");
327 element->appendChild(materialsElement);
328
329 isotopeList.clear();
330 elementList.clear();
331 materialList.clear();
332 propertyList.clear();
333}
std::vector< const G4PhysicsFreeVector * > propertyList

References G4GDMLWriteMaterials::elementList, G4cout, G4endl, G4GDMLWriteMaterials::isotopeList, G4GDMLWriteMaterials::materialList, G4GDMLWriteMaterials::materialsElement, G4GDMLWrite::NewElement(), and G4GDMLWriteMaterials::propertyList.

◆ MaterialWrite()

void G4GDMLWriteMaterials::MaterialWrite ( const G4Material * const  materialPtr)
protectedinherited

Definition at line 152 of file G4GDMLWriteMaterials.cc.

153{
154 G4String state_str("undefined");
155 const G4State state = materialPtr->GetState();
156 if(state == kStateSolid)
157 {
158 state_str = "solid";
159 }
160 else if(state == kStateLiquid)
161 {
162 state_str = "liquid";
163 }
164 else if(state == kStateGas)
165 {
166 state_str = "gas";
167 }
168
169 const G4String name = GenerateName(materialPtr->GetName(), materialPtr);
170
171 xercesc::DOMElement* materialElement = NewElement("material");
172 materialElement->setAttributeNode(NewAttribute("name", name));
173 materialElement->setAttributeNode(NewAttribute("state", state_str));
174
175 // Write any property attached to the material...
176 //
177 if(materialPtr->GetMaterialPropertiesTable())
178 {
179 PropertyWrite(materialElement, materialPtr);
180 }
181
182 if(materialPtr->GetTemperature() != STP_Temperature)
183 {
184 TWrite(materialElement, materialPtr->GetTemperature());
185 }
186
187 if(materialPtr->GetPressure() != STP_Pressure)
188 {
189 PWrite(materialElement, materialPtr->GetPressure());
190 }
191
192 // Write Ionisation potential (mean excitation energy)
193 MEEWrite(materialElement,
194 materialPtr->GetIonisation()->GetMeanExcitationEnergy());
195
196 DWrite(materialElement, materialPtr->GetDensity());
197
198 const std::size_t NumberOfElements = materialPtr->GetNumberOfElements();
199
200 if((NumberOfElements > 1) ||
201 (materialPtr->GetElement(0) != nullptr &&
202 materialPtr->GetElement(0)->GetNumberOfIsotopes() > 1))
203 {
204 const G4double* MassFractionVector = materialPtr->GetFractionVector();
205
206 for(std::size_t i = 0; i < NumberOfElements; ++i)
207 {
208 const G4String fractionref = GenerateName(
209 materialPtr->GetElement(i)->GetName(), materialPtr->GetElement(i));
210 xercesc::DOMElement* fractionElement = NewElement("fraction");
211 fractionElement->setAttributeNode(
212 NewAttribute("n", MassFractionVector[i]));
213 fractionElement->setAttributeNode(NewAttribute("ref", fractionref));
214 materialElement->appendChild(fractionElement);
215 AddElement(materialPtr->GetElement(i));
216 }
217 }
218 else
219 {
220 materialElement->setAttributeNode(NewAttribute("Z", materialPtr->GetZ()));
221 AtomWrite(materialElement, materialPtr->GetA());
222 }
223
224 // Append the material AFTER all the possible components are appended!
225 //
226 materialsElement->appendChild(materialElement);
227}
G4State
Definition: G4Material.hh:111
@ kStateSolid
Definition: G4Material.hh:111
@ kStateLiquid
Definition: G4Material.hh:111
@ kStateGas
Definition: G4Material.hh:111
void TWrite(xercesc::DOMElement *, const G4double &)
void MEEWrite(xercesc::DOMElement *, const G4double &)
void DWrite(xercesc::DOMElement *, const G4double &)
void PWrite(xercesc::DOMElement *, const G4double &)
void AddElement(const G4Element *const)
void PropertyWrite(xercesc::DOMElement *, const G4Material *const)
G4double GetMeanExcitationEnergy() const
G4double GetPressure() const
Definition: G4Material.hh:179
G4double GetDensity() const
Definition: G4Material.hh:176
G4MaterialPropertiesTable * GetMaterialPropertiesTable() const
Definition: G4Material.hh:252
G4State GetState() const
Definition: G4Material.hh:177
G4double GetTemperature() const
Definition: G4Material.hh:178
const G4Element * GetElement(G4int iel) const
Definition: G4Material.hh:198
G4double GetZ() const
Definition: G4Material.cc:736
const G4double * GetFractionVector() const
Definition: G4Material.hh:190
G4IonisParamMat * GetIonisation() const
Definition: G4Material.hh:222
size_t GetNumberOfElements() const
Definition: G4Material.hh:182
G4double GetA() const
Definition: G4Material.cc:750
const G4String & GetName() const
Definition: G4Material.hh:173
float STP_Temperature
Definition: hepunit.py:301
int STP_Pressure
Definition: hepunit.py:302

References G4GDMLWriteMaterials::AddElement(), G4GDMLWriteMaterials::AtomWrite(), G4GDMLWriteMaterials::DWrite(), G4GDMLWrite::GenerateName(), G4Material::GetA(), G4Material::GetDensity(), G4Material::GetElement(), G4Material::GetFractionVector(), G4Material::GetIonisation(), G4Material::GetMaterialPropertiesTable(), G4IonisParamMat::GetMeanExcitationEnergy(), G4Element::GetName(), G4Material::GetName(), G4Material::GetNumberOfElements(), G4Element::GetNumberOfIsotopes(), G4Material::GetPressure(), G4Material::GetState(), G4Material::GetTemperature(), G4Material::GetZ(), kStateGas, kStateLiquid, kStateSolid, G4GDMLWriteMaterials::materialsElement, G4GDMLWriteMaterials::MEEWrite(), G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), G4GDMLWriteMaterials::PropertyWrite(), G4GDMLWriteMaterials::PWrite(), source.hepunit::STP_Pressure, source.hepunit::STP_Temperature, and G4GDMLWriteMaterials::TWrite().

Referenced by G4GDMLWriteMaterials::AddMaterial().

◆ MEEWrite()

void G4GDMLWriteMaterials::MEEWrite ( xercesc::DOMElement *  element,
const G4double MEE 
)
protectedinherited

Definition at line 93 of file G4GDMLWriteMaterials.cc.

95{
96 xercesc::DOMElement* PElement = NewElement("MEE");
97 PElement->setAttributeNode(NewAttribute("unit", "eV"));
98 PElement->setAttributeNode(NewAttribute("value", MEE / electronvolt));
99 element->appendChild(PElement);
100}
static constexpr double electronvolt
Definition: G4SIunits.hh:191

References electronvolt, G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by G4GDMLWriteMaterials::MaterialWrite().

◆ Modularize()

G4String G4GDMLWrite::Modularize ( const G4VPhysicalVolume *const  topvol,
const G4int  depth 
)
protectedinherited

Definition at line 377 of file G4GDMLWrite.cc.

379{
380 if(PvolumeMap().find(physvol) != PvolumeMap().cend())
381 {
382 return PvolumeMap()[physvol]; // Modularize via physvol
383 }
384
385 if(DepthMap().find(depth) != DepthMap().cend()) // Modularize via depth
386 {
387 std::stringstream stream;
388 stream << "depth" << depth << "_module" << DepthMap()[depth] << ".gdml";
389 DepthMap()[depth]++; // There can be more modules at this depth!
390 return G4String(stream.str());
391 }
392
393 return G4String(""); // Empty string for module name = no modularization
394 // was requested at that level/physvol!
395}

References G4GDMLWrite::DepthMap(), and G4GDMLWrite::PvolumeMap().

Referenced by G4GDMLWriteStructure::TraverseVolumeTree().

◆ MultiUnionWrite()

void G4GDMLWriteSolids::MultiUnionWrite ( xercesc::DOMElement *  solElement,
const G4MultiUnion * const  munionSolid 
)
protected

Definition at line 82 of file G4GDMLWriteSolids.cc.

84{
85 G4int numSolids = munionSolid->GetNumberOfSolids();
86 G4String tag("multiUnion");
87
88 G4VSolid* solid;
90
91 const G4String& name = GenerateName(munionSolid->GetName(), munionSolid);
92 xercesc::DOMElement* multiUnionElement = NewElement(tag);
93 multiUnionElement->setAttributeNode(NewAttribute("name", name));
94
95 for(G4int i = 0; i < numSolids; ++i)
96 {
97 solid = munionSolid->GetSolid(i);
98 transform = munionSolid->GetTransformation(i);
99
100 HepGeom::Rotate3D rot3d;
102 HepGeom::Scale3D scale;
103 transform.getDecomposition(scale, rot3d, transl);
104
106 G4RotationMatrix rotm(CLHEP::HepRep3x3(rot3d.xx(), rot3d.xy(), rot3d.xz(),
107 rot3d.yx(), rot3d.yy(), rot3d.yz(),
108 rot3d.zx(), rot3d.zy(), rot3d.zz()));
109 G4ThreeVector rot = GetAngles(rotm);
110
111 AddSolid(solid);
112 const G4String& solidref = GenerateName(solid->GetName(), solid);
113 std::ostringstream os;
114 os << i + 1;
115 const G4String& nodeName = "Node-" + G4String(os.str());
116 xercesc::DOMElement* solidElement = NewElement("solid");
117 solidElement->setAttributeNode(NewAttribute("ref", solidref));
118 xercesc::DOMElement* multiUnionNodeElement = NewElement("multiUnionNode");
119 multiUnionNodeElement->setAttributeNode(NewAttribute("name", name+"_"+nodeName));
120 multiUnionNodeElement->appendChild(solidElement); // Append solid to node
121 if((std::fabs(pos.x()) > kLinearPrecision) ||
122 (std::fabs(pos.y()) > kLinearPrecision) ||
123 (std::fabs(pos.z()) > kLinearPrecision))
124 {
125 PositionWrite(multiUnionNodeElement,name+"_"+nodeName+"_pos",pos);
126 }
127 if((std::fabs(rot.x()) > kAngularPrecision) ||
128 (std::fabs(rot.y()) > kAngularPrecision) ||
129 (std::fabs(rot.z()) > kAngularPrecision))
130 {
131 RotationWrite(multiUnionNodeElement,name+"_"+nodeName+"_rot",rot);
132 }
133 multiUnionElement->appendChild(multiUnionNodeElement); // Append node
134 }
135
136 solElement->appendChild(multiUnionElement);
137 // Add the multiUnion solid AFTER the constituent nodes!
138}
const G4Transform3D & GetTransformation(G4int index) const
G4int GetNumberOfSolids() const
G4VSolid * GetSolid(G4int index) const
double zz() const
Definition: Transform3D.h:281
double yz() const
Definition: Transform3D.h:272
double xy() const
Definition: Transform3D.h:260
CLHEP::Hep3Vector getTranslation() const
double zx() const
Definition: Transform3D.h:275
double yx() const
Definition: Transform3D.h:266
double zy() const
Definition: Transform3D.h:278
double xx() const
Definition: Transform3D.h:257
double yy() const
Definition: Transform3D.h:269
double xz() const
Definition: Transform3D.h:263
G4bool transform(G4String &input, const G4String &type)

References AddSolid(), G4GDMLWrite::GenerateName(), G4GDMLWriteDefine::GetAngles(), G4VSolid::GetName(), G4MultiUnion::GetNumberOfSolids(), G4MultiUnion::GetSolid(), G4MultiUnion::GetTransformation(), HepGeom::Transform3D::getTranslation(), G4GDMLWriteDefine::kAngularPrecision, G4GDMLWriteDefine::kLinearPrecision, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), pos, G4GDMLWriteDefine::PositionWrite(), G4GDMLWriteDefine::RotationWrite(), G4coutFormatters::anonymous_namespace{G4coutFormatters.cc}::transform(), CLHEP::Hep3Vector::x(), HepGeom::Transform3D::xx(), HepGeom::Transform3D::xy(), HepGeom::Transform3D::xz(), CLHEP::Hep3Vector::y(), HepGeom::Transform3D::yx(), HepGeom::Transform3D::yy(), HepGeom::Transform3D::yz(), CLHEP::Hep3Vector::z(), HepGeom::Transform3D::zx(), HepGeom::Transform3D::zy(), and HepGeom::Transform3D::zz().

Referenced by AddSolid().

◆ NewAttribute() [1/2]

xercesc::DOMAttr * G4GDMLWrite::NewAttribute ( const G4String name,
const G4double value 
)
protectedinherited

Definition at line 171 of file G4GDMLWrite.cc.

173{
174 XMLCh* tempStr = NULL;
175 tempStr = xercesc::XMLString::transcode(name);
176 xercesc::DOMAttr* att = doc->createAttribute(tempStr);
177 xercesc::XMLString::release(&tempStr);
178
179 std::ostringstream ostream;
180 ostream.precision(15);
181 ostream << value;
182 G4String str = ostream.str();
183
184 tempStr = xercesc::XMLString::transcode(str);
185 att->setValue(tempStr);
186 xercesc::XMLString::release(&tempStr);
187
188 return att;
189}
xercesc::DOMDocument * doc
Definition: G4GDMLWrite.hh:132

References G4GDMLWrite::doc, and G4InuclParticleNames::name().

◆ NewAttribute() [2/2]

xercesc::DOMAttr * G4GDMLWrite::NewAttribute ( const G4String name,
const G4String value 
)
protectedinherited

Definition at line 155 of file G4GDMLWrite.cc.

157{
158 XMLCh* tempStr = NULL;
159 tempStr = xercesc::XMLString::transcode(name);
160 xercesc::DOMAttr* att = doc->createAttribute(tempStr);
161 xercesc::XMLString::release(&tempStr);
162
163 tempStr = xercesc::XMLString::transcode(value);
164 att->setValue(tempStr);
165 xercesc::XMLString::release(&tempStr);
166
167 return att;
168}

References G4GDMLWrite::doc, and G4InuclParticleNames::name().

Referenced by G4GDMLWrite::AddAuxInfo(), G4GDMLWriteStructure::AssemblyWrite(), G4GDMLWriteMaterials::AtomWrite(), BooleanWrite(), G4GDMLWriteStructure::BorderSurfaceCache(), G4GDMLWriteParamvol::Box_dimensionsWrite(), BoxWrite(), G4GDMLWriteParamvol::Cone_dimensionsWrite(), ConeWrite(), CutTubeWrite(), G4GDMLWriteStructure::DivisionvolWrite(), G4GDMLWriteMaterials::DWrite(), ElconeWrite(), G4GDMLWriteMaterials::ElementWrite(), G4GDMLWriteParamvol::Ellipsoid_dimensionsWrite(), EllipsoidWrite(), EltubeWrite(), GenericPolyconeWrite(), GenTrapWrite(), G4GDMLWriteParamvol::Hype_dimensionsWrite(), HypeWrite(), G4GDMLWriteMaterials::IsotopeWrite(), G4GDMLWriteMaterials::MaterialWrite(), G4GDMLWriteMaterials::MEEWrite(), MultiUnionWrite(), OpticalSurfaceWrite(), G4GDMLWriteParamvol::Orb_dimensionsWrite(), OrbWrite(), G4GDMLWriteParamvol::Para_dimensionsWrite(), ParaboloidWrite(), G4GDMLWriteParamvol::ParametersWrite(), G4GDMLWriteParamvol::ParamvolWrite(), ParaWrite(), G4GDMLWriteStructure::PhysvolWrite(), G4GDMLWriteParamvol::Polycone_dimensionsWrite(), PolyconeWrite(), G4GDMLWriteParamvol::Polyhedra_dimensionsWrite(), PolyhedraWrite(), G4GDMLWriteDefine::Position_vectorWrite(), G4GDMLWriteMaterials::PropertyConstWrite(), G4GDMLWriteMaterials::PropertyVectorWrite(), G4GDMLWriteMaterials::PropertyWrite(), PropertyWrite(), G4GDMLWriteMaterials::PWrite(), G4GDMLWriteStructure::ReplicavolWrite(), G4GDMLWriteDefine::Rotation_vectorWrite(), RZPointWrite(), G4GDMLWriteDefine::Scale_vectorWrite(), ScaledWrite(), G4GDMLWriteSetup::SetupWrite(), G4GDMLWriteStructure::SkinSurfaceCache(), G4GDMLWriteParamvol::Sphere_dimensionsWrite(), SphereWrite(), TessellatedWrite(), TetWrite(), G4GDMLWriteParamvol::Torus_dimensionsWrite(), TorusWrite(), G4GDMLWriteParamvol::Trap_dimensionsWrite(), TrapWrite(), G4GDMLWriteStructure::TraverseVolumeTree(), G4GDMLWriteParamvol::Trd_dimensionsWrite(), TrdWrite(), G4GDMLWriteParamvol::Tube_dimensionsWrite(), TubeWrite(), TwistedboxWrite(), TwistedtrapWrite(), TwistedtrdWrite(), TwistedtubsWrite(), G4GDMLWriteMaterials::TWrite(), G4GDMLWrite::Write(), XtruWrite(), and ZplaneWrite().

◆ NewElement()

xercesc::DOMElement * G4GDMLWrite::NewElement ( const G4String name)
protectedinherited

Definition at line 192 of file G4GDMLWrite.cc.

193{
194 XMLCh* tempStr = NULL;
195 tempStr = xercesc::XMLString::transcode(name);
196 xercesc::DOMElement* elem = doc->createElement(tempStr);
197 xercesc::XMLString::release(&tempStr);
198
199 return elem;
200}
#define elem(i, j)

References G4GDMLWrite::doc, elem, and G4InuclParticleNames::name().

Referenced by G4GDMLWrite::AddAuxInfo(), G4GDMLWriteStructure::AssemblyWrite(), G4GDMLWriteMaterials::AtomWrite(), BooleanWrite(), G4GDMLWriteStructure::BorderSurfaceCache(), G4GDMLWriteParamvol::Box_dimensionsWrite(), BoxWrite(), G4GDMLWriteParamvol::Cone_dimensionsWrite(), ConeWrite(), CutTubeWrite(), G4GDMLWriteDefine::DefineWrite(), G4GDMLWriteStructure::DivisionvolWrite(), G4GDMLWriteMaterials::DWrite(), ElconeWrite(), G4GDMLWriteMaterials::ElementWrite(), G4GDMLWriteParamvol::Ellipsoid_dimensionsWrite(), EllipsoidWrite(), EltubeWrite(), GenericPolyconeWrite(), GenTrapWrite(), G4GDMLWriteParamvol::Hype_dimensionsWrite(), HypeWrite(), G4GDMLWriteMaterials::IsotopeWrite(), G4GDMLWriteMaterials::MaterialsWrite(), G4GDMLWriteMaterials::MaterialWrite(), G4GDMLWriteMaterials::MEEWrite(), MultiUnionWrite(), OpticalSurfaceWrite(), G4GDMLWriteParamvol::Orb_dimensionsWrite(), OrbWrite(), G4GDMLWriteParamvol::Para_dimensionsWrite(), ParaboloidWrite(), G4GDMLWriteParamvol::ParametersWrite(), G4GDMLWriteParamvol::ParamvolWrite(), ParaWrite(), G4GDMLWriteStructure::PhysvolWrite(), G4GDMLWriteParamvol::Polycone_dimensionsWrite(), PolyconeWrite(), G4GDMLWriteParamvol::Polyhedra_dimensionsWrite(), PolyhedraWrite(), G4GDMLWriteDefine::Position_vectorWrite(), G4GDMLWriteMaterials::PropertyConstWrite(), G4GDMLWriteMaterials::PropertyVectorWrite(), G4GDMLWriteMaterials::PropertyWrite(), PropertyWrite(), G4GDMLWriteMaterials::PWrite(), G4GDMLWriteStructure::ReplicavolWrite(), G4GDMLWriteDefine::Rotation_vectorWrite(), RZPointWrite(), G4GDMLWriteDefine::Scale_vectorWrite(), ScaledWrite(), G4GDMLWriteSetup::SetupWrite(), G4GDMLWriteStructure::SkinSurfaceCache(), SolidsWrite(), G4GDMLWriteParamvol::Sphere_dimensionsWrite(), SphereWrite(), G4GDMLWriteStructure::StructureWrite(), TessellatedWrite(), TetWrite(), G4GDMLWriteParamvol::Torus_dimensionsWrite(), TorusWrite(), G4GDMLWriteParamvol::Trap_dimensionsWrite(), TrapWrite(), G4GDMLWriteStructure::TraverseVolumeTree(), G4GDMLWriteParamvol::Trd_dimensionsWrite(), TrdWrite(), G4GDMLWriteParamvol::Tube_dimensionsWrite(), TubeWrite(), TwistedboxWrite(), TwistedtrapWrite(), TwistedtrdWrite(), TwistedtubsWrite(), G4GDMLWriteMaterials::TWrite(), G4GDMLWrite::UserinfoWrite(), XtruWrite(), and ZplaneWrite().

◆ OpticalSurfaceWrite()

void G4GDMLWriteSolids::OpticalSurfaceWrite ( xercesc::DOMElement *  solElement,
const G4OpticalSurface * const  surf 
)
protected

Definition at line 1058 of file G4GDMLWriteSolids.cc.

1060{
1061 xercesc::DOMElement* optElement = NewElement("opticalsurface");
1062 G4OpticalSurfaceModel smodel = surf->GetModel();
1063 G4double sval =
1064 (smodel == glisur) ? surf->GetPolish() : surf->GetSigmaAlpha();
1065 const G4String& name = GenerateName(surf->GetName(), surf);
1066
1067 optElement->setAttributeNode(NewAttribute("name", name));
1068 optElement->setAttributeNode(NewAttribute("model", smodel));
1069 optElement->setAttributeNode(NewAttribute("finish", surf->GetFinish()));
1070 optElement->setAttributeNode(NewAttribute("type", surf->GetType()));
1071 optElement->setAttributeNode(NewAttribute("value", sval));
1072
1073 // Write any property attached to the optical surface...
1074 //
1075 if(surf->GetMaterialPropertiesTable())
1076 {
1077 PropertyWrite(optElement, surf);
1078 }
1079
1080 solElement->appendChild(optElement);
1081}
G4OpticalSurfaceModel
@ glisur
void PropertyWrite(xercesc::DOMElement *, const G4OpticalSurface *const)
G4OpticalSurfaceModel GetModel() const
G4double GetSigmaAlpha() const
G4OpticalSurfaceFinish GetFinish() const
G4double GetPolish() const
G4MaterialPropertiesTable * GetMaterialPropertiesTable() const
const G4String & GetName() const
const G4SurfaceType & GetType() const

References G4GDMLWrite::GenerateName(), G4OpticalSurface::GetFinish(), G4OpticalSurface::GetMaterialPropertiesTable(), G4OpticalSurface::GetModel(), G4SurfaceProperty::GetName(), G4OpticalSurface::GetPolish(), G4OpticalSurface::GetSigmaAlpha(), G4SurfaceProperty::GetType(), glisur, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), and PropertyWrite().

Referenced by G4GDMLWriteStructure::BorderSurfaceCache(), and G4GDMLWriteStructure::SkinSurfaceCache().

◆ OrbWrite()

void G4GDMLWriteSolids::OrbWrite ( xercesc::DOMElement *  solElement,
const G4Orb * const  orb 
)
protected

Definition at line 461 of file G4GDMLWriteSolids.cc.

463{
464 const G4String& name = GenerateName(orb->GetName(), orb);
465
466 xercesc::DOMElement* orbElement = NewElement("orb");
467 orbElement->setAttributeNode(NewAttribute("name", name));
468 orbElement->setAttributeNode(NewAttribute("r", orb->GetRadius() / mm));
469 orbElement->setAttributeNode(NewAttribute("lunit", "mm"));
470 solElement->appendChild(orbElement);
471}
G4double GetRadius() const

References G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4Orb::GetRadius(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ ParaboloidWrite()

void G4GDMLWriteSolids::ParaboloidWrite ( xercesc::DOMElement *  solElement,
const G4Paraboloid * const  paraboloid 
)
protected

Definition at line 501 of file G4GDMLWriteSolids.cc.

503{
504 const G4String& name = GenerateName(paraboloid->GetName(), paraboloid);
505
506 xercesc::DOMElement* paraboloidElement = NewElement("paraboloid");
507 paraboloidElement->setAttributeNode(NewAttribute("name", name));
508 paraboloidElement->setAttributeNode(
509 NewAttribute("rlo", paraboloid->GetRadiusMinusZ() / mm));
510 paraboloidElement->setAttributeNode(
511 NewAttribute("rhi", paraboloid->GetRadiusPlusZ() / mm));
512 paraboloidElement->setAttributeNode(
513 NewAttribute("dz", paraboloid->GetZHalfLength() / mm));
514 paraboloidElement->setAttributeNode(NewAttribute("lunit", "mm"));
515 solElement->appendChild(paraboloidElement);
516}
G4double GetRadiusPlusZ() const
G4double GetRadiusMinusZ() const
G4double GetZHalfLength() const

References G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4Paraboloid::GetRadiusMinusZ(), G4Paraboloid::GetRadiusPlusZ(), G4Paraboloid::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ ParaWrite()

void G4GDMLWriteSolids::ParaWrite ( xercesc::DOMElement *  solElement,
const G4Para * const  para 
)
protected

Definition at line 474 of file G4GDMLWriteSolids.cc.

476{
477 const G4String& name = GenerateName(para->GetName(), para);
478
479 const G4ThreeVector simaxis = para->GetSymAxis();
480 const G4double alpha = std::atan(para->GetTanAlpha());
481 const G4double phi = simaxis.phi();
482 const G4double theta = simaxis.theta();
483
484 xercesc::DOMElement* paraElement = NewElement("para");
485 paraElement->setAttributeNode(NewAttribute("name", name));
486 paraElement->setAttributeNode(
487 NewAttribute("x", 2.0 * para->GetXHalfLength() / mm));
488 paraElement->setAttributeNode(
489 NewAttribute("y", 2.0 * para->GetYHalfLength() / mm));
490 paraElement->setAttributeNode(
491 NewAttribute("z", 2.0 * para->GetZHalfLength() / mm));
492 paraElement->setAttributeNode(NewAttribute("alpha", alpha / degree));
493 paraElement->setAttributeNode(NewAttribute("theta", theta / degree));
494 paraElement->setAttributeNode(NewAttribute("phi", phi / degree));
495 paraElement->setAttributeNode(NewAttribute("aunit", "deg"));
496 paraElement->setAttributeNode(NewAttribute("lunit", "mm"));
497 solElement->appendChild(paraElement);
498}
static const G4double alpha
double phi() const
double theta() const
G4double GetTanAlpha() const
G4ThreeVector GetSymAxis() const
G4double GetYHalfLength() const
G4double GetZHalfLength() const
G4double GetXHalfLength() const

References alpha, degree, G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4Para::GetSymAxis(), G4Para::GetTanAlpha(), G4Para::GetXHalfLength(), G4Para::GetYHalfLength(), G4Para::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), CLHEP::Hep3Vector::phi(), and CLHEP::Hep3Vector::theta().

Referenced by AddSolid().

◆ PolyconeWrite()

void G4GDMLWriteSolids::PolyconeWrite ( xercesc::DOMElement *  solElement,
const G4Polycone * const  polycone 
)
protected

Definition at line 519 of file G4GDMLWriteSolids.cc.

521{
522 const G4String& name = GenerateName(polycone->GetName(), polycone);
523
524 xercesc::DOMElement* polyconeElement = NewElement("polycone");
525 polyconeElement->setAttributeNode(NewAttribute("name", name));
526 polyconeElement->setAttributeNode(NewAttribute(
527 "startphi", polycone->GetOriginalParameters()->Start_angle / degree));
528 polyconeElement->setAttributeNode(NewAttribute(
529 "deltaphi", polycone->GetOriginalParameters()->Opening_angle / degree));
530 polyconeElement->setAttributeNode(NewAttribute("aunit", "deg"));
531 polyconeElement->setAttributeNode(NewAttribute("lunit", "mm"));
532 solElement->appendChild(polyconeElement);
533
534 const std::size_t num_zplanes
535 = polycone->GetOriginalParameters()->Num_z_planes;
536 const G4double* z_array = polycone->GetOriginalParameters()->Z_values;
537 const G4double* rmin_array = polycone->GetOriginalParameters()->Rmin;
538 const G4double* rmax_array = polycone->GetOriginalParameters()->Rmax;
539
540 for(std::size_t i = 0; i < num_zplanes; ++i)
541 {
542 ZplaneWrite(polyconeElement, z_array[i], rmin_array[i], rmax_array[i]);
543 }
544}
void ZplaneWrite(xercesc::DOMElement *, const G4double &, const G4double &, const G4double &)
G4PolyconeHistorical * GetOriginalParameters() const

References degree, G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4Polycone::GetOriginalParameters(), G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), G4PolyconeHistorical::Num_z_planes, G4PolyconeHistorical::Opening_angle, G4PolyconeHistorical::Rmax, G4PolyconeHistorical::Rmin, G4PolyconeHistorical::Start_angle, G4PolyconeHistorical::Z_values, and ZplaneWrite().

Referenced by AddSolid().

◆ PolyhedraWrite()

void G4GDMLWriteSolids::PolyhedraWrite ( xercesc::DOMElement *  solElement,
const G4Polyhedra * const  polyhedra 
)
protected

Definition at line 572 of file G4GDMLWriteSolids.cc.

574{
575 const G4String& name = GenerateName(polyhedra->GetName(), polyhedra);
576 if(polyhedra->IsGeneric() == false)
577 {
578 xercesc::DOMElement* polyhedraElement = NewElement("polyhedra");
579 polyhedraElement->setAttributeNode(NewAttribute("name", name));
580 polyhedraElement->setAttributeNode(NewAttribute(
581 "startphi", polyhedra->GetOriginalParameters()->Start_angle / degree));
582 polyhedraElement->setAttributeNode(NewAttribute(
583 "deltaphi", polyhedra->GetOriginalParameters()->Opening_angle / degree));
584 polyhedraElement->setAttributeNode(
585 NewAttribute("numsides", polyhedra->GetOriginalParameters()->numSide));
586 polyhedraElement->setAttributeNode(NewAttribute("aunit", "deg"));
587 polyhedraElement->setAttributeNode(NewAttribute("lunit", "mm"));
588 solElement->appendChild(polyhedraElement);
589
590 const std::size_t num_zplanes
591 = polyhedra->GetOriginalParameters()->Num_z_planes;
592 const G4double* z_array = polyhedra->GetOriginalParameters()->Z_values;
593 const G4double* rmin_array = polyhedra->GetOriginalParameters()->Rmin;
594 const G4double* rmax_array = polyhedra->GetOriginalParameters()->Rmax;
595
596 const G4double convertRad =
597 std::cos(0.5 * polyhedra->GetOriginalParameters()->Opening_angle /
598 polyhedra->GetOriginalParameters()->numSide);
599
600 for(std::size_t i = 0; i < num_zplanes; ++i)
601 {
602 ZplaneWrite(polyhedraElement, z_array[i], rmin_array[i] * convertRad,
603 rmax_array[i] * convertRad);
604 }
605 }
606 else // generic polyhedra
607 {
608 xercesc::DOMElement* polyhedraElement = NewElement("genericPolyhedra");
609 polyhedraElement->setAttributeNode(NewAttribute("name", name));
610 polyhedraElement->setAttributeNode(NewAttribute(
611 "startphi", polyhedra->GetOriginalParameters()->Start_angle / degree));
612 polyhedraElement->setAttributeNode(NewAttribute(
613 "deltaphi", polyhedra->GetOriginalParameters()->Opening_angle / degree));
614 polyhedraElement->setAttributeNode(
615 NewAttribute("numsides", polyhedra->GetOriginalParameters()->numSide));
616 polyhedraElement->setAttributeNode(NewAttribute("aunit", "deg"));
617 polyhedraElement->setAttributeNode(NewAttribute("lunit", "mm"));
618 solElement->appendChild(polyhedraElement);
619
620 const std::size_t num_rzpoints = polyhedra->GetNumRZCorner();
621
622 for(std::size_t i = 0; i < num_rzpoints; ++i)
623 {
624 const G4double r_point = polyhedra->GetCorner(i).r;
625 const G4double z_point = polyhedra->GetCorner(i).z;
626 RZPointWrite(polyhedraElement, r_point, z_point);
627 }
628 }
629}
G4int GetNumRZCorner() const
G4PolyhedraHistorical * GetOriginalParameters() const
G4PolyhedraSideRZ GetCorner(const G4int index) const
G4bool IsGeneric() const

References degree, G4GDMLWrite::GenerateName(), G4Polyhedra::GetCorner(), G4VSolid::GetName(), G4Polyhedra::GetNumRZCorner(), G4Polyhedra::GetOriginalParameters(), G4Polyhedra::IsGeneric(), G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), G4PolyhedraHistorical::Num_z_planes, G4PolyhedraHistorical::numSide, G4PolyhedraHistorical::Opening_angle, G4PolyhedraSideRZ::r, G4PolyhedraHistorical::Rmax, G4PolyhedraHistorical::Rmin, RZPointWrite(), G4PolyhedraHistorical::Start_angle, G4PolyhedraSideRZ::z, G4PolyhedraHistorical::Z_values, and ZplaneWrite().

Referenced by AddSolid().

◆ Position_vectorWrite()

void G4GDMLWriteDefine::Position_vectorWrite ( xercesc::DOMElement *  element,
const G4String tag,
const G4String name,
const G4ThreeVector pos 
)
protectedinherited

Definition at line 118 of file G4GDMLWriteDefine.cc.

122{
123 const G4double x = (std::fabs(pos.x()) < kLinearPrecision) ? 0.0 : pos.x();
124 const G4double y = (std::fabs(pos.y()) < kLinearPrecision) ? 0.0 : pos.y();
125 const G4double z = (std::fabs(pos.z()) < kLinearPrecision) ? 0.0 : pos.z();
126
127 xercesc::DOMElement* positionElement = NewElement(tag);
128 positionElement->setAttributeNode(NewAttribute("name", name));
129 positionElement->setAttributeNode(NewAttribute("x", x / mm));
130 positionElement->setAttributeNode(NewAttribute("y", y / mm));
131 positionElement->setAttributeNode(NewAttribute("z", z / mm));
132 positionElement->setAttributeNode(NewAttribute("unit", "mm"));
133 element->appendChild(positionElement);
134}

References G4GDMLWriteDefine::kLinearPrecision, mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), and pos.

Referenced by G4GDMLWriteDefine::AddPosition(), G4GDMLWriteDefine::FirstpositionWrite(), and G4GDMLWriteDefine::PositionWrite().

◆ PositionWrite()

void G4GDMLWriteDefine::PositionWrite ( xercesc::DOMElement *  element,
const G4String name,
const G4ThreeVector pos 
)
inlineinherited

◆ PropertyConstWrite()

void G4GDMLWriteMaterials::PropertyConstWrite ( const G4String key,
const G4double  pval,
const G4MaterialPropertiesTable ptable 
)
protectedinherited

Definition at line 261 of file G4GDMLWriteMaterials.cc.

264{
265 const G4String matrixref = GenerateName(key, ptable);
266 xercesc::DOMElement* matrixElement = NewElement("matrix");
267 matrixElement->setAttributeNode(NewAttribute("name", matrixref));
268 matrixElement->setAttributeNode(NewAttribute("coldim", "1"));
269 std::ostringstream pvalues;
270
271 pvalues << pval;
272 matrixElement->setAttributeNode(NewAttribute("values", pvalues.str()));
273
274 defineElement->appendChild(matrixElement);
275}

References G4GDMLWriteDefine::defineElement, G4GDMLWrite::GenerateName(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by G4GDMLWriteMaterials::PropertyWrite().

◆ PropertyVectorWrite()

void G4GDMLWriteMaterials::PropertyVectorWrite ( const G4String key,
const G4PhysicsFreeVector * const  pvec 
)
protectedinherited

Definition at line 230 of file G4GDMLWriteMaterials.cc.

232{
233 for(std::size_t i = 0; i < propertyList.size(); ++i) // Check if property is
234 { // already in the list!
235 if(propertyList[i] == pvec)
236 {
237 return;
238 }
239 }
240 propertyList.push_back(pvec);
241
242 const G4String matrixref = GenerateName(key, pvec);
243 xercesc::DOMElement* matrixElement = NewElement("matrix");
244 matrixElement->setAttributeNode(NewAttribute("name", matrixref));
245 matrixElement->setAttributeNode(NewAttribute("coldim", "2"));
246 std::ostringstream pvalues;
247 for(std::size_t i = 0; i < pvec->GetVectorLength(); ++i)
248 {
249 if(i != 0)
250 {
251 pvalues << " ";
252 }
253 pvalues << pvec->Energy(i) << " " << (*pvec)[i];
254 }
255 matrixElement->setAttributeNode(NewAttribute("values", pvalues.str()));
256
257 defineElement->appendChild(matrixElement);
258}
G4double Energy(const std::size_t index) const
std::size_t GetVectorLength() const

References G4GDMLWriteDefine::defineElement, G4PhysicsVector::Energy(), G4GDMLWrite::GenerateName(), G4PhysicsVector::GetVectorLength(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), and G4GDMLWriteMaterials::propertyList.

Referenced by G4GDMLWriteMaterials::PropertyWrite(), and PropertyWrite().

◆ PropertyWrite() [1/2]

void G4GDMLWriteMaterials::PropertyWrite ( xercesc::DOMElement *  matElement,
const G4Material * const  mat 
)
protectedinherited

Definition at line 278 of file G4GDMLWriteMaterials.cc.

280{
281 xercesc::DOMElement* propElement;
283
284 auto pvec = ptable->GetProperties();
285 auto cvec = ptable->GetConstProperties();
286
287 for(size_t i = 0; i < pvec.size(); ++i)
288 {
289 if (pvec[i] != nullptr)
290 {
291 propElement = NewElement("property");
292 propElement->setAttributeNode(
293 NewAttribute("name", ptable->GetMaterialPropertyNames()[i]));
294 propElement->setAttributeNode(NewAttribute(
295 "ref", GenerateName(ptable->GetMaterialPropertyNames()[i],
296 pvec[i])));
298 pvec[i]);
299 matElement->appendChild(propElement);
300 }
301 }
302
303 for(size_t i = 0; i < cvec.size(); ++i)
304 {
305 if (cvec[i].second == true)
306 {
307 propElement = NewElement("property");
308 propElement->setAttributeNode(NewAttribute(
309 "name", ptable->GetMaterialConstPropertyNames()[i]));
310 propElement->setAttributeNode(NewAttribute(
312 ptable)));
314 cvec[i].first, ptable);
315 matElement->appendChild(propElement);
316 }
317 }
318}
static constexpr double second
Definition: G4SIunits.hh:137
void PropertyVectorWrite(const G4String &, const G4PhysicsFreeVector *const)
void PropertyConstWrite(const G4String &, const G4double, const G4MaterialPropertiesTable *)
const std::vector< G4String > & GetMaterialConstPropertyNames() const
const std::vector< std::pair< G4double, G4bool > > & GetConstProperties() const
const std::vector< G4MaterialPropertyVector * > & GetProperties() const
const std::vector< G4String > & GetMaterialPropertyNames() const

References G4GDMLWrite::GenerateName(), G4MaterialPropertiesTable::GetConstProperties(), G4MaterialPropertiesTable::GetMaterialConstPropertyNames(), G4Material::GetMaterialPropertiesTable(), G4MaterialPropertiesTable::GetMaterialPropertyNames(), G4MaterialPropertiesTable::GetProperties(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), G4GDMLWriteMaterials::PropertyConstWrite(), G4GDMLWriteMaterials::PropertyVectorWrite(), and second.

Referenced by G4GDMLWriteMaterials::MaterialWrite().

◆ PropertyWrite() [2/2]

void G4GDMLWriteSolids::PropertyWrite ( xercesc::DOMElement *  optElement,
const G4OpticalSurface * const  surf 
)
protected

Definition at line 1084 of file G4GDMLWriteSolids.cc.

1086{
1087 xercesc::DOMElement* propElement;
1089 auto pvec = ptable->GetProperties();
1090 auto cvec = ptable->GetConstProperties();
1091
1092 for(size_t i = 0; i < pvec.size(); ++i)
1093 {
1094 if(pvec[i] != nullptr) {
1095 propElement = NewElement("property");
1096 propElement->setAttributeNode(
1097 NewAttribute("name", ptable->GetMaterialPropertyNames()[i]));
1098 propElement->setAttributeNode(NewAttribute(
1099 "ref", GenerateName(ptable->GetMaterialPropertyNames()[i],
1100 pvec[i])));
1102 pvec[i]);
1103 optElement->appendChild(propElement);
1104 }
1105 }
1106 for(size_t i = 0; i < cvec.size(); ++i)
1107 {
1108 if (cvec[i].second == true) {
1109 propElement = NewElement("property");
1110 propElement->setAttributeNode(NewAttribute(
1111 "name", ptable->GetMaterialConstPropertyNames()[i]));
1112 propElement->setAttributeNode(NewAttribute(
1113 "ref", ptable->GetMaterialConstPropertyNames()[i]));
1114 xercesc::DOMElement* constElement = NewElement("constant");
1115 constElement->setAttributeNode(NewAttribute(
1116 "name", ptable->GetMaterialConstPropertyNames()[i]));
1117 constElement->setAttributeNode(NewAttribute("value", cvec[i].first));
1118 defineElement->appendChild(constElement);
1119 optElement->appendChild(propElement);
1120 }
1121 }
1122}

References G4GDMLWriteDefine::defineElement, G4GDMLWrite::GenerateName(), G4MaterialPropertiesTable::GetConstProperties(), G4MaterialPropertiesTable::GetMaterialConstPropertyNames(), G4OpticalSurface::GetMaterialPropertiesTable(), G4MaterialPropertiesTable::GetMaterialPropertyNames(), G4MaterialPropertiesTable::GetProperties(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), G4GDMLWriteMaterials::PropertyVectorWrite(), and second.

Referenced by OpticalSurfaceWrite().

◆ PvolumeMap()

G4GDMLWrite::PhysVolumeMapType & G4GDMLWrite::PvolumeMap ( )
protectedinherited

Definition at line 66 of file G4GDMLWrite.cc.

67{
68 static PhysVolumeMapType instance;
69 return instance;
70}
std::map< const G4VPhysicalVolume *, G4String > PhysVolumeMapType
Definition: G4GDMLWrite.hh:54

Referenced by G4GDMLWrite::AddModule(), and G4GDMLWrite::Modularize().

◆ PWrite()

void G4GDMLWriteMaterials::PWrite ( xercesc::DOMElement *  element,
const G4double P 
)
protectedinherited

Definition at line 73 of file G4GDMLWriteMaterials.cc.

75{
76 xercesc::DOMElement* PElement = NewElement("P");
77 PElement->setAttributeNode(NewAttribute("unit", "pascal"));
78 PElement->setAttributeNode(NewAttribute("value", P / hep_pascal));
79 element->appendChild(PElement);
80}
static constexpr double hep_pascal
Definition: G4SIunits.hh:222
static double P[]

References hep_pascal, G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), and P.

Referenced by G4GDMLWriteMaterials::MaterialWrite().

◆ Rotation_vectorWrite()

void G4GDMLWriteDefine::Rotation_vectorWrite ( xercesc::DOMElement *  element,
const G4String tag,
const G4String name,
const G4ThreeVector rot 
)
protectedinherited

Definition at line 99 of file G4GDMLWriteDefine.cc.

103{
104 const G4double x = (std::fabs(rot.x()) < kAngularPrecision) ? 0.0 : rot.x();
105 const G4double y = (std::fabs(rot.y()) < kAngularPrecision) ? 0.0 : rot.y();
106 const G4double z = (std::fabs(rot.z()) < kAngularPrecision) ? 0.0 : rot.z();
107
108 xercesc::DOMElement* rotationElement = NewElement(tag);
109 rotationElement->setAttributeNode(NewAttribute("name", name));
110 rotationElement->setAttributeNode(NewAttribute("x", x / degree));
111 rotationElement->setAttributeNode(NewAttribute("y", y / degree));
112 rotationElement->setAttributeNode(NewAttribute("z", z / degree));
113 rotationElement->setAttributeNode(NewAttribute("unit", "deg"));
114 element->appendChild(rotationElement);
115}

References degree, G4GDMLWriteDefine::kAngularPrecision, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by G4GDMLWriteDefine::FirstrotationWrite(), and G4GDMLWriteDefine::RotationWrite().

◆ RotationWrite()

void G4GDMLWriteDefine::RotationWrite ( xercesc::DOMElement *  element,
const G4String name,
const G4ThreeVector rot 
)
inlineinherited

◆ RZPointWrite()

void G4GDMLWriteSolids::RZPointWrite ( xercesc::DOMElement *  element,
const G4double r,
const G4double z 
)
protected

Definition at line 1048 of file G4GDMLWriteSolids.cc.

1050{
1051 xercesc::DOMElement* rzpointElement = NewElement("rzpoint");
1052 rzpointElement->setAttributeNode(NewAttribute("r", r / mm));
1053 rzpointElement->setAttributeNode(NewAttribute("z", z / mm));
1054 element->appendChild(rzpointElement);
1055}

References mm, G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by GenericPolyconeWrite(), and PolyhedraWrite().

◆ Scale_vectorWrite()

void G4GDMLWriteDefine::Scale_vectorWrite ( xercesc::DOMElement *  element,
const G4String tag,
const G4String name,
const G4ThreeVector scl 
)
protectedinherited

Definition at line 78 of file G4GDMLWriteDefine.cc.

82{
83 const G4double x =
84 (std::fabs(scl.x() - 1.0) < kRelativePrecision) ? 1.0 : scl.x();
85 const G4double y =
86 (std::fabs(scl.y() - 1.0) < kRelativePrecision) ? 1.0 : scl.y();
87 const G4double z =
88 (std::fabs(scl.z() - 1.0) < kRelativePrecision) ? 1.0 : scl.z();
89
90 xercesc::DOMElement* scaleElement = NewElement(tag);
91 scaleElement->setAttributeNode(NewAttribute("name", name));
92 scaleElement->setAttributeNode(NewAttribute("x", x));
93 scaleElement->setAttributeNode(NewAttribute("y", y));
94 scaleElement->setAttributeNode(NewAttribute("z", z));
95 element->appendChild(scaleElement);
96}
static const G4double kRelativePrecision

References G4GDMLWriteDefine::kRelativePrecision, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by G4GDMLWriteDefine::ScaleWrite().

◆ ScaledWrite()

void G4GDMLWriteSolids::ScaledWrite ( xercesc::DOMElement *  solElement,
const G4ScaledSolid * const  scaled 
)
protected

Definition at line 259 of file G4GDMLWriteSolids.cc.

261{
262 G4String tag("scaledSolid");
263
264 G4VSolid* solid = const_cast<G4VSolid*>(scaled->GetUnscaledSolid());
265 G4Scale3D scale = scaled->GetScaleTransform();
266 G4ThreeVector sclVector = G4ThreeVector(scale.xx(), scale.yy(), scale.zz());
267
268 AddSolid(solid); // Add the constituent solid!
269
270 const G4String& name = GenerateName(scaled->GetName(), scaled);
271 const G4String& solidref = GenerateName(solid->GetName(), solid);
272
273 xercesc::DOMElement* scaledElement = NewElement(tag);
274 scaledElement->setAttributeNode(NewAttribute("name", name));
275
276 xercesc::DOMElement* solidElement = NewElement("solidref");
277 solidElement->setAttributeNode(NewAttribute("ref", solidref));
278 scaledElement->appendChild(solidElement);
279
280 if((std::fabs(scale.xx()) > kLinearPrecision) &&
281 (std::fabs(scale.yy()) > kLinearPrecision) &&
282 (std::fabs(scale.zz()) > kLinearPrecision))
283 {
284 ScaleWrite(scaledElement, name + "_scl", sclVector);
285 }
286
287 solElement->appendChild(scaledElement);
288 // Add the scaled solid AFTER its constituent solid!
289}
void ScaleWrite(xercesc::DOMElement *element, const G4String &name, const G4ThreeVector &scl)
G4VSolid * GetUnscaledSolid() const
G4Scale3D GetScaleTransform() const

References AddSolid(), G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4ScaledSolid::GetScaleTransform(), G4ScaledSolid::GetUnscaledSolid(), G4GDMLWriteDefine::kLinearPrecision, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), G4GDMLWriteDefine::ScaleWrite(), HepGeom::Transform3D::xx(), HepGeom::Transform3D::yy(), and HepGeom::Transform3D::zz().

Referenced by AddSolid().

◆ ScaleWrite()

void G4GDMLWriteDefine::ScaleWrite ( xercesc::DOMElement *  element,
const G4String name,
const G4ThreeVector scl 
)
inlineinherited

Definition at line 48 of file G4GDMLWriteDefine.hh.

50 {
51 Scale_vectorWrite(element, "scale", name, scl);
52 }
void Scale_vectorWrite(xercesc::DOMElement *, const G4String &, const G4String &, const G4ThreeVector &)

References G4InuclParticleNames::name(), and G4GDMLWriteDefine::Scale_vectorWrite().

Referenced by G4GDMLWriteStructure::PhysvolWrite(), ScaledWrite(), and G4GDMLWriteStructure::TraverseVolumeTree().

◆ SetAddPointerToName()

void G4GDMLWrite::SetAddPointerToName ( G4bool  set)
staticinherited

Definition at line 410 of file G4GDMLWrite.cc.

411{
412 addPointerToName = set;
413}

References G4GDMLWrite::addPointerToName.

◆ SetOutputFileOverwrite()

void G4GDMLWrite::SetOutputFileOverwrite ( G4bool  flag)
inherited

Definition at line 404 of file G4GDMLWrite.cc.

405{
406 overwriteOutputFile = flag;
407}
G4bool overwriteOutputFile
Definition: G4GDMLWrite.hh:137

References G4GDMLWrite::overwriteOutputFile.

◆ SetupWrite()

virtual void G4GDMLWrite::SetupWrite ( xercesc::DOMElement *  ,
const G4LogicalVolume * const   
)
pure virtualinherited

Implemented in G4GDMLWriteSetup.

Referenced by G4GDMLWrite::Write().

◆ SolidsWrite()

void G4GDMLWriteSolids::SolidsWrite ( xercesc::DOMElement *  gdmlElement)
virtual

Implements G4GDMLWrite.

Definition at line 1125 of file G4GDMLWriteSolids.cc.

1126{
1127#ifdef G4VERBOSE
1128 G4cout << "G4GDML: Writing solids..." << G4endl;
1129#endif
1130 solidsElement = NewElement("solids");
1131 gdmlElement->appendChild(solidsElement);
1132
1133 solidList.clear();
1134}

References G4cout, G4endl, G4GDMLWrite::NewElement(), solidList, and solidsElement.

◆ SphereWrite()

void G4GDMLWriteSolids::SphereWrite ( xercesc::DOMElement *  solElement,
const G4Sphere * const  sphere 
)
protected

Definition at line 632 of file G4GDMLWriteSolids.cc.

634{
635 const G4String& name = GenerateName(sphere->GetName(), sphere);
636
637 xercesc::DOMElement* sphereElement = NewElement("sphere");
638 sphereElement->setAttributeNode(NewAttribute("name", name));
639 sphereElement->setAttributeNode(
640 NewAttribute("rmin", sphere->GetInnerRadius() / mm));
641 sphereElement->setAttributeNode(
642 NewAttribute("rmax", sphere->GetOuterRadius() / mm));
643 sphereElement->setAttributeNode(
644 NewAttribute("startphi", sphere->GetStartPhiAngle() / degree));
645 sphereElement->setAttributeNode(
646 NewAttribute("deltaphi", sphere->GetDeltaPhiAngle() / degree));
647 sphereElement->setAttributeNode(
648 NewAttribute("starttheta", sphere->GetStartThetaAngle() / degree));
649 sphereElement->setAttributeNode(
650 NewAttribute("deltatheta", sphere->GetDeltaThetaAngle() / degree));
651 sphereElement->setAttributeNode(NewAttribute("aunit", "deg"));
652 sphereElement->setAttributeNode(NewAttribute("lunit", "mm"));
653 solElement->appendChild(sphereElement);
654}
G4double GetStartPhiAngle() const
G4double GetDeltaPhiAngle() const
G4double GetInnerRadius() const
G4double GetOuterRadius() const
G4double GetDeltaThetaAngle() const
G4double GetStartThetaAngle() const

References degree, G4GDMLWrite::GenerateName(), G4Sphere::GetDeltaPhiAngle(), G4Sphere::GetDeltaThetaAngle(), G4Sphere::GetInnerRadius(), G4VSolid::GetName(), G4Sphere::GetOuterRadius(), G4Sphere::GetStartPhiAngle(), G4Sphere::GetStartThetaAngle(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ StructureWrite()

virtual void G4GDMLWrite::StructureWrite ( xercesc::DOMElement *  )
pure virtualinherited

Implemented in G4GDMLWriteStructure.

Referenced by G4GDMLWrite::Write().

◆ SurfacesWrite()

virtual void G4GDMLWrite::SurfacesWrite ( )
pure virtualinherited

Implemented in G4GDMLWriteStructure.

Referenced by G4GDMLWrite::Write().

◆ TessellatedWrite()

void G4GDMLWriteSolids::TessellatedWrite ( xercesc::DOMElement *  solElement,
const G4TessellatedSolid * const  tessellated 
)
protected

Definition at line 657 of file G4GDMLWriteSolids.cc.

659{
660 const G4String& solid_name = tessellated->GetName();
661 const G4String& name = GenerateName(solid_name, tessellated);
662
663 xercesc::DOMElement* tessellatedElement = NewElement("tessellated");
664 tessellatedElement->setAttributeNode(NewAttribute("name", name));
665 tessellatedElement->setAttributeNode(NewAttribute("aunit", "deg"));
666 tessellatedElement->setAttributeNode(NewAttribute("lunit", "mm"));
667 solElement->appendChild(tessellatedElement);
668
669 std::map<G4ThreeVector, G4String, G4ThreeVectorCompare> vertexMap;
670
671 const std::size_t NumFacets = tessellated->GetNumberOfFacets();
672 std::size_t NumVertex = 0;
673
674 for(std::size_t i = 0; i < NumFacets; ++i)
675 {
676 const G4VFacet* facet = tessellated->GetFacet(i);
677 const size_t NumVertexPerFacet = facet->GetNumberOfVertices();
678
679 G4String FacetTag;
680
681 if(NumVertexPerFacet == 3)
682 {
683 FacetTag = "triangular";
684 }
685 else if(NumVertexPerFacet == 4)
686 {
687 FacetTag = "quadrangular";
688 }
689 else
690 {
691 G4Exception("G4GDMLWriteSolids::TessellatedWrite()", "InvalidSetup",
692 FatalException, "Facet should contain 3 or 4 vertices!");
693 }
694
695 xercesc::DOMElement* facetElement = NewElement(FacetTag);
696 tessellatedElement->appendChild(facetElement);
697
698 for(std::size_t j = 0; j < NumVertexPerFacet; ++j)
699 {
700 std::stringstream name_stream;
701 std::stringstream ref_stream;
702
703 name_stream << "vertex" << (j + 1);
704 ref_stream << solid_name << "_v" << NumVertex;
705
706 const G4String& fname = name_stream.str(); // facet's tag variable
707 G4String ref = ref_stream.str(); // vertex tag to be associated
708
709 // Now search for the existance of the current vertex in the
710 // map of cached vertices. If existing, do NOT store it as
711 // position in the GDML file, so avoiding duplication; otherwise
712 // cache it in the local map and add it as position in the
713 // "define" section of the GDML file.
714
715 const G4ThreeVector& vertex = facet->GetVertex(j);
716
717 if(vertexMap.find(vertex) != vertexMap.cend()) // Vertex is cached
718 {
719 ref = vertexMap[vertex]; // Set the proper tag for it
720 }
721 else // Vertex not found
722 {
723 vertexMap.insert(std::make_pair(vertex, ref)); // Cache vertex and ...
724 AddPosition(ref, vertex); // ... add it to define section!
725 ++NumVertex;
726 }
727
728 // Now create association of the vertex with its facet
729 //
730 facetElement->setAttributeNode(NewAttribute(fname, ref));
731 }
732 }
733}
void AddPosition(const G4String &name, const G4ThreeVector &pos)
G4int GetNumberOfFacets() const
G4VFacet * GetFacet(G4int i) const
virtual G4ThreeVector GetVertex(G4int i) const =0
virtual G4int GetNumberOfVertices() const =0

References G4GDMLWriteDefine::AddPosition(), FatalException, test::fname, G4Exception(), G4GDMLWrite::GenerateName(), G4TessellatedSolid::GetFacet(), G4VSolid::GetName(), G4TessellatedSolid::GetNumberOfFacets(), G4VFacet::GetNumberOfVertices(), G4VFacet::GetVertex(), G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ TetWrite()

void G4GDMLWriteSolids::TetWrite ( xercesc::DOMElement *  solElement,
const G4Tet * const  tet 
)
protected

Definition at line 736 of file G4GDMLWriteSolids.cc.

738{
739 const G4String& solid_name = tet->GetName();
740 const G4String& name = GenerateName(solid_name, tet);
741
742 std::vector<G4ThreeVector> vertexList = tet->GetVertices();
743
744 xercesc::DOMElement* tetElement = NewElement("tet");
745 tetElement->setAttributeNode(NewAttribute("name", name));
746 tetElement->setAttributeNode(NewAttribute("vertex1", solid_name + "_v1"));
747 tetElement->setAttributeNode(NewAttribute("vertex2", solid_name + "_v2"));
748 tetElement->setAttributeNode(NewAttribute("vertex3", solid_name + "_v3"));
749 tetElement->setAttributeNode(NewAttribute("vertex4", solid_name + "_v4"));
750 tetElement->setAttributeNode(NewAttribute("lunit", "mm"));
751 solElement->appendChild(tetElement);
752
753 AddPosition(solid_name + "_v1", vertexList[0]);
754 AddPosition(solid_name + "_v2", vertexList[1]);
755 AddPosition(solid_name + "_v3", vertexList[2]);
756 AddPosition(solid_name + "_v4", vertexList[3]);
757}
void GetVertices(G4ThreeVector &anchor, G4ThreeVector &p1, G4ThreeVector &p2, G4ThreeVector &p3) const
Definition: G4Tet.cc:286

References G4GDMLWriteDefine::AddPosition(), G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4Tet::GetVertices(), G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ TorusWrite()

void G4GDMLWriteSolids::TorusWrite ( xercesc::DOMElement *  solElement,
const G4Torus * const  torus 
)
protected

Definition at line 760 of file G4GDMLWriteSolids.cc.

762{
763 const G4String& name = GenerateName(torus->GetName(), torus);
764
765 xercesc::DOMElement* torusElement = NewElement("torus");
766 torusElement->setAttributeNode(NewAttribute("name", name));
767 torusElement->setAttributeNode(NewAttribute("rmin", torus->GetRmin() / mm));
768 torusElement->setAttributeNode(NewAttribute("rmax", torus->GetRmax() / mm));
769 torusElement->setAttributeNode(NewAttribute("rtor", torus->GetRtor() / mm));
770 torusElement->setAttributeNode(
771 NewAttribute("startphi", torus->GetSPhi() / degree));
772 torusElement->setAttributeNode(
773 NewAttribute("deltaphi", torus->GetDPhi() / degree));
774 torusElement->setAttributeNode(NewAttribute("aunit", "deg"));
775 torusElement->setAttributeNode(NewAttribute("lunit", "mm"));
776 solElement->appendChild(torusElement);
777}
G4double GetDPhi() const
G4double GetRmin() const
G4double GetRtor() const
G4double GetRmax() const
G4double GetSPhi() const

References degree, G4GDMLWrite::GenerateName(), G4Torus::GetDPhi(), G4VSolid::GetName(), G4Torus::GetRmax(), G4Torus::GetRmin(), G4Torus::GetRtor(), G4Torus::GetSPhi(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ TrapWrite()

void G4GDMLWriteSolids::TrapWrite ( xercesc::DOMElement *  solElement,
const G4Trap * const  trap 
)
protected

Definition at line 812 of file G4GDMLWriteSolids.cc.

814{
815 const G4String& name = GenerateName(trap->GetName(), trap);
816
817 const G4ThreeVector& simaxis = trap->GetSymAxis();
818 const G4double phi = simaxis.phi();
819 const G4double theta = simaxis.theta();
820 const G4double alpha1 = std::atan(trap->GetTanAlpha1());
821 const G4double alpha2 = std::atan(trap->GetTanAlpha2());
822
823 xercesc::DOMElement* trapElement = NewElement("trap");
824 trapElement->setAttributeNode(NewAttribute("name", name));
825 trapElement->setAttributeNode(
826 NewAttribute("z", 2.0 * trap->GetZHalfLength() / mm));
827 trapElement->setAttributeNode(NewAttribute("theta", theta / degree));
828 trapElement->setAttributeNode(NewAttribute("phi", phi / degree));
829 trapElement->setAttributeNode(
830 NewAttribute("y1", 2.0 * trap->GetYHalfLength1() / mm));
831 trapElement->setAttributeNode(
832 NewAttribute("x1", 2.0 * trap->GetXHalfLength1() / mm));
833 trapElement->setAttributeNode(
834 NewAttribute("x2", 2.0 * trap->GetXHalfLength2() / mm));
835 trapElement->setAttributeNode(NewAttribute("alpha1", alpha1 / degree));
836 trapElement->setAttributeNode(
837 NewAttribute("y2", 2.0 * trap->GetYHalfLength2() / mm));
838 trapElement->setAttributeNode(
839 NewAttribute("x3", 2.0 * trap->GetXHalfLength3() / mm));
840 trapElement->setAttributeNode(
841 NewAttribute("x4", 2.0 * trap->GetXHalfLength4() / mm));
842 trapElement->setAttributeNode(NewAttribute("alpha2", alpha2 / degree));
843 trapElement->setAttributeNode(NewAttribute("aunit", "deg"));
844 trapElement->setAttributeNode(NewAttribute("lunit", "mm"));
845 solElement->appendChild(trapElement);
846}
const G4double alpha2
G4double GetYHalfLength1() const
G4double GetTanAlpha2() const
G4double GetXHalfLength2() const
G4ThreeVector GetSymAxis() const
G4double GetXHalfLength4() const
G4double GetZHalfLength() const
G4double GetYHalfLength2() const
G4double GetTanAlpha1() const
G4double GetXHalfLength3() const
G4double GetXHalfLength1() const

References alpha2, degree, G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4Trap::GetSymAxis(), G4Trap::GetTanAlpha1(), G4Trap::GetTanAlpha2(), G4Trap::GetXHalfLength1(), G4Trap::GetXHalfLength2(), G4Trap::GetXHalfLength3(), G4Trap::GetXHalfLength4(), G4Trap::GetYHalfLength1(), G4Trap::GetYHalfLength2(), G4Trap::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), CLHEP::Hep3Vector::phi(), and CLHEP::Hep3Vector::theta().

Referenced by AddSolid().

◆ TraverseVolumeTree()

virtual G4Transform3D G4GDMLWrite::TraverseVolumeTree ( const G4LogicalVolume * const  ,
const  G4int 
)
pure virtualinherited

Implemented in G4GDMLWriteStructure.

Referenced by G4GDMLWrite::Write().

◆ TrdWrite()

void G4GDMLWriteSolids::TrdWrite ( xercesc::DOMElement *  solElement,
const G4Trd * const  trd 
)
protected

Definition at line 849 of file G4GDMLWriteSolids.cc.

851{
852 const G4String& name = GenerateName(trd->GetName(), trd);
853
854 xercesc::DOMElement* trdElement = NewElement("trd");
855 trdElement->setAttributeNode(NewAttribute("name", name));
856 trdElement->setAttributeNode(
857 NewAttribute("x1", 2.0 * trd->GetXHalfLength1() / mm));
858 trdElement->setAttributeNode(
859 NewAttribute("x2", 2.0 * trd->GetXHalfLength2() / mm));
860 trdElement->setAttributeNode(
861 NewAttribute("y1", 2.0 * trd->GetYHalfLength1() / mm));
862 trdElement->setAttributeNode(
863 NewAttribute("y2", 2.0 * trd->GetYHalfLength2() / mm));
864 trdElement->setAttributeNode(
865 NewAttribute("z", 2.0 * trd->GetZHalfLength() / mm));
866 trdElement->setAttributeNode(NewAttribute("lunit", "mm"));
867 solElement->appendChild(trdElement);
868}
G4double GetXHalfLength2() const
G4double GetYHalfLength2() const
G4double GetXHalfLength1() const
G4double GetYHalfLength1() const
G4double GetZHalfLength() const

References G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4Trd::GetXHalfLength1(), G4Trd::GetXHalfLength2(), G4Trd::GetYHalfLength1(), G4Trd::GetYHalfLength2(), G4Trd::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ TubeWrite()

void G4GDMLWriteSolids::TubeWrite ( xercesc::DOMElement *  solElement,
const G4Tubs * const  tube 
)
protected

Definition at line 871 of file G4GDMLWriteSolids.cc.

873{
874 const G4String& name = GenerateName(tube->GetName(), tube);
875
876 xercesc::DOMElement* tubeElement = NewElement("tube");
877 tubeElement->setAttributeNode(NewAttribute("name", name));
878 tubeElement->setAttributeNode(
879 NewAttribute("rmin", tube->GetInnerRadius() / mm));
880 tubeElement->setAttributeNode(
881 NewAttribute("rmax", tube->GetOuterRadius() / mm));
882 tubeElement->setAttributeNode(
883 NewAttribute("z", 2.0 * tube->GetZHalfLength() / mm));
884 tubeElement->setAttributeNode(
885 NewAttribute("startphi", tube->GetStartPhiAngle() / degree));
886 tubeElement->setAttributeNode(
887 NewAttribute("deltaphi", tube->GetDeltaPhiAngle() / degree));
888 tubeElement->setAttributeNode(NewAttribute("aunit", "deg"));
889 tubeElement->setAttributeNode(NewAttribute("lunit", "mm"));
890 solElement->appendChild(tubeElement);
891}
G4double GetZHalfLength() const
G4double GetInnerRadius() const
G4double GetOuterRadius() const
G4double GetStartPhiAngle() const
G4double GetDeltaPhiAngle() const

References degree, G4GDMLWrite::GenerateName(), G4Tubs::GetDeltaPhiAngle(), G4Tubs::GetInnerRadius(), G4VSolid::GetName(), G4Tubs::GetOuterRadius(), G4Tubs::GetStartPhiAngle(), G4Tubs::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ TwistedboxWrite()

void G4GDMLWriteSolids::TwistedboxWrite ( xercesc::DOMElement *  solElement,
const G4TwistedBox * const  twistedbox 
)
protected

Definition at line 929 of file G4GDMLWriteSolids.cc.

931{
932 const G4String& name = GenerateName(twistedbox->GetName(), twistedbox);
933
934 xercesc::DOMElement* twistedboxElement = NewElement("twistedbox");
935 twistedboxElement->setAttributeNode(NewAttribute("name", name));
936 twistedboxElement->setAttributeNode(
937 NewAttribute("x", 2.0 * twistedbox->GetXHalfLength() / mm));
938 twistedboxElement->setAttributeNode(
939 NewAttribute("y", 2.0 * twistedbox->GetYHalfLength() / mm));
940 twistedboxElement->setAttributeNode(
941 NewAttribute("z", 2.0 * twistedbox->GetZHalfLength() / mm));
942 twistedboxElement->setAttributeNode(
943 NewAttribute("PhiTwist", twistedbox->GetPhiTwist() / degree));
944 twistedboxElement->setAttributeNode(NewAttribute("aunit", "deg"));
945 twistedboxElement->setAttributeNode(NewAttribute("lunit", "mm"));
946 solElement->appendChild(twistedboxElement);
947}
G4double GetPhiTwist() const
Definition: G4TwistedBox.hh:65
G4double GetXHalfLength() const
Definition: G4TwistedBox.hh:62
G4double GetZHalfLength() const
Definition: G4TwistedBox.hh:64
G4double GetYHalfLength() const
Definition: G4TwistedBox.hh:63

References degree, G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4TwistedBox::GetPhiTwist(), G4TwistedBox::GetXHalfLength(), G4TwistedBox::GetYHalfLength(), G4TwistedBox::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ TwistedtrapWrite()

void G4GDMLWriteSolids::TwistedtrapWrite ( xercesc::DOMElement *  solElement,
const G4TwistedTrap * const  twistedtrap 
)
protected

Definition at line 950 of file G4GDMLWriteSolids.cc.

952{
953 const G4String& name = GenerateName(twistedtrap->GetName(), twistedtrap);
954
955 xercesc::DOMElement* twistedtrapElement = NewElement("twistedtrap");
956 twistedtrapElement->setAttributeNode(NewAttribute("name", name));
957 twistedtrapElement->setAttributeNode(
958 NewAttribute("y1", 2.0 * twistedtrap->GetY1HalfLength() / mm));
959 twistedtrapElement->setAttributeNode(
960 NewAttribute("x1", 2.0 * twistedtrap->GetX1HalfLength() / mm));
961 twistedtrapElement->setAttributeNode(
962 NewAttribute("x2", 2.0 * twistedtrap->GetX2HalfLength() / mm));
963 twistedtrapElement->setAttributeNode(
964 NewAttribute("y2", 2.0 * twistedtrap->GetY2HalfLength() / mm));
965 twistedtrapElement->setAttributeNode(
966 NewAttribute("x3", 2.0 * twistedtrap->GetX3HalfLength() / mm));
967 twistedtrapElement->setAttributeNode(
968 NewAttribute("x4", 2.0 * twistedtrap->GetX4HalfLength() / mm));
969 twistedtrapElement->setAttributeNode(
970 NewAttribute("z", 2.0 * twistedtrap->GetZHalfLength() / mm));
971 twistedtrapElement->setAttributeNode(
972 NewAttribute("Alph", twistedtrap->GetTiltAngleAlpha() / degree));
973 twistedtrapElement->setAttributeNode(
974 NewAttribute("Theta", twistedtrap->GetPolarAngleTheta() / degree));
975 twistedtrapElement->setAttributeNode(
976 NewAttribute("Phi", twistedtrap->GetAzimuthalAnglePhi() / degree));
977 twistedtrapElement->setAttributeNode(
978 NewAttribute("PhiTwist", twistedtrap->GetPhiTwist() / degree));
979 twistedtrapElement->setAttributeNode(NewAttribute("aunit", "deg"));
980 twistedtrapElement->setAttributeNode(NewAttribute("lunit", "mm"));
981
982 solElement->appendChild(twistedtrapElement);
983}
G4double GetPolarAngleTheta() const
G4double GetAzimuthalAnglePhi() const
G4double GetTiltAngleAlpha() const
G4double GetZHalfLength() const
G4double GetX1HalfLength() const
G4double GetX2HalfLength() const
G4double GetX3HalfLength() const
G4double GetX4HalfLength() const
G4double GetY2HalfLength() const
G4double GetPhiTwist() const
G4double GetY1HalfLength() const

References degree, G4GDMLWrite::GenerateName(), G4TwistedTrap::GetAzimuthalAnglePhi(), G4VSolid::GetName(), G4TwistedTrap::GetPhiTwist(), G4TwistedTrap::GetPolarAngleTheta(), G4TwistedTrap::GetTiltAngleAlpha(), G4TwistedTrap::GetX1HalfLength(), G4TwistedTrap::GetX2HalfLength(), G4TwistedTrap::GetX3HalfLength(), G4TwistedTrap::GetX4HalfLength(), G4TwistedTrap::GetY1HalfLength(), G4TwistedTrap::GetY2HalfLength(), G4TwistedTrap::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ TwistedtrdWrite()

void G4GDMLWriteSolids::TwistedtrdWrite ( xercesc::DOMElement *  solElement,
const G4TwistedTrd * const  twistedtrd 
)
protected

Definition at line 986 of file G4GDMLWriteSolids.cc.

988{
989 const G4String& name = GenerateName(twistedtrd->GetName(), twistedtrd);
990
991 xercesc::DOMElement* twistedtrdElement = NewElement("twistedtrd");
992 twistedtrdElement->setAttributeNode(NewAttribute("name", name));
993 twistedtrdElement->setAttributeNode(
994 NewAttribute("x1", 2.0 * twistedtrd->GetX1HalfLength() / mm));
995 twistedtrdElement->setAttributeNode(
996 NewAttribute("x2", 2.0 * twistedtrd->GetX2HalfLength() / mm));
997 twistedtrdElement->setAttributeNode(
998 NewAttribute("y1", 2.0 * twistedtrd->GetY1HalfLength() / mm));
999 twistedtrdElement->setAttributeNode(
1000 NewAttribute("y2", 2.0 * twistedtrd->GetY2HalfLength() / mm));
1001 twistedtrdElement->setAttributeNode(
1002 NewAttribute("z", 2.0 * twistedtrd->GetZHalfLength() / mm));
1003 twistedtrdElement->setAttributeNode(
1004 NewAttribute("PhiTwist", twistedtrd->GetPhiTwist() / degree));
1005 twistedtrdElement->setAttributeNode(NewAttribute("aunit", "deg"));
1006 twistedtrdElement->setAttributeNode(NewAttribute("lunit", "mm"));
1007 solElement->appendChild(twistedtrdElement);
1008}
G4double GetX2HalfLength() const
Definition: G4TwistedTrd.hh:67
G4double GetY2HalfLength() const
Definition: G4TwistedTrd.hh:69
G4double GetY1HalfLength() const
Definition: G4TwistedTrd.hh:68
G4double GetZHalfLength() const
Definition: G4TwistedTrd.hh:70
G4double GetPhiTwist() const
Definition: G4TwistedTrd.hh:71
G4double GetX1HalfLength() const
Definition: G4TwistedTrd.hh:66

References degree, G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4TwistedTrd::GetPhiTwist(), G4TwistedTrd::GetX1HalfLength(), G4TwistedTrd::GetX2HalfLength(), G4TwistedTrd::GetY1HalfLength(), G4TwistedTrd::GetY2HalfLength(), G4TwistedTrd::GetZHalfLength(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ TwistedtubsWrite()

void G4GDMLWriteSolids::TwistedtubsWrite ( xercesc::DOMElement *  solElement,
const G4TwistedTubs * const  twistedtubs 
)
protected

Definition at line 1011 of file G4GDMLWriteSolids.cc.

1013{
1014 const G4String& name = GenerateName(twistedtubs->GetName(), twistedtubs);
1015
1016 xercesc::DOMElement* twistedtubsElement = NewElement("twistedtubs");
1017 twistedtubsElement->setAttributeNode(NewAttribute("name", name));
1018 twistedtubsElement->setAttributeNode(
1019 NewAttribute("twistedangle", twistedtubs->GetPhiTwist() / degree));
1020 twistedtubsElement->setAttributeNode(
1021 NewAttribute("midinnerrad", twistedtubs->GetInnerRadius() / mm));
1022 twistedtubsElement->setAttributeNode(
1023 NewAttribute("midouterrad", twistedtubs->GetOuterRadius() / mm));
1024 twistedtubsElement->setAttributeNode(
1025 NewAttribute("negativeEndz", twistedtubs->GetEndZ(0) / mm));
1026 twistedtubsElement->setAttributeNode(
1027 NewAttribute("positiveEndz", twistedtubs->GetEndZ(1) / mm));
1028 twistedtubsElement->setAttributeNode(
1029 NewAttribute("phi", twistedtubs->GetDPhi() / degree));
1030 twistedtubsElement->setAttributeNode(NewAttribute("aunit", "deg"));
1031 twistedtubsElement->setAttributeNode(NewAttribute("lunit", "mm"));
1032 solElement->appendChild(twistedtubsElement);
1033}
G4double GetOuterRadius() const
G4double GetPhiTwist() const
G4double GetEndZ(G4int i) const
G4double GetInnerRadius() const
G4double GetDPhi() const

References degree, G4GDMLWrite::GenerateName(), G4TwistedTubs::GetDPhi(), G4TwistedTubs::GetEndZ(), G4TwistedTubs::GetInnerRadius(), G4VSolid::GetName(), G4TwistedTubs::GetOuterRadius(), G4TwistedTubs::GetPhiTwist(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by AddSolid().

◆ TWrite()

void G4GDMLWriteMaterials::TWrite ( xercesc::DOMElement *  element,
const G4double T 
)
protectedinherited

Definition at line 83 of file G4GDMLWriteMaterials.cc.

85{
86 xercesc::DOMElement* TElement = NewElement("T");
87 TElement->setAttributeNode(NewAttribute("unit", "K"));
88 TElement->setAttributeNode(NewAttribute("value", T / kelvin));
89 element->appendChild(TElement);
90}
static constexpr double kelvin
Definition: G4SIunits.hh:274

References kelvin, G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by G4GDMLWriteMaterials::MaterialWrite().

◆ UserinfoWrite()

void G4GDMLWrite::UserinfoWrite ( xercesc::DOMElement *  gdmlElement)
virtualinherited

Definition at line 118 of file G4GDMLWrite.cc.

119{
120 if(auxList.size() > 0)
121 {
122#ifdef G4VERBOSE
123 G4cout << "G4GDML: Writing userinfo..." << G4endl;
124#endif
125 userinfoElement = NewElement("userinfo");
126 gdmlElement->appendChild(userinfoElement);
128 }
129}
xercesc::DOMElement * userinfoElement
Definition: G4GDMLWrite.hh:134

References G4GDMLWrite::AddAuxInfo(), G4GDMLWrite::auxList, G4cout, G4endl, G4GDMLWrite::NewElement(), and G4GDMLWrite::userinfoElement.

Referenced by G4GDMLWrite::Write().

◆ VolumeMap()

G4GDMLWrite::VolumeMapType & G4GDMLWrite::VolumeMap ( )
protectedinherited

Definition at line 60 of file G4GDMLWrite.cc.

61{
62 static VolumeMapType instance;
63 return instance;
64}
std::map< const G4LogicalVolume *, G4Transform3D > VolumeMapType
Definition: G4GDMLWrite.hh:53

Referenced by G4GDMLWriteStructure::TraverseVolumeTree(), and G4GDMLWrite::Write().

◆ Write()

G4Transform3D G4GDMLWrite::Write ( const G4String filename,
const G4LogicalVolume *const  topLog,
const G4String schemaPath,
const G4int  depth,
G4bool  storeReferences = true 
)
inherited

Definition at line 203 of file G4GDMLWrite.cc.

207{
208 SchemaLocation = setSchemaLocation;
209 addPointerToName = refs;
210#ifdef G4VERBOSE
211 if(depth == 0)
212 {
213 G4cout << "G4GDML: Writing '" << fname << "'..." << G4endl;
214 }
215 else
216 {
217 G4cout << "G4GDML: Writing module '" << fname << "'..." << G4endl;
218 }
219#endif
221 {
222 G4String ErrorMessage = "File '" + fname + "' already exists!";
223 G4Exception("G4GDMLWrite::Write()", "InvalidSetup", FatalException,
224 ErrorMessage);
225 }
226
227 VolumeMap().clear(); // The module map is global for all modules,
228 // so clear it only at once!
229
230 XMLCh* tempStr = NULL;
231 tempStr = xercesc::XMLString::transcode("LS");
232 xercesc::DOMImplementationRegistry::getDOMImplementation(tempStr);
233 xercesc::XMLString::release(&tempStr);
234 tempStr = xercesc::XMLString::transcode("Range");
235 xercesc::DOMImplementation* impl =
236 xercesc::DOMImplementationRegistry::getDOMImplementation(tempStr);
237 xercesc::XMLString::release(&tempStr);
238 tempStr = xercesc::XMLString::transcode("gdml");
239 doc = impl->createDocument(0, tempStr, 0);
240 xercesc::XMLString::release(&tempStr);
241 xercesc::DOMElement* gdml = doc->getDocumentElement();
242
243#if XERCES_VERSION_MAJOR >= 3
244 // DOM L3 as per Xerces 3.0 API
245 xercesc::DOMLSSerializer* writer =
246 ((xercesc::DOMImplementationLS*) impl)->createLSSerializer();
247
248 xercesc::DOMConfiguration* dc = writer->getDomConfig();
249 dc->setParameter(xercesc::XMLUni::fgDOMWRTFormatPrettyPrint, true);
250
251#else
252
253 xercesc::DOMWriter* writer =
254 ((xercesc::DOMImplementationLS*) impl)->createDOMWriter();
255
256 if(writer->canSetFeature(xercesc::XMLUni::fgDOMWRTFormatPrettyPrint, true))
257 writer->setFeature(xercesc::XMLUni::fgDOMWRTFormatPrettyPrint, true);
258
259#endif
260
261 gdml->setAttributeNode(
262 NewAttribute("xmlns:xsi", "http://www.w3.org/2001/XMLSchema-instance"));
263 gdml->setAttributeNode(
264 NewAttribute("xsi:noNamespaceSchemaLocation", SchemaLocation));
265
266 ExtensionWrite(gdml);
267 DefineWrite(gdml);
268 MaterialsWrite(gdml);
269 SolidsWrite(gdml);
270 StructureWrite(gdml);
271 UserinfoWrite(gdml);
272 SetupWrite(gdml, logvol);
273
274 G4Transform3D R = TraverseVolumeTree(logvol, depth);
275
277 xercesc::XMLFormatTarget* myFormTarget =
278 new xercesc::LocalFileFormatTarget(fname.c_str());
279
280 try
281 {
282#if XERCES_VERSION_MAJOR >= 3
283 // DOM L3 as per Xerces 3.0 API
284 xercesc::DOMLSOutput* theOutput =
285 ((xercesc::DOMImplementationLS*) impl)->createLSOutput();
286 theOutput->setByteStream(myFormTarget);
287 writer->write(doc, theOutput);
288#else
289 writer->writeNode(myFormTarget, *doc);
290#endif
291 } catch(const xercesc::XMLException& toCatch)
292 {
293 char* message = xercesc::XMLString::transcode(toCatch.getMessage());
294 G4cout << "G4GDML: Exception message is: " << message << G4endl;
295 xercesc::XMLString::release(&message);
297 } catch(const xercesc::DOMException& toCatch)
298 {
299 char* message = xercesc::XMLString::transcode(toCatch.msg);
300 G4cout << "G4GDML: Exception message is: " << message << G4endl;
301 xercesc::XMLString::release(&message);
303 } catch(...)
304 {
305 G4cout << "G4GDML: Unexpected Exception!" << G4endl;
307 }
308
309 delete myFormTarget;
310 writer->release();
311
312 if(depth == 0)
313 {
314 G4cout << "G4GDML: Writing '" << fname << "' done !" << G4endl;
315 }
316 else
317 {
318#ifdef G4VERBOSE
319 G4cout << "G4GDML: Writing module '" << fname << "' done !" << G4endl;
320#endif
321 }
322
323 return R;
324}
virtual void SurfacesWrite()=0
virtual void MaterialsWrite(xercesc::DOMElement *)=0
G4bool FileExists(const G4String &) const
Definition: G4GDMLWrite.cc:53
G4String SchemaLocation
Definition: G4GDMLWrite.hh:130
virtual void SetupWrite(xercesc::DOMElement *, const G4LogicalVolume *const)=0
virtual G4Transform3D TraverseVolumeTree(const G4LogicalVolume *const, const G4int)=0
virtual void DefineWrite(xercesc::DOMElement *)=0
virtual void UserinfoWrite(xercesc::DOMElement *)
Definition: G4GDMLWrite.cc:118
virtual void ExtensionWrite(xercesc::DOMElement *)
Definition: G4GDMLWrite.cc:88
virtual void SolidsWrite(xercesc::DOMElement *)=0
VolumeMapType & VolumeMap()
Definition: G4GDMLWrite.cc:60
virtual void StructureWrite(xercesc::DOMElement *)=0
static DLL_API const Transform3D Identity
Definition: Transform3D.h:196

References G4GDMLWrite::addPointerToName, G4GDMLWrite::DefineWrite(), G4GDMLWrite::doc, G4GDMLWrite::ExtensionWrite(), FatalException, G4GDMLWrite::FileExists(), test::fname, G4cout, G4endl, G4Exception(), HepGeom::Transform3D::Identity, G4GDMLWrite::MaterialsWrite(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::overwriteOutputFile, G4GDMLWrite::SchemaLocation, G4GDMLWrite::SetupWrite(), G4GDMLWrite::SolidsWrite(), G4GDMLWrite::StructureWrite(), G4GDMLWrite::SurfacesWrite(), G4GDMLWrite::TraverseVolumeTree(), G4GDMLWrite::UserinfoWrite(), and G4GDMLWrite::VolumeMap().

Referenced by G4GDMLWriteStructure::TraverseVolumeTree().

◆ XtruWrite()

void G4GDMLWriteSolids::XtruWrite ( xercesc::DOMElement *  solElement,
const G4ExtrudedSolid * const  xtru 
)
protected

Definition at line 393 of file G4GDMLWriteSolids.cc.

395{
396 const G4String& name = GenerateName(xtru->GetName(), xtru);
397
398 xercesc::DOMElement* xtruElement = NewElement("xtru");
399 xtruElement->setAttributeNode(NewAttribute("name", name));
400 xtruElement->setAttributeNode(NewAttribute("lunit", "mm"));
401 solElement->appendChild(xtruElement);
402
403 const G4int NumVertex = xtru->GetNofVertices();
404
405 for(G4int i = 0; i < NumVertex; ++i)
406 {
407 xercesc::DOMElement* twoDimVertexElement = NewElement("twoDimVertex");
408 xtruElement->appendChild(twoDimVertexElement);
409
410 const G4TwoVector& vertex = xtru->GetVertex(i);
411
412 twoDimVertexElement->setAttributeNode(NewAttribute("x", vertex.x() / mm));
413 twoDimVertexElement->setAttributeNode(NewAttribute("y", vertex.y() / mm));
414 }
415
416 const G4int NumSection = xtru->GetNofZSections();
417
418 for(G4int i = 0; i < NumSection; ++i)
419 {
420 xercesc::DOMElement* sectionElement = NewElement("section");
421 xtruElement->appendChild(sectionElement);
422
423 const G4ExtrudedSolid::ZSection section = xtru->GetZSection(i);
424
425 sectionElement->setAttributeNode(NewAttribute("zOrder", i));
426 sectionElement->setAttributeNode(
427 NewAttribute("zPosition", section.fZ / mm));
428 sectionElement->setAttributeNode(
429 NewAttribute("xOffset", section.fOffset.x() / mm));
430 sectionElement->setAttributeNode(
431 NewAttribute("yOffset", section.fOffset.y() / mm));
432 sectionElement->setAttributeNode(
433 NewAttribute("scalingFactor", section.fScale));
434 }
435}
double x() const
double y() const
ZSection GetZSection(G4int index) const
G4int GetNofZSections() const
G4int GetNofVertices() const
G4TwoVector GetVertex(G4int index) const

References G4ExtrudedSolid::ZSection::fOffset, G4ExtrudedSolid::ZSection::fScale, G4ExtrudedSolid::ZSection::fZ, G4GDMLWrite::GenerateName(), G4VSolid::GetName(), G4ExtrudedSolid::GetNofVertices(), G4ExtrudedSolid::GetNofZSections(), G4ExtrudedSolid::GetVertex(), G4ExtrudedSolid::GetZSection(), mm, G4InuclParticleNames::name(), G4GDMLWrite::NewAttribute(), G4GDMLWrite::NewElement(), CLHEP::Hep2Vector::x(), and CLHEP::Hep2Vector::y().

Referenced by AddSolid().

◆ ZplaneWrite()

void G4GDMLWriteSolids::ZplaneWrite ( xercesc::DOMElement *  element,
const G4double z,
const G4double rmin,
const G4double rmax 
)
protected

Definition at line 1036 of file G4GDMLWriteSolids.cc.

1039{
1040 xercesc::DOMElement* zplaneElement = NewElement("zplane");
1041 zplaneElement->setAttributeNode(NewAttribute("z", z / mm));
1042 zplaneElement->setAttributeNode(NewAttribute("rmin", rmin / mm));
1043 zplaneElement->setAttributeNode(NewAttribute("rmax", rmax / mm));
1044 element->appendChild(zplaneElement);
1045}

References mm, G4GDMLWrite::NewAttribute(), and G4GDMLWrite::NewElement().

Referenced by G4GDMLWriteParamvol::Polycone_dimensionsWrite(), PolyconeWrite(), G4GDMLWriteParamvol::Polyhedra_dimensionsWrite(), and PolyhedraWrite().

Field Documentation

◆ addPointerToName

G4bool G4GDMLWrite::addPointerToName = true
staticprotectedinherited

◆ auxList

G4GDMLAuxListType G4GDMLWrite::auxList
protectedinherited

Definition at line 136 of file G4GDMLWrite.hh.

Referenced by G4GDMLWrite::AddAuxiliary(), and G4GDMLWrite::UserinfoWrite().

◆ defineElement

xercesc::DOMElement* G4GDMLWriteDefine::defineElement = nullptr
protectedinherited

◆ doc

xercesc::DOMDocument* G4GDMLWrite::doc = nullptr
protectedinherited

◆ elementList

std::vector<const G4Element*> G4GDMLWriteMaterials::elementList
protectedinherited

◆ extElement

xercesc::DOMElement* G4GDMLWrite::extElement = nullptr
protectedinherited

Definition at line 133 of file G4GDMLWrite.hh.

◆ isotopeList

std::vector<const G4Isotope*> G4GDMLWriteMaterials::isotopeList
protectedinherited

◆ kAngularPrecision

const G4double G4GDMLWriteDefine::kAngularPrecision = DBL_EPSILON
staticprotectedinherited

◆ kLinearPrecision

const G4double G4GDMLWriteDefine::kLinearPrecision = DBL_EPSILON
staticprotectedinherited

◆ kRelativePrecision

const G4double G4GDMLWriteDefine::kRelativePrecision = DBL_EPSILON
staticprotectedinherited

◆ materialList

std::vector<const G4Material*> G4GDMLWriteMaterials::materialList
protectedinherited

◆ materialsElement

xercesc::DOMElement* G4GDMLWriteMaterials::materialsElement = nullptr
protectedinherited

◆ maxTransforms

const G4int G4GDMLWriteSolids::maxTransforms = 8
staticprotected

Definition at line 147 of file G4GDMLWriteSolids.hh.

Referenced by BooleanWrite(), and G4GDMLWriteStructure::TraverseVolumeTree().

◆ overwriteOutputFile

G4bool G4GDMLWrite::overwriteOutputFile = false
protectedinherited

Definition at line 137 of file G4GDMLWrite.hh.

Referenced by G4GDMLWrite::SetOutputFileOverwrite(), and G4GDMLWrite::Write().

◆ propertyList

std::vector<const G4PhysicsFreeVector*> G4GDMLWriteMaterials::propertyList
protectedinherited

◆ SchemaLocation

G4String G4GDMLWrite::SchemaLocation
protectedinherited

Definition at line 130 of file G4GDMLWrite.hh.

Referenced by G4GDMLWriteStructure::TraverseVolumeTree(), and G4GDMLWrite::Write().

◆ solidList

std::vector<const G4VSolid*> G4GDMLWriteSolids::solidList
protected

Definition at line 145 of file G4GDMLWriteSolids.hh.

Referenced by AddSolid(), and SolidsWrite().

◆ solidsElement

xercesc::DOMElement* G4GDMLWriteSolids::solidsElement = nullptr
protected

◆ userinfoElement

xercesc::DOMElement* G4GDMLWrite::userinfoElement = nullptr
protectedinherited

Definition at line 134 of file G4GDMLWrite.hh.

Referenced by G4GDMLWrite::UserinfoWrite().


The documentation for this class was generated from the following files: