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 // G4tgrLineProcessor 27 // 28 // Author: P.Arce, CIEMAT (November 2007) 29 // ------------------------------------------- 30 31 #include "G4tgrLineProcessor.hh" 32 33 #include "G4SystemOfUnits.hh" 34 #include "G4tgrParameterMgr.hh" 35 #include "G4tgrFileIn.hh" 36 #include "G4tgrElementSimple.hh" 37 #include "G4tgrElementFromIsotopes.hh" 38 #include "G4tgrVolume.hh" 39 #include "G4tgrVolumeDivision.hh" 40 #include "G4tgrVolumeAssembly.hh" 41 #include "G4tgrPlaceDivRep.hh" 42 #include "G4tgrPlaceParameterisation.hh" 43 #include "G4tgrVolumeMgr.hh" 44 #include "G4tgrUtils.hh" 45 #include "G4tgrMaterialFactory.hh" 46 #include "G4tgrRotationMatrixFactory.hh" 47 #include "G4tgrMessenger.hh" 48 49 // ------------------------------------------- 50 G4tgrLineProcessor::G4tgrLineProcessor() 51 { 52 volmgr = G4tgrVolumeMgr::GetInstance(); 53 } 54 55 // ------------------------------------------- 56 G4tgrLineProcessor::~G4tgrLineProcessor() 57 { 58 } 59 60 // ------------------------------------------- 61 G4bool G4tgrLineProcessor::ProcessLine(const s 62 { 63 #ifdef G4VERBOSE 64 if(G4tgrMessenger::GetVerboseLevel() >= 1) 65 { 66 G4tgrUtils::DumpVS(wl, "@@@ Processing inp 67 } 68 #endif 69 70 G4String wl0 = wl[0]; 71 for(G4int ii = 0; ii < (G4int)wl0.length(); 72 { 73 wl0[ii] = (char)std::toupper(wl0[ii]); 74 } 75 76 //------------------------------- parameter 77 if(wl0 == ":P") 78 { 79 G4tgrParameterMgr::GetInstance()->AddParam 80 81 //------------------------------- paramete 82 } 83 else if(wl0 == ":PS") 84 { 85 G4tgrParameterMgr::GetInstance()->AddParam 86 87 //------------------------------- isotope 88 } 89 else if(wl0 == ":ISOT") 90 { 91 G4tgrIsotope* isot = G4tgrMaterialFactory: 92 volmgr->RegisterMe(isot); 93 94 //------------------------------- element 95 } 96 else if(wl0 == ":ELEM") 97 { 98 G4tgrElementSimple* elem = 99 G4tgrMaterialFactory::GetInstance()->Add 100 volmgr->RegisterMe(elem); 101 102 //------------------------------- element 103 } 104 else if(wl0 == ":ELEM_FROM_ISOT") 105 { 106 //:ELEM_FROM_ISOT NAME SYMBOL N_ISOT (ISOT 107 G4tgrElementFromIsotopes* elem = 108 G4tgrMaterialFactory::GetInstance()->Add 109 volmgr->RegisterMe(elem); 110 111 //------------------------------- material 112 } 113 else if(wl0 == ":MATE") 114 { 115 G4tgrMaterialSimple* mate = 116 G4tgrMaterialFactory::GetInstance()->Add 117 volmgr->RegisterMe(mate); 118 119 //------------------------------- material 120 } 121 else if((wl0 == ":MIXT") || (wl0 == ":MIXT_B 122 { 123 G4tgrMaterialMixture* mate = 124 G4tgrMaterialFactory::GetInstance()->Add 125 wl, "MaterialMixtureByWeight"); 126 volmgr->RegisterMe(mate); 127 128 //------------------------------- material 129 } 130 else if(wl0 == ":MIXT_BY_NATOMS") 131 { 132 G4tgrMaterialMixture* mate = 133 G4tgrMaterialFactory::GetInstance()->Add 134 wl, "MaterialMixtureByNoAtoms"); 135 volmgr->RegisterMe(mate); 136 137 //------------------------------- material 138 } 139 else if(wl0 == ":MIXT_BY_VOLUME") 140 { 141 G4tgrMaterialMixture* mate = 142 G4tgrMaterialFactory::GetInstance()->Add 143 wl, "MaterialMixtureByVolume"); 144 volmgr->RegisterMe(mate); 145 146 //------------------------------- material 147 // Ionisati 148 } 149 else if(wl0 == ":MATE_MEE") 150 { 151 G4tgrMaterial* mate = G4tgrMaterialFactory 152 G4tgrUtils::GetString(wl[1])); 153 if(mate == 0) 154 { 155 G4Exception("G4tgrLineProcessor::Process 156 FatalException, G4tgrUtils:: 157 return false; 158 } 159 mate->SetIonisationMeanExcitationEnergy(G4 160 161 //------------------------------- material 162 } 163 else if(wl0 == ":MATE_STATE") 164 { 165 G4tgrMaterial* mate = G4tgrMaterialFactory 166 G4tgrUtils::GetString(wl[1])); 167 if(mate == 0) 168 { 169 G4Exception("G4tgrLineProcessor::Process 170 FatalException, wl[1]); 171 } 172 mate->SetState(wl[2]); 173 174 //------------------------------- material 175 } 176 else if(wl0 == ":MATE_TEMPERATURE") 177 { 178 G4tgrMaterial* mate = G4tgrMaterialFactory 179 G4tgrUtils::GetString(wl[1])); 180 if(mate == 0) 181 { 182 G4Exception("G4tgrLineProcessor::Process 183 FatalException, wl[1]); 184 } 185 mate->SetTemperature(G4tgrUtils::GetDouble 186 187 //------------------------------- material 188 } 189 else if(wl0 == ":MATE_PRESSURE") 190 { 191 G4tgrMaterial* mate = G4tgrMaterialFactory 192 G4tgrUtils::GetString(wl[1])); 193 if(mate == 0) 194 { 195 G4Exception("G4tgrLineProcessor::Process 196 FatalException, wl[1]); 197 } 198 mate->SetPressure(G4tgrUtils::GetDouble(wl 199 200 //------------------------------- solid 201 } 202 else if(wl0 == ":SOLID") 203 { // called from here or from G4tgrVolume:: 204 volmgr->CreateSolid(wl, 0); 205 206 //------------------------------- volume 207 } 208 else if(wl0 == ":VOLU") 209 { 210 G4tgrVolume* vol = new G4tgrVolume(wl); 211 volmgr->RegisterMe(vol); 212 213 //--------------------------------- single 214 } 215 else if(wl0 == ":PLACE") 216 { 217 G4tgrVolume* vol = FindVolume(G4tgrUtils:: 218 G4tgrPlace* vpl = vol->AddPlace(wl); 219 volmgr->RegisterMe(vpl); 220 221 //--------------------------------- parame 222 } 223 else if(wl0 == ":PLACE_PARAM") 224 { 225 G4tgrVolume* vol = FindVolu 226 G4tgrPlaceParameterisation* vpl = vol->Add 227 volmgr->RegisterMe(vpl); 228 229 //--------------------------------- divisi 230 } 231 else if((wl0 == ":DIV_NDIV") || (wl0 == ":DI 232 (wl0 == ":DIV_NDIV_WIDTH")) 233 { 234 //---------- Create G4tgrVolumeDivision an 235 G4tgrVolumeDivision* vol = new G4tgrVolume 236 volmgr->RegisterMe(vol); 237 238 //--------------------------------- replic 239 } 240 else if(wl0 == ":REPL") 241 { 242 G4tgrVolume* vol = FindVolume(G4tgrUt 243 G4tgrPlaceDivRep* vpl = vol->AddPlaceRepli 244 volmgr->RegisterMe(vpl); 245 246 //----------------------------- assembly v 247 } 248 else if(wl0 == ":VOLU_ASSEMBLY") 249 { 250 G4tgrVolumeAssembly* vol = new G4tgrVolume 251 volmgr->RegisterMe(vol); 252 253 //----------------------------- assembly v 254 } 255 else if(wl0 == ":PLACE_ASSEMBLY") 256 { 257 G4tgrVolume* vol = FindVolume(G4tgrUtils:: 258 G4tgrPlace* vpl = vol->AddPlace(wl); 259 volmgr->RegisterMe(vpl); 260 261 //--------------------------------- rotat 262 } 263 else if(wl0 == ":ROTM") 264 { 265 //---------- When second word is ':NEXT/:M 266 // rotation matrix that will be 267 G4tgrRotationMatrix* rm = 268 G4tgrRotationMatrixFactory::GetInstance( 269 volmgr->RegisterMe(rm); 270 271 //------------------------------- visualis 272 } 273 else if(wl0 == ":VIS") 274 { 275 std::vector<G4tgrVolume*> vols = 276 volmgr->FindVolumes(G4tgrUtils::GetStrin 277 for(std::size_t ii = 0; ii < vols.size(); 278 { 279 vols[ii]->AddVisibility(wl); 280 } 281 282 //--------------------------------- colour 283 } 284 else if((wl0 == ":COLOUR") || (wl0 == ":COLO 285 { 286 std::vector<G4tgrVolume*> vols = 287 volmgr->FindVolumes(G4tgrUtils::GetStrin 288 for(std::size_t ii = 0; ii < vols.size(); 289 { 290 vols[ii]->AddRGBColour(wl); 291 } 292 293 //--------------------------------- check 294 } 295 else if(wl0 == ":CHECK_OVERLAPS") 296 { 297 std::vector<G4tgrVolume*> vols = 298 volmgr->FindVolumes(G4tgrUtils::GetStrin 299 for(std::size_t ii = 0; ii < vols.size(); 300 { 301 vols[ii]->AddCheckOverlaps(wl); 302 } 303 //--------------------------------- ERROR 304 } 305 else 306 { 307 return false; 308 } 309 310 return true; 311 } 312 313 // ------------------------------------------- 314 G4tgrVolume* G4tgrLineProcessor::FindVolume(co 315 { 316 G4tgrVolume* vol = volmgr->FindVolume(volnam 317 318 if(vol->GetType() == "VOLDivision") 319 { 320 G4Exception("G4tgrLineProcessor::FindVolum 321 FatalException, 322 "Using 'PLACE' for a volume cr 323 } 324 325 return vol; 326 } 327