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 // 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(const std::vector<G4String>& wl) 40 { 41 // :SOLID/:VOLU SOL MULTIUNION NSOL SOL1 ROTM1 POSX1 POSY1 POSZ1 42 // SOL2 ROTM2 POSX2 POSY2 POSZ2 ... 43 44 //---------- Set name 45 theName = G4tgrUtils::GetString(wl[1]); 46 nSolid = G4tgrUtils::GetInt(wl[3]); 47 G4tgrVolumeMgr* volmgr = G4tgrVolumeMgr::GetInstance(); 48 const G4tgrSolid* sol1; 49 50 if(G4int(wl.size()) != 4 + 5*nSolid) 51 { 52 G4String Err1 = "Solid type MULTIUNION with "; 53 G4String Err2 = std::to_string(nSolid); 54 G4String Err3 = " Solids, should have "; 55 G4String Err4 = std::to_string(4 + 5*nSolid); 56 G4String Err5 = " parameters."; 57 G4String Err6 = " It has "+G4UIcommand::ConvertToString(int(wl.size())); 58 59 G4String ErrMessage = Err1 + Err2 + Err3 + Err4 + Err5 + Err6; 60 61 G4tgrUtils::DumpVS(wl, "G4tgrSolidMultiUnion::G4tgrSolidMultiUnion()"); 62 G4Exception("G4tgrSolidMultiUnion::G4tgrSolidMultiUnion()", "InvalidInput", 63 FatalException, ErrMessage); 64 } 65 66 for (G4int i = 4; i< nSolid*5; i+=5) 67 { 68 sol1 = volmgr->FindSolid(G4tgrUtils::GetString(wl[i])); 69 if(sol1 == nullptr) 70 { 71 sol1 = volmgr->FindVolume(G4tgrUtils::GetString(wl[i]), 1)->GetSolid(); 72 } 73 theSolids.push_back(sol1); 74 75 //---------- Populate transformation matrix 76 theRotMatName = G4tgrUtils::GetString(wl[i + 1]); 77 theRotMat = 78 G4tgbRotationMatrixMgr::GetInstance()->FindOrBuildG4RotMatrix(theRotMatName); 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(); ++ii) 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(this); 103 } 104 105 // -------------------------------------------------------------------- 106 G4tgrSolidMultiUnion::~G4tgrSolidMultiUnion() 107 { 108 } 109 110 // -------------------------------------------------------------------- 111 std::ostream& operator<<(std::ostream& os, const G4tgrSolidMultiUnion& sol) 112 { 113 os << "G4tgrSolidMultiUnion= " << sol.theName << " of type " << sol.theType 114 << " PARAMS: "; 115 if(sol.theSolidParams.size() != 0) 116 { 117 std::vector<G4double> solpar = *(sol.theSolidParams[0]); 118 for(std::size_t ii = 0; ii < solpar.size(); ++ii) 119 { 120 os << solpar[ii] << " "; 121 } 122 } 123 os << G4endl; 124 125 return os; 126 } 127