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 // G4tgrSolidMultiUnion implementation 27 // 28 // Author: P.Heidary, AEOI - November 2021 29 // ------------------------------------------- 30 31 #include "G4tgrSolidMultiUnion.hh" 32 #include "G4tgrUtils.hh" 33 #include "G4tgrVolume.hh" 34 #include "G4tgrVolumeMgr.hh" 35 #include "G4tgrMessenger.hh" 36 #include "G4tgbRotationMatrixMgr.hh" 37 38 // ------------------------------------------- 39 G4tgrSolidMultiUnion::G4tgrSolidMultiUnion(con 40 { 41 // :SOLID/:VOLU SOL MULTIUNION NSOL SOL1 ROT 42 // SOL2 ROTM2 P 43 44 //---------- Set name 45 theName = G4tgrUtils::GetString(wl[1]); 46 nSolid = G4tgrUtils::GetInt(wl[3]); 47 G4tgrVolumeMgr* volmgr = G4tgrVolumeMgr::Get 48 const G4tgrSolid* sol1; 49 50 if(G4int(wl.size()) != 4 + 5*nSolid) 51 { 52 G4String Err1 = "Solid type MULTIUNION wit 53 G4String Err2 = std::to_string(nSolid); 54 G4String Err3 = " Solids, should have "; 55 G4String Err4 = std::to_string(4 + 5*nSoli 56 G4String Err5 = " parameters."; 57 G4String Err6 = " It has "+G4UIcommand::Co 58 59 G4String ErrMessage = Err1 + Err2 + Err3 + 60 61 G4tgrUtils::DumpVS(wl, "G4tgrSolidMultiUni 62 G4Exception("G4tgrSolidMultiUnion::G4tgrSo 63 FatalException, ErrMessage); 64 } 65 66 for (G4int i = 4; i< nSolid*5; i+=5) 67 { 68 sol1 = volmgr->FindSolid(G4tgrUtils::GetSt 69 if(sol1 == nullptr) 70 { 71 sol1 = volmgr->FindVolume(G4tgrUtils::Ge 72 } 73 theSolids.push_back(sol1); 74 75 //---------- Populate transformation matri 76 theRotMatName = G4tgrUtils::GetString(wl[i 77 theRotMat = 78 G4tgbRotationMatrixMgr::GetInstance()->Fin 79 thePosition = 80 G4ThreeVector(G4tgrUtils::GetDouble(wl[i+2 81 G4tgrUtils::GetDouble(wl[i+3 82 G4tgrUtils::GetDouble(wl[i+4 83 tr1 = G4Transform3D(*theRotMat, thePosition) 84 theTransformations.push_back(tr1); 85 } 86 87 // ---------- Set solid type 88 G4String wl2 = wl[2]; 89 for(G4int ii = 0; ii < (G4int)wl2.length(); 90 { 91 wl2[ii] = (char)std::toupper(wl2[ii]); 92 } 93 theType = std::move(wl2); 94 95 #ifdef G4VERBOSE 96 if(G4tgrMessenger::GetVerboseLevel() >= 1) 97 { 98 G4cout << " Created " << *this << G4endl; 99 } 100 #endif 101 102 G4tgrVolumeMgr::GetInstance()->RegisterMe(th 103 } 104 105 // ------------------------------------------- 106 G4tgrSolidMultiUnion::~G4tgrSolidMultiUnion() 107 { 108 } 109 110 // ------------------------------------------- 111 std::ostream& operator<<(std::ostream& os, con 112 { 113 os << "G4tgrSolidMultiUnion= " << sol.theNam 114 << " PARAMS: "; 115 if(sol.theSolidParams.size() != 0) 116 { 117 std::vector<G4double> solpar = *(sol.theSo 118 for(std::size_t ii = 0; ii < solpar.size() 119 { 120 os << solpar[ii] << " "; 121 } 122 } 123 os << G4endl; 124 125 return os; 126 } 127