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