Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // G4GDMLWriteDefine implementation 27 // 28 // Author: Zoltan Torzsok, November 2007 29 // ------------------------------------------- 30 31 #include "G4GDMLWriteDefine.hh" 32 #include "G4SystemOfUnits.hh" 33 34 const G4double G4GDMLWriteDefine::kRelativePre 35 const G4double G4GDMLWriteDefine::kAngularPrec 36 const G4double G4GDMLWriteDefine::kLinearPreci 37 38 // ------------------------------------------- 39 G4GDMLWriteDefine::G4GDMLWriteDefine() 40 : G4GDMLWrite() 41 { 42 } 43 44 // ------------------------------------------- 45 G4GDMLWriteDefine::~G4GDMLWriteDefine() 46 { 47 } 48 49 // ------------------------------------------- 50 G4ThreeVector G4GDMLWriteDefine::GetAngles(con 51 { 52 G4double x, y, z; 53 G4RotationMatrix mat = mtx; 54 mat.rectify(); // Rectify matrix from possi 55 56 // Direction of rotation given by left-hand 57 58 static const G4double kMatrixPrecision = 10E 59 const G4double cosb = std::sqrt(mtx.xx() * m 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 } 76 77 // ------------------------------------------- 78 void G4GDMLWriteDefine::Scale_vectorWrite(xerc 79 cons 80 cons 81 cons 82 { 83 const G4double x = 84 (std::fabs(scl.x() - 1.0) < kRelativePreci 85 const G4double y = 86 (std::fabs(scl.y() - 1.0) < kRelativePreci 87 const G4double z = 88 (std::fabs(scl.z() - 1.0) < kRelativePreci 89 90 xercesc::DOMElement* scaleElement = NewEleme 91 scaleElement->setAttributeNode(NewAttribute( 92 scaleElement->setAttributeNode(NewAttribute( 93 scaleElement->setAttributeNode(NewAttribute( 94 scaleElement->setAttributeNode(NewAttribute( 95 element->appendChild(scaleElement); 96 } 97 98 // ------------------------------------------- 99 void G4GDMLWriteDefine::Rotation_vectorWrite(x 100 c 101 c 102 c 103 { 104 const G4double x = (std::fabs(rot.x()) < kAn 105 const G4double y = (std::fabs(rot.y()) < kAn 106 const G4double z = (std::fabs(rot.z()) < kAn 107 108 xercesc::DOMElement* rotationElement = NewEl 109 rotationElement->setAttributeNode(NewAttribu 110 rotationElement->setAttributeNode(NewAttribu 111 rotationElement->setAttributeNode(NewAttribu 112 rotationElement->setAttributeNode(NewAttribu 113 rotationElement->setAttributeNode(NewAttribu 114 element->appendChild(rotationElement); 115 } 116 117 // ------------------------------------------- 118 void G4GDMLWriteDefine::Position_vectorWrite(x 119 c 120 c 121 c 122 { 123 const G4double x = (std::fabs(pos.x()) < kLi 124 const G4double y = (std::fabs(pos.y()) < kLi 125 const G4double z = (std::fabs(pos.z()) < kLi 126 127 xercesc::DOMElement* positionElement = NewEl 128 positionElement->setAttributeNode(NewAttribu 129 positionElement->setAttributeNode(NewAttribu 130 positionElement->setAttributeNode(NewAttribu 131 positionElement->setAttributeNode(NewAttribu 132 positionElement->setAttributeNode(NewAttribu 133 element->appendChild(positionElement); 134 } 135 136 // ------------------------------------------- 137 void G4GDMLWriteDefine::DefineWrite(xercesc::D 138 { 139 #ifdef G4VERBOSE 140 G4cout << "G4GDML: Writing definitions..." < 141 #endif 142 defineElement = NewElement("define"); 143 element->appendChild(defineElement); 144 } 145