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