ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
PHG4GDMLWriteDefine.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PHG4GDMLWriteDefine.cc
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: PHG4GDMLWriteDefine.cc 68053 2013-03-13 14:39:51Z gcosmo $
28 //
29 // class PHG4GDMLWriteDefine Implementation
30 //
31 // Original author: Zoltan Torzsok, November 2007
32 //
33 // --------------------------------------------------------------------
34 
35 #include "PHG4GDMLWriteDefine.hh"
36 #include <Geant4/G4SystemOfUnits.hh>
37 
38 const G4double PHG4GDMLWriteDefine::kRelativePrecision = DBL_EPSILON;
39 const G4double PHG4GDMLWriteDefine::kAngularPrecision = DBL_EPSILON;
40 const G4double PHG4GDMLWriteDefine::kLinearPrecision = DBL_EPSILON;
41 
42 PHG4GDMLWriteDefine::PHG4GDMLWriteDefine()
43  : PHG4GDMLWrite(), defineElement(0)
44 {
45 }
46 
47 PHG4GDMLWriteDefine::~PHG4GDMLWriteDefine()
48 {
49 }
50 
51 G4ThreeVector PHG4GDMLWriteDefine::GetAngles(const G4RotationMatrix& mtx)
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 
78 void PHG4GDMLWriteDefine::
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 
97 void PHG4GDMLWriteDefine::
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 
114 void PHG4GDMLWriteDefine::
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 PHG4GDMLWriteDefine::DefineWrite(xercesc::DOMElement* element)
132 {
133  std::cout << "PHG4GDML: Writing definitions..." << std::endl;
134 
135  defineElement = NewElement("define");
136  element->appendChild(defineElement);
137 }