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 // G4tgrVolumeAssembly implementation 26 // G4tgrVolumeAssembly implementation 27 // 27 // 28 // Author: P.Arce, CIEMAT (November 2007) 28 // Author: P.Arce, CIEMAT (November 2007) 29 // ------------------------------------------- 29 // -------------------------------------------------------------------- 30 30 31 #include "G4tgrVolumeAssembly.hh" 31 #include "G4tgrVolumeAssembly.hh" 32 #include "G4tgrUtils.hh" 32 #include "G4tgrUtils.hh" 33 #include "G4tgrVolumeMgr.hh" 33 #include "G4tgrVolumeMgr.hh" 34 #include "G4tgrPlaceSimple.hh" 34 #include "G4tgrPlaceSimple.hh" 35 #include "G4tgrPlaceDivRep.hh" 35 #include "G4tgrPlaceDivRep.hh" 36 #include "G4tgrPlaceParameterisation.hh" 36 #include "G4tgrPlaceParameterisation.hh" 37 #include "G4tgrFileReader.hh" 37 #include "G4tgrFileReader.hh" 38 #include "G4tgrMessenger.hh" 38 #include "G4tgrMessenger.hh" 39 39 40 // ------------------------------------------- 40 // -------------------------------------------------------------------- 41 G4tgrVolumeAssembly::G4tgrVolumeAssembly() 41 G4tgrVolumeAssembly::G4tgrVolumeAssembly() 42 { 42 { 43 } 43 } 44 44 45 // ------------------------------------------- 45 // -------------------------------------------------------------------- 46 G4tgrVolumeAssembly::~G4tgrVolumeAssembly() 46 G4tgrVolumeAssembly::~G4tgrVolumeAssembly() 47 { 47 { 48 } 48 } 49 49 50 // ------------------------------------------- 50 // -------------------------------------------------------------------- 51 G4tgrVolumeAssembly::G4tgrVolumeAssembly(const 51 G4tgrVolumeAssembly::G4tgrVolumeAssembly(const std::vector<G4String>& wl) 52 { 52 { 53 theType = "VOLAssembly"; 53 theType = "VOLAssembly"; 54 54 55 //---------- set name 55 //---------- set name 56 theName = G4tgrUtils::GetString(wl[1]); 56 theName = G4tgrUtils::GetString(wl[1]); 57 57 58 G4int nVol = G4tgrUtils::GetInt(wl[2]); 58 G4int nVol = G4tgrUtils::GetInt(wl[2]); 59 59 60 G4tgrUtils::CheckWLsize(wl, 3 + nVol * 5, WL 60 G4tgrUtils::CheckWLsize(wl, 3 + nVol * 5, WLSIZE_GE, 61 "G4tgrVolumeAssembly 61 "G4tgrVolumeAssembly::G4tgrVolumeAssembly"); 62 62 63 for(G4int ii = 0; ii < nVol * 5; ii += 5) 63 for(G4int ii = 0; ii < nVol * 5; ii += 5) 64 { 64 { 65 #ifdef G4VERBOSE 65 #ifdef G4VERBOSE 66 if(G4tgrMessenger::GetVerboseLevel() >= 2) 66 if(G4tgrMessenger::GetVerboseLevel() >= 2) 67 { 67 { 68 G4cout << " G4tgrVolumeAssembly::G4tgrVo 68 G4cout << " G4tgrVolumeAssembly::G4tgrVolumeAssembly() -" 69 << " Adding component: " << ii << 69 << " Adding component: " << ii << " - " << wl[ii + 3] << G4endl; 70 } 70 } 71 #endif 71 #endif 72 theComponentNames.push_back(G4tgrUtils::Ge 72 theComponentNames.push_back(G4tgrUtils::GetString(wl[3 + ii + 0])); 73 theComponentRMs.push_back(G4tgrUtils::GetS 73 theComponentRMs.push_back(G4tgrUtils::GetString(wl[3 + ii + 1])); 74 theComponentPos.push_back( 74 theComponentPos.push_back( 75 G4ThreeVector(G4tgrUtils::GetDouble(wl[3 75 G4ThreeVector(G4tgrUtils::GetDouble(wl[3 + ii + 2]), 76 G4tgrUtils::GetDouble(wl[3 76 G4tgrUtils::GetDouble(wl[3 + ii + 3]), 77 G4tgrUtils::GetDouble(wl[3 77 G4tgrUtils::GetDouble(wl[3 + ii + 4]))); 78 } 78 } 79 theVisibility = 1; 79 theVisibility = 1; 80 theRGBColour = new G4double[4]; 80 theRGBColour = new G4double[4]; 81 for(std::size_t ii = 0; ii < 4; ++ii) 81 for(std::size_t ii = 0; ii < 4; ++ii) 82 { 82 { 83 theRGBColour[ii] = -1.; 83 theRGBColour[ii] = -1.; 84 } 84 } 85 85 86 theSolid = nullptr; 86 theSolid = nullptr; 87 87 88 #ifdef G4VERBOSE 88 #ifdef G4VERBOSE 89 if(G4tgrMessenger::GetVerboseLevel() >= 1) 89 if(G4tgrMessenger::GetVerboseLevel() >= 1) 90 { 90 { 91 G4cout << " Created " << *this << G4endl; 91 G4cout << " Created " << *this << G4endl; 92 } 92 } 93 #endif 93 #endif 94 } 94 } 95 95 96 // ------------------------------------------- 96 // -------------------------------------------------------------------- 97 G4tgrPlace* G4tgrVolumeAssembly::AddPlace(cons 97 G4tgrPlace* G4tgrVolumeAssembly::AddPlace(const std::vector<G4String>& wl) 98 { 98 { 99 //---------- Check for exact number of words 99 //---------- Check for exact number of words read 100 G4tgrUtils::CheckWLsize(wl, 7, WLSIZE_EQ, " 100 G4tgrUtils::CheckWLsize(wl, 7, WLSIZE_EQ, " G4tgrVolumeAssembly::AddPlace"); 101 101 102 //---------- set G4tgrPlace 102 //---------- set G4tgrPlace 103 G4tgrPlaceSimple* pl = new G4tgrPlaceSimple( 103 G4tgrPlaceSimple* pl = new G4tgrPlaceSimple(wl); 104 104 105 pl->SetVolume(this); 105 pl->SetVolume(this); 106 thePlacements.push_back(pl); 106 thePlacements.push_back(pl); 107 107 108 #ifdef G4VERBOSE 108 #ifdef G4VERBOSE 109 if(G4tgrMessenger::GetVerboseLevel() >= 2) 109 if(G4tgrMessenger::GetVerboseLevel() >= 2) 110 { 110 { 111 G4cout << " New placement: " << thePlace 111 G4cout << " New placement: " << thePlacements.size() 112 << " added for Volume " << theName 112 << " added for Volume " << theName << " inside " 113 << pl->GetParentName() << " type " 113 << pl->GetParentName() << " type " << pl->GetType() << G4endl; 114 } 114 } 115 #endif 115 #endif 116 //---------- register parent - child 116 //---------- register parent - child 117 G4tgrVolumeMgr::GetInstance()->RegisterParen 117 G4tgrVolumeMgr::GetInstance()->RegisterParentChild(pl->GetParentName(), pl); 118 118 119 return pl; 119 return pl; 120 } 120 } 121 121 122 // ------------------------------------------- 122 // -------------------------------------------------------------------- 123 std::ostream& operator<<(std::ostream& os, con 123 std::ostream& operator<<(std::ostream& os, const G4tgrVolumeAssembly& obj) 124 { 124 { 125 os << "G4tgrVolumeAssembly= " << obj.theName 125 os << "G4tgrVolumeAssembly= " << obj.theName; 126 126 127 for(std::size_t ii = 0; ii < obj.theComponen 127 for(std::size_t ii = 0; ii < obj.theComponentNames.size(); ++ii) 128 { 128 { 129 os << obj.theComponentNames[ii] 129 os << obj.theComponentNames[ii] 130 << " RotMatName= " << obj.theComponentR 130 << " RotMatName= " << obj.theComponentRMs[ii] 131 << " Position= " << obj.theComponentPos 131 << " Position= " << obj.theComponentPos[ii].x() << " " 132 << obj.theComponentPos[ii].y() << " " < 132 << obj.theComponentPos[ii].y() << " " << obj.theComponentPos[ii].z(); 133 } 133 } 134 os << G4endl; 134 os << G4endl; 135 135 136 return os; 136 return os; 137 } 137 } 138 138