| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // G4tgrSolidBoolean implementation <<
27 // 26 //
28 // Author: P.Arce, CIEMAT (November 2007) << 27 // $Id: G4tgrSolidBoolean.cc 66363 2012-12-18 09:12:54Z gcosmo $
29 // ------------------------------------------- << 28 //
>> 29 //
>> 30 // class G4tgrSolidBoolean
>> 31
>> 32 // History:
>> 33 // - Created. P.Arce, CIEMAT (November 2007)
>> 34 // -------------------------------------------------------------------------
30 35
31 #include "G4tgrSolidBoolean.hh" 36 #include "G4tgrSolidBoolean.hh"
32 #include "G4tgrUtils.hh" 37 #include "G4tgrUtils.hh"
33 #include "G4tgrVolume.hh" 38 #include "G4tgrVolume.hh"
34 #include "G4tgrVolumeMgr.hh" 39 #include "G4tgrVolumeMgr.hh"
35 #include "G4tgrMessenger.hh" 40 #include "G4tgrMessenger.hh"
36 #include "G4tgrFileReader.hh" 41 #include "G4tgrFileReader.hh"
37 42
38 // ------------------------------------------- << 43 // -------------------------------------------------------------------------
39 G4tgrSolidBoolean::G4tgrSolidBoolean(const std << 44 G4tgrSolidBoolean::G4tgrSolidBoolean( const std::vector<G4String>& wl )
40 { << 45 {
41 // :SOLID/:VOLU VOLU UNION/SUBS/INTERS VOLU1 46 // :SOLID/:VOLU VOLU UNION/SUBS/INTERS VOLU1 VOLU2 ROTM POSX POSY POSZ
42 47
43 if(wl.size() != 9) << 48 if( wl.size() != 9 )
44 { 49 {
45 G4tgrUtils::DumpVS(wl, "G4tgrSolidBoolean: 50 G4tgrUtils::DumpVS(wl, "G4tgrSolidBoolean::G4tgrSolidBoolean()");
46 G4Exception("G4tgrSolidBoolean::G4tgrSolid 51 G4Exception("G4tgrSolidBoolean::G4tgrSolidBoolean()", "InvalidInput",
47 FatalException, "Line read wit 52 FatalException, "Line read with less or more than 9 words.");
48 } 53 }
49 54
50 //---------- Set name << 55 //---------- Set name
51 theName = G4tgrUtils::GetString(wl[1]); << 56 theName = G4tgrUtils::GetString( wl[1] );
52 57
53 G4tgrVolumeMgr* volmgr = G4tgrVolumeMgr::Get 58 G4tgrVolumeMgr* volmgr = G4tgrVolumeMgr::GetInstance();
54 const G4tgrSolid* sol1 = volmgr->FindSolid(G << 59 const G4tgrSolid* sol1 = volmgr->FindSolid( G4tgrUtils::GetString( wl[3] ));
55 if(sol1 == nullptr) << 60 if( !sol1 )
56 { 61 {
57 sol1 = volmgr->FindVolume(G4tgrUtils::GetS << 62 sol1 = volmgr->FindVolume( G4tgrUtils::GetString( wl[3] ), 1)->GetSolid();
58 } 63 }
59 const G4tgrSolid* sol2 = volmgr->FindSolid(G << 64 const G4tgrSolid* sol2 = volmgr->FindSolid( G4tgrUtils::GetString( wl[4] ));
60 if(sol2 == nullptr) << 65 if( !sol2 )
61 { 66 {
62 sol2 = volmgr->FindVolume(G4tgrUtils::GetS << 67 sol2 = volmgr->FindVolume( G4tgrUtils::GetString( wl[4] ), 1)->GetSolid();
63 } 68 }
64 theSolids.push_back(sol1); << 69 theSolids.push_back( sol1 );
65 theSolids.push_back(sol2); << 70 theSolids.push_back( sol2 );
66 71
67 //---------- Set relative placement and rota 72 //---------- Set relative placement and rotation matrix
68 theRelativeRotMatName = G4tgrUtils::GetStrin << 73 theRelativeRotMatName = G4tgrUtils::GetString( wl[5] );
69 theRelativePlace = << 74 theRelativePlace = G4ThreeVector( G4tgrUtils::GetDouble(wl[6]),
70 G4ThreeVector(G4tgrUtils::GetDouble(wl[6]) << 75 G4tgrUtils::GetDouble(wl[7]),
71 G4tgrUtils::GetDouble(wl[8]) << 76 G4tgrUtils::GetDouble(wl[8]) );
72 //---------- Set solid type 77 //---------- Set solid type
73 G4String wl2 = wl[2]; 78 G4String wl2 = wl[2];
74 for(G4int ii = 0; ii < (G4int)wl2.length(); << 79 for( size_t ii = 0; ii < wl2.length(); ii++ )
75 { 80 {
76 wl2[ii] = (char)std::toupper(wl2[ii]); << 81 wl2[ii] = toupper( wl2[ii] );
77 } 82 }
78 theType = "Boolean_" + wl2; 83 theType = "Boolean_" + wl2;
79 84
80 #ifdef G4VERBOSE 85 #ifdef G4VERBOSE
81 if(G4tgrMessenger::GetVerboseLevel() >= 1) << 86 if( G4tgrMessenger::GetVerboseLevel() >= 1 )
82 { 87 {
83 G4cout << " Created " << *this << G4endl; << 88 G4cout << " Created " << *this << G4endl;
84 } 89 }
85 #endif 90 #endif
86 91
87 G4tgrVolumeMgr::GetInstance()->RegisterMe(th << 92 G4tgrVolumeMgr::GetInstance()->RegisterMe( this );
88 } 93 }
89 94
90 // ------------------------------------------- << 95
>> 96 // -------------------------------------------------------------------------
91 G4tgrSolidBoolean::~G4tgrSolidBoolean() 97 G4tgrSolidBoolean::~G4tgrSolidBoolean()
92 { 98 {
93 } 99 }
94 100
95 // ------------------------------------------- << 101
>> 102 // -------------------------------------------------------------------------
96 const G4String& G4tgrSolidBoolean::GetRelative 103 const G4String& G4tgrSolidBoolean::GetRelativeRotMatName() const
97 { 104 {
98 return theRelativeRotMatName; 105 return theRelativeRotMatName;
99 } 106 }
100 107
101 // ------------------------------------------- << 108
>> 109 // -------------------------------------------------------------------------
102 G4ThreeVector G4tgrSolidBoolean::GetRelativePl 110 G4ThreeVector G4tgrSolidBoolean::GetRelativePlace() const
103 { 111 {
104 return theRelativePlace; 112 return theRelativePlace;
105 } 113 }
106 114
107 // ------------------------------------------- << 115
>> 116 // -------------------------------------------------------------------------
108 std::ostream& operator<<(std::ostream& os, con 117 std::ostream& operator<<(std::ostream& os, const G4tgrSolidBoolean& sol)
109 { 118 {
110 os << "G4tgrSolidBoolean= " << sol.theName < 119 os << "G4tgrSolidBoolean= " << sol.theName << " of type " << sol.theType
111 << " PARAMS: "; 120 << " PARAMS: ";
112 if(sol.theSolidParams.size() != 0) << 121 if( sol.theSolidParams.size() != 0 )
113 { 122 {
114 std::vector<G4double> solpar = *(sol.theSo << 123 std::vector<G4double> solpar = *(sol.theSolidParams[0]);
115 for(std::size_t ii = 0; ii < solpar.size() << 124 for( size_t ii = 0; ii < solpar.size(); ii++)
116 { 125 {
117 os << solpar[ii] << " "; << 126 os << solpar[ii] << " " ;
118 } 127 }
119 } 128 }
120 os << G4endl; 129 os << G4endl;
121 130
122 return os; 131 return os;
123 } 132 }
124 133