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 // Gorad (Geant4 Open-source Radiation Analys 26 // Gorad (Geant4 Open-source Radiation Analysis and Design) 27 // 27 // 28 // Author : Makoto Asai (SLAC National Accele 28 // Author : Makoto Asai (SLAC National Accelerator Laboratory) 29 // 29 // 30 // Development of Gorad is funded by NASA Joh 30 // Development of Gorad is funded by NASA Johnson Space Center (JSC) 31 // under the contract NNJ15HK11B. 31 // under the contract NNJ15HK11B. 32 // 32 // 33 // ******************************************* 33 // ******************************************************************** 34 // 34 // 35 // GRPrimGenActionMessenger.hh 35 // GRPrimGenActionMessenger.hh 36 // Header file of a messenger that handles p 36 // Header file of a messenger that handles primary generator action. 37 // Input radiation spectrum file should be i 37 // Input radiation spectrum file should be in ASCII format and each 38 // row should have low-end kinetic energy an 38 // row should have low-end kinetic energy and differential flux 39 // separated by a space. 39 // separated by a space. 40 // 40 // 41 // History 41 // History 42 // September 8th, 2020 : first implementatio 42 // September 8th, 2020 : first implementation 43 // 43 // 44 // ******************************************* 44 // ******************************************************************** 45 45 46 #include "GRPrimGenActionMessenger.hh" 46 #include "GRPrimGenActionMessenger.hh" 47 47 48 #include "G4UIdirectory.hh" 48 #include "G4UIdirectory.hh" 49 #include "G4UIcommand.hh" 49 #include "G4UIcommand.hh" 50 #include "G4UIparameter.hh" 50 #include "G4UIparameter.hh" 51 51 52 GRPrimGenActionMessenger::GRPrimGenActionMesse 52 GRPrimGenActionMessenger::GRPrimGenActionMessenger(G4int verboseLvl) 53 : verboseLevel(verboseLvl) 53 : verboseLevel(verboseLvl) 54 { 54 { 55 G4UIparameter* param = nullptr; 55 G4UIparameter* param = nullptr; 56 56 57 srcDir = new G4UIdirectory("/gorad/source/") 57 srcDir = new G4UIdirectory("/gorad/source/"); 58 srcDir->SetGuidance("Define primary particle 58 srcDir->SetGuidance("Define primary particle spectrum."); 59 59 60 defCmd = new G4UIcommand("/gorad/source/defi 60 defCmd = new G4UIcommand("/gorad/source/define",this); 61 defCmd->SetGuidance("Define primary particle 61 defCmd->SetGuidance("Define primary particle spectrum."); 62 defCmd->SetGuidance("[usage] /gorad/source/d 62 defCmd->SetGuidance("[usage] /gorad/source/define pName fName srcType radius unit (x0 y0 z0)"); 63 defCmd->SetGuidance(" pName : Particle name 63 defCmd->SetGuidance(" pName : Particle name."); 64 defCmd->SetGuidance(" fName : File name of 64 defCmd->SetGuidance(" fName : File name of the spectrum. Directory may be preceded to the file name."); 65 defCmd->SetGuidance(" srcType : defines how 65 defCmd->SetGuidance(" srcType : defines how the primaries are generated."); 66 defCmd->SetGuidance(" Arb : Gener 66 defCmd->SetGuidance(" Arb : Generate primaries along with the defined spectrum."); 67 defCmd->SetGuidance(" LW : Gener 67 defCmd->SetGuidance(" LW : Generate primaries in flat ditribution with track weight representing the spectrum."); 68 defCmd->SetGuidance(" radius : Radius of th 68 defCmd->SetGuidance(" radius : Radius of the source sphere. If it is set to -1 (default), radius is set to 98% of the world volume."); 69 defCmd->SetGuidance(" unit : Unit of radius 69 defCmd->SetGuidance(" unit : Unit of radius value."); 70 defCmd->SetGuidance(" x0,y0,z0 : (Optional) 70 defCmd->SetGuidance(" x0,y0,z0 : (Optional) location of the center of the sphere (same unit is applied as radius)."); 71 defCmd->SetGuidance(" By default 71 defCmd->SetGuidance(" By default the sphere is located at the center of the world volume."); 72 defCmd->SetGuidance("[note] This command is 72 defCmd->SetGuidance("[note] This command is not mandatory but dwialternatively granular /gps/ commands can be used for defining primary particles."); 73 particlePara = new G4UIparameter("pName",'s' 73 particlePara = new G4UIparameter("pName",'s',false); 74 defCmd->SetParameter(particlePara); 74 defCmd->SetParameter(particlePara); 75 param = new G4UIparameter("fName",'s',false) 75 param = new G4UIparameter("fName",'s',false); 76 defCmd->SetParameter(param); 76 defCmd->SetParameter(param); 77 param = new G4UIparameter("srcType",'s',fals 77 param = new G4UIparameter("srcType",'s',false); 78 param->SetParameterCandidates("Arb LW"); 78 param->SetParameterCandidates("Arb LW"); 79 defCmd->SetParameter(param); 79 defCmd->SetParameter(param); 80 param = new G4UIparameter("radius",'d',true) 80 param = new G4UIparameter("radius",'d',true); 81 param->SetParameterRange("radius == -1. || r 81 param->SetParameterRange("radius == -1. || radius>0."); 82 param->SetDefaultValue(-1.); 82 param->SetDefaultValue(-1.); 83 defCmd->SetParameter(param); 83 defCmd->SetParameter(param); 84 param = new G4UIparameter("unit",'s',true); 84 param = new G4UIparameter("unit",'s',true); 85 param->SetDefaultUnit("m"); 85 param->SetDefaultUnit("m"); 86 defCmd->SetParameter(param); 86 defCmd->SetParameter(param); 87 param = new G4UIparameter("x0",'d',true); 87 param = new G4UIparameter("x0",'d',true); 88 param->SetDefaultValue(0.); 88 param->SetDefaultValue(0.); 89 defCmd->SetParameter(param); 89 defCmd->SetParameter(param); 90 param = new G4UIparameter("y0",'d',true); 90 param = new G4UIparameter("y0",'d',true); 91 param->SetDefaultValue(0.); 91 param->SetDefaultValue(0.); 92 defCmd->SetParameter(param); 92 defCmd->SetParameter(param); 93 param = new G4UIparameter("z0",'d',true); 93 param = new G4UIparameter("z0",'d',true); 94 param->SetDefaultValue(0.); 94 param->SetDefaultValue(0.); 95 defCmd->SetParameter(param); 95 defCmd->SetParameter(param); 96 defCmd->AvailableForStates(G4State_Idle); 96 defCmd->AvailableForStates(G4State_Idle); 97 defCmd->SetToBeBroadcasted(false); 97 defCmd->SetToBeBroadcasted(false); 98 } 98 } 99 99 100 GRPrimGenActionMessenger::~GRPrimGenActionMess 100 GRPrimGenActionMessenger::~GRPrimGenActionMessenger() 101 { 101 { 102 delete defCmd; 102 delete defCmd; 103 delete srcDir; 103 delete srcDir; 104 } 104 } 105 105 106 #include "G4UImanager.hh" 106 #include "G4UImanager.hh" 107 #include "G4Threading.hh" 107 #include "G4Threading.hh" 108 #include "G4Tokenizer.hh" 108 #include "G4Tokenizer.hh" 109 #include "G4UnitsTable.hh" 109 #include "G4UnitsTable.hh" 110 #include <fstream> 110 #include <fstream> 111 #include "GRDetectorConstruction.hh" 111 #include "GRDetectorConstruction.hh" 112 112 113 void GRPrimGenActionMessenger::SetNewValue(G4U 113 void GRPrimGenActionMessenger::SetNewValue(G4UIcommand * command,G4String newVal) 114 { 114 { 115 auto UI = G4UImanager::GetUIpointer(); 115 auto UI = G4UImanager::GetUIpointer(); 116 if(verboseLevel>0) G4cout << newVal << G4end 116 if(verboseLevel>0) G4cout << newVal << G4endl; 117 117 118 G4Tokenizer next(newVal); 118 G4Tokenizer next(newVal); 119 G4String pName = next(); 119 G4String pName = next(); 120 G4String fName = next(); 120 G4String fName = next(); 121 std::ifstream specFile; 121 std::ifstream specFile; 122 specFile.open(fName,std::ios::in); 122 specFile.open(fName,std::ios::in); 123 if(specFile.fail()) 123 if(specFile.fail()) 124 { 124 { 125 G4ExceptionDescription ed; 125 G4ExceptionDescription ed; 126 ed << "ERROR : File <" << fName << "> is n 126 ed << "ERROR : File <" << fName << "> is not found."; 127 command->CommandFailed(ed); 127 command->CommandFailed(ed); 128 return; 128 return; 129 } 129 } 130 130 131 G4String srcType = next(); 131 G4String srcType = next(); 132 G4String radius = next(); 132 G4String radius = next(); 133 G4String unit = next(); 133 G4String unit = next(); 134 G4double r = StoD(radius); 134 G4double r = StoD(radius); 135 if(r<0.) 135 if(r<0.) 136 { 136 { 137 r = GRDetectorConstruction::GetWorldSize() 137 r = GRDetectorConstruction::GetWorldSize() * 0.98; 138 radius = DtoS(r); 138 radius = DtoS(r); 139 unit = "mm"; 139 unit = "mm"; 140 G4cout<<"################ Radius of the GP 140 G4cout<<"################ Radius of the GPS sphere is set to "<<r<<" (mm)"<<G4endl; 141 } 141 } 142 G4String pos = next() + " " + next() + " " + 142 G4String pos = next() + " " + next() + " " + next(); 143 143 144 G4String cmd; 144 G4String cmd; 145 G4int ec = 0; 145 G4int ec = 0; 146 ec = std::max(UI->ApplyCommand("/gps/pos/sha 146 ec = std::max(UI->ApplyCommand("/gps/pos/shape Sphere"),ec); 147 cmd = "/gps/pos/centre " + pos + " " + unit; 147 cmd = "/gps/pos/centre " + pos + " " + unit; 148 ec = std::max(UI->ApplyCommand(cmd),ec); 148 ec = std::max(UI->ApplyCommand(cmd),ec); 149 cmd = "/gps/pos/radius " + radius + " " + un 149 cmd = "/gps/pos/radius " + radius + " " + unit; 150 ec = std::max(UI->ApplyCommand(cmd),ec); 150 ec = std::max(UI->ApplyCommand(cmd),ec); 151 ec = std::max(UI->ApplyCommand("/gps/pos/typ 151 ec = std::max(UI->ApplyCommand("/gps/pos/type Surface"),ec); 152 ec = std::max(UI->ApplyCommand("/gps/number 152 ec = std::max(UI->ApplyCommand("/gps/number 1"),ec); 153 cmd = "/gps/particle " + pName; 153 cmd = "/gps/particle " + pName; 154 ec = std::max(UI->ApplyCommand(cmd),ec); 154 ec = std::max(UI->ApplyCommand(cmd),ec); 155 ec = std::max(UI->ApplyCommand("/gps/ang/typ 155 ec = std::max(UI->ApplyCommand("/gps/ang/type cos"),ec); 156 ec = std::max(UI->ApplyCommand("/gps/ang/max 156 ec = std::max(UI->ApplyCommand("/gps/ang/maxtheta 90.0 deg"),ec); 157 ec = std::max(UI->ApplyCommand("/gps/ang/min 157 ec = std::max(UI->ApplyCommand("/gps/ang/mintheta 0.0 deg"),ec); 158 cmd = "/gps/ene/type " + srcType; 158 cmd = "/gps/ene/type " + srcType; 159 ec = std::max(UI->ApplyCommand(cmd),ec); 159 ec = std::max(UI->ApplyCommand(cmd),ec); 160 ec = std::max(UI->ApplyCommand("/gps/hist/ty 160 ec = std::max(UI->ApplyCommand("/gps/hist/type arb"),ec); 161 161 162 enum { bufsize = 128 }; 162 enum { bufsize = 128 }; 163 static char* line = new char[bufsize]; 163 static char* line = new char[bufsize]; 164 while(specFile.good()) 164 while(specFile.good()) 165 { 165 { 166 specFile.getline(line,bufsize); 166 specFile.getline(line,bufsize); 167 if(line[(size_t)0]=='#') continue; 167 if(line[(size_t)0]=='#') continue; 168 G4String valStr = line; 168 G4String valStr = line; 169 G4StrUtil::strip(valStr); 169 G4StrUtil::strip(valStr); 170 G4StrUtil::rstrip(valStr, 0x0d); 170 G4StrUtil::rstrip(valStr, 0x0d); 171 if(valStr.size()==0) continue; 171 if(valStr.size()==0) continue; 172 cmd = "/gps/hist/point " + valStr; 172 cmd = "/gps/hist/point " + valStr; 173 ec = std::max(UI->ApplyCommand(cmd),ec); 173 ec = std::max(UI->ApplyCommand(cmd),ec); 174 if(specFile.eof()) break; 174 if(specFile.eof()) break; 175 } 175 } 176 ec = std::max(UI->ApplyCommand("/gps/hist/in 176 ec = std::max(UI->ApplyCommand("/gps/hist/inter Lin"),ec); 177 177 178 if(ec>0) 178 if(ec>0) 179 { 179 { 180 G4ExceptionDescription ed; 180 G4ExceptionDescription ed; 181 ed << "ERROR : Internal error while proces 181 ed << "ERROR : Internal error while processing /gorad/source/define command."; 182 command->CommandFailed(ec,ed); 182 command->CommandFailed(ec,ed); 183 } 183 } 184 } 184 } 185 185 186 G4String GRPrimGenActionMessenger::GetCurrentV 186 G4String GRPrimGenActionMessenger::GetCurrentValue(G4UIcommand* /*command*/) 187 { 187 { 188 G4String val; 188 G4String val; 189 return val; 189 return val; 190 } 190 } 191 191 192 #include "G4ParticleTable.hh" 192 #include "G4ParticleTable.hh" 193 void GRPrimGenActionMessenger::UpdateParticleL 193 void GRPrimGenActionMessenger::UpdateParticleList() 194 { 194 { 195 auto particleTable = G4ParticleTable::GetPar 195 auto particleTable = G4ParticleTable::GetParticleTable(); 196 G4String candList; 196 G4String candList; 197 for(G4int i=0;i<particleTable->entries();i++ 197 for(G4int i=0;i<particleTable->entries();i++) 198 { candList += particleTable->GetParticleName 198 { candList += particleTable->GetParticleName(i) + " "; } 199 candList += "ion"; 199 candList += "ion"; 200 particlePara->SetParameterCandidates(candLis 200 particlePara->SetParameterCandidates(candList); 201 } 201 } 202 202 203 203