Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4GDMLWriteDefine.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 // $Id: G4GDMLWriteDefine.cc 93077 2015-10-02 14:38:25Z gcosmo $
28 //
29 // class G4GDMLWriteDefine Implementation
30 //
31 // Original author: Zoltan Torzsok, November 2007
32 //
33 // --------------------------------------------------------------------
34 
35 #include "G4GDMLWriteDefine.hh"
36 #include "G4SystemOfUnits.hh"
37 
41 
43  : G4GDMLWrite(), defineElement(0)
44 {
45 }
46 
48 {
49 }
50 
52 {
53  G4double x,y,z;
54  G4RotationMatrix mat = mtx;
55  mat.rectify(); // Rectify matrix from possible roundoff errors
56 
57  // Direction of rotation given by left-hand rule; clockwise rotation
58 
59  static const G4double kMatrixPrecision = 10E-10;
60  const G4double cosb = std::sqrt(mtx.xx()*mtx.xx()+mtx.yx()*mtx.yx());
61 
62  if (cosb > kMatrixPrecision)
63  {
64  x = std::atan2(mtx.zy(),mtx.zz());
65  y = std::atan2(-mtx.zx(),cosb);
66  z = std::atan2(mtx.yx(),mtx.xx());
67  }
68  else
69  {
70  x = std::atan2(-mtx.yz(),mtx.yy());
71  y = std::atan2(-mtx.zx(),cosb);
72  z = 0.0;
73  }
74 
75  return G4ThreeVector(x,y,z);
76 }
77 
79 Scale_vectorWrite(xercesc::DOMElement* element, const G4String& tag,
80  const G4String& name, const G4ThreeVector& scl)
81 {
82  const G4double x = (std::fabs(scl.x()-1.0) < kRelativePrecision)
83  ? 1.0 : scl.x();
84  const G4double y = (std::fabs(scl.y()-1.0) < kRelativePrecision)
85  ? 1.0 : scl.y();
86  const G4double z = (std::fabs(scl.z()-1.0) < kRelativePrecision)
87  ? 1.0 : scl.z();
88 
89  xercesc::DOMElement* scaleElement = NewElement(tag);
90  scaleElement->setAttributeNode(NewAttribute("name",name));
91  scaleElement->setAttributeNode(NewAttribute("x",x));
92  scaleElement->setAttributeNode(NewAttribute("y",y));
93  scaleElement->setAttributeNode(NewAttribute("z",z));
94  element->appendChild(scaleElement);
95 }
96 
98 Rotation_vectorWrite(xercesc::DOMElement* element, const G4String& tag,
99  const G4String& name, const G4ThreeVector& rot)
100 {
101  const G4double x = (std::fabs(rot.x()) < kAngularPrecision) ? 0.0 : rot.x();
102  const G4double y = (std::fabs(rot.y()) < kAngularPrecision) ? 0.0 : rot.y();
103  const G4double z = (std::fabs(rot.z()) < kAngularPrecision) ? 0.0 : rot.z();
104 
105  xercesc::DOMElement* rotationElement = NewElement(tag);
106  rotationElement->setAttributeNode(NewAttribute("name",name));
107  rotationElement->setAttributeNode(NewAttribute("x",x/degree));
108  rotationElement->setAttributeNode(NewAttribute("y",y/degree));
109  rotationElement->setAttributeNode(NewAttribute("z",z/degree));
110  rotationElement->setAttributeNode(NewAttribute("unit","deg"));
111  element->appendChild(rotationElement);
112 }
113 
115 Position_vectorWrite(xercesc::DOMElement* element, const G4String& tag,
116  const G4String& name, const G4ThreeVector& pos)
117 {
118  const G4double x = (std::fabs(pos.x()) < kLinearPrecision) ? 0.0 : pos.x();
119  const G4double y = (std::fabs(pos.y()) < kLinearPrecision) ? 0.0 : pos.y();
120  const G4double z = (std::fabs(pos.z()) < kLinearPrecision) ? 0.0 : pos.z();
121 
122  xercesc::DOMElement* positionElement = NewElement(tag);
123  positionElement->setAttributeNode(NewAttribute("name",name));
124  positionElement->setAttributeNode(NewAttribute("x",x/mm));
125  positionElement->setAttributeNode(NewAttribute("y",y/mm));
126  positionElement->setAttributeNode(NewAttribute("z",z/mm));
127  positionElement->setAttributeNode(NewAttribute("unit","mm"));
128  element->appendChild(positionElement);
129 }
130 
131 void G4GDMLWriteDefine::DefineWrite(xercesc::DOMElement* element)
132 {
133  G4cout << "G4GDML: Writing definitions..." << G4endl;
134 
135  defineElement = NewElement("define");
136  element->appendChild(defineElement);
137 }
static const G4double kAngularPrecision
void Scale_vectorWrite(xercesc::DOMElement *, const G4String &, const G4String &, const G4ThreeVector &)
const XML_Char * name
Definition: expat.h:151
double xx() const
static constexpr double mm
Definition: G4SIunits.hh:115
CLHEP::Hep3Vector G4ThreeVector
double x() const
Definition: xmlparse.cc:187
double yy() const
xercesc::DOMElement * defineElement
tuple x
Definition: test.py:50
xercesc::DOMElement * NewElement(const G4String &)
Definition: G4GDMLWrite.cc:160
double zx() const
double z() const
void Position_vectorWrite(xercesc::DOMElement *, const G4String &, const G4String &, const G4ThreeVector &)
double yz() const
virtual void DefineWrite(xercesc::DOMElement *)
G4GLOB_DLL std::ostream G4cout
static constexpr double degree
Definition: G4SIunits.hh:144
static const G4double kRelativePrecision
#define DBL_EPSILON
Definition: templates.hh:87
xercesc::DOMAttr * NewAttribute(const G4String &, const G4String &)
Definition: G4GDMLWrite.cc:136
static const G4double kLinearPrecision
double zy() const
G4ThreeVector GetAngles(const G4RotationMatrix &)
double y() const
void Rotation_vectorWrite(xercesc::DOMElement *, const G4String &, const G4String &, const G4ThreeVector &)
tuple z
Definition: test.py:28
#define G4endl
Definition: G4ios.hh:61
double yx() const
double G4double
Definition: G4Types.hh:76
double zz() const
virtual ~G4GDMLWriteDefine()
static const G4double pos