| Geant4 Cross Reference |
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 // G4tgrSolidBoolean implementation
27 //
28 // Author: P.Arce, CIEMAT (November 2007)
29 // --------------------------------------------------------------------
30
31 #include "G4tgrSolidBoolean.hh"
32 #include "G4tgrUtils.hh"
33 #include "G4tgrVolume.hh"
34 #include "G4tgrVolumeMgr.hh"
35 #include "G4tgrMessenger.hh"
36 #include "G4tgrFileReader.hh"
37
38 // --------------------------------------------------------------------
39 G4tgrSolidBoolean::G4tgrSolidBoolean(const std::vector<G4String>& wl)
40 {
41 // :SOLID/:VOLU VOLU UNION/SUBS/INTERS VOLU1 VOLU2 ROTM POSX POSY POSZ
42
43 if(wl.size() != 9)
44 {
45 G4tgrUtils::DumpVS(wl, "G4tgrSolidBoolean::G4tgrSolidBoolean()");
46 G4Exception("G4tgrSolidBoolean::G4tgrSolidBoolean()", "InvalidInput",
47 FatalException, "Line read with less or more than 9 words.");
48 }
49
50 //---------- Set name
51 theName = G4tgrUtils::GetString(wl[1]);
52
53 G4tgrVolumeMgr* volmgr = G4tgrVolumeMgr::GetInstance();
54 const G4tgrSolid* sol1 = volmgr->FindSolid(G4tgrUtils::GetString(wl[3]));
55 if(sol1 == nullptr)
56 {
57 sol1 = volmgr->FindVolume(G4tgrUtils::GetString(wl[3]), 1)->GetSolid();
58 }
59 const G4tgrSolid* sol2 = volmgr->FindSolid(G4tgrUtils::GetString(wl[4]));
60 if(sol2 == nullptr)
61 {
62 sol2 = volmgr->FindVolume(G4tgrUtils::GetString(wl[4]), 1)->GetSolid();
63 }
64 theSolids.push_back(sol1);
65 theSolids.push_back(sol2);
66
67 //---------- Set relative placement and rotation matrix
68 theRelativeRotMatName = G4tgrUtils::GetString(wl[5]);
69 theRelativePlace =
70 G4ThreeVector(G4tgrUtils::GetDouble(wl[6]), G4tgrUtils::GetDouble(wl[7]),
71 G4tgrUtils::GetDouble(wl[8]));
72 //---------- Set solid type
73 G4String wl2 = wl[2];
74 for(G4int ii = 0; ii < (G4int)wl2.length(); ++ii)
75 {
76 wl2[ii] = (char)std::toupper(wl2[ii]);
77 }
78 theType = "Boolean_" + wl2;
79
80 #ifdef G4VERBOSE
81 if(G4tgrMessenger::GetVerboseLevel() >= 1)
82 {
83 G4cout << " Created " << *this << G4endl;
84 }
85 #endif
86
87 G4tgrVolumeMgr::GetInstance()->RegisterMe(this);
88 }
89
90 // --------------------------------------------------------------------
91 G4tgrSolidBoolean::~G4tgrSolidBoolean()
92 {
93 }
94
95 // --------------------------------------------------------------------
96 const G4String& G4tgrSolidBoolean::GetRelativeRotMatName() const
97 {
98 return theRelativeRotMatName;
99 }
100
101 // --------------------------------------------------------------------
102 G4ThreeVector G4tgrSolidBoolean::GetRelativePlace() const
103 {
104 return theRelativePlace;
105 }
106
107 // --------------------------------------------------------------------
108 std::ostream& operator<<(std::ostream& os, const G4tgrSolidBoolean& sol)
109 {
110 os << "G4tgrSolidBoolean= " << sol.theName << " of type " << sol.theType
111 << " PARAMS: ";
112 if(sol.theSolidParams.size() != 0)
113 {
114 std::vector<G4double> solpar = *(sol.theSolidParams[0]);
115 for(std::size_t ii = 0; ii < solpar.size(); ++ii)
116 {
117 os << solpar[ii] << " ";
118 }
119 }
120 os << G4endl;
121
122 return os;
123 }
124