| 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 // ------------------------------------------- 26 // -------------------------------------------------------------------
27 // 27 //
28 // GEANT4 Class file 28 // GEANT4 Class file
29 // 29 //
30 // File name: G4EmExtraParametersMessenger 30 // File name: G4EmExtraParametersMessenger
31 // 31 //
32 // Author: Vladimir Ivanchenko 32 // Author: Vladimir Ivanchenko
33 // 33 //
34 // Creation date: 07-05-2019 34 // Creation date: 07-05-2019
35 // 35 //
36 // ------------------------------------------- 36 // -------------------------------------------------------------------
37 // 37 //
38 38
39 //....oooOO0OOooo........oooOO0OOooo........oo 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
40 //....oooOO0OOooo........oooOO0OOooo........oo 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
41 41
42 #include "G4EmExtraParametersMessenger.hh" 42 #include "G4EmExtraParametersMessenger.hh"
43 #include "G4UIcommand.hh" 43 #include "G4UIcommand.hh"
44 #include "G4UIparameter.hh" 44 #include "G4UIparameter.hh"
45 #include "G4UIcmdWithABool.hh" 45 #include "G4UIcmdWithABool.hh"
46 #include "G4UIcmdWithAnInteger.hh" 46 #include "G4UIcmdWithAnInteger.hh"
47 #include "G4UIcmdWithADouble.hh" 47 #include "G4UIcmdWithADouble.hh"
48 #include "G4UIcmdWithADoubleAndUnit.hh" 48 #include "G4UIcmdWithADoubleAndUnit.hh"
49 #include "G4UIcmdWithAString.hh" 49 #include "G4UIcmdWithAString.hh"
50 #include "G4UIcmdWith3VectorAndUnit.hh" 50 #include "G4UIcmdWith3VectorAndUnit.hh"
51 #include "G4UImanager.hh" 51 #include "G4UImanager.hh"
52 #include "G4EmExtraParameters.hh" 52 #include "G4EmExtraParameters.hh"
53 53
54 #include <sstream> 54 #include <sstream>
55 55
56 //....oooOO0OOooo........oooOO0OOooo........oo 56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
57 57
58 G4EmExtraParametersMessenger::G4EmExtraParamet 58 G4EmExtraParametersMessenger::G4EmExtraParametersMessenger(G4EmExtraParameters* ptr)
59 : theParameters(ptr) 59 : theParameters(ptr)
60 { 60 {
61 paiCmd = new G4UIcommand("/process/em/AddPAI 61 paiCmd = new G4UIcommand("/process/em/AddPAIRegion",this);
62 paiCmd->SetGuidance("Activate PAI in the G4R 62 paiCmd->SetGuidance("Activate PAI in the G4Region.");
63 paiCmd->SetGuidance(" partName : particle 63 paiCmd->SetGuidance(" partName : particle name (default - all)");
64 paiCmd->SetGuidance(" regName : G4Region 64 paiCmd->SetGuidance(" regName : G4Region name");
65 paiCmd->SetGuidance(" paiType : PAI, PAIp 65 paiCmd->SetGuidance(" paiType : PAI, PAIphoton");
66 paiCmd->AvailableForStates(G4State_PreInit); 66 paiCmd->AvailableForStates(G4State_PreInit);
67 paiCmd->SetToBeBroadcasted(false); 67 paiCmd->SetToBeBroadcasted(false);
68 68
69 auto part = new G4UIparameter("partName",'s' 69 auto part = new G4UIparameter("partName",'s',false);
70 paiCmd->SetParameter(part); 70 paiCmd->SetParameter(part);
71 71
72 auto pregName = new G4UIparameter("regName", 72 auto pregName = new G4UIparameter("regName",'s',false);
73 paiCmd->SetParameter(pregName); 73 paiCmd->SetParameter(pregName);
74 74
75 auto ptype = new G4UIparameter("type",'s',fa 75 auto ptype = new G4UIparameter("type",'s',false);
76 paiCmd->SetParameter(ptype); 76 paiCmd->SetParameter(ptype);
77 ptype->SetParameterCandidates("pai PAI PAIph 77 ptype->SetParameterCandidates("pai PAI PAIphoton");
78 78
79 mscoCmd = new G4UIcommand("/process/em/AddEm 79 mscoCmd = new G4UIcommand("/process/em/AddEmRegion",this);
80 mscoCmd->SetGuidance("Add optional EM config 80 mscoCmd->SetGuidance("Add optional EM configuration for a G4Region.");
81 mscoCmd->SetGuidance(" regName : G4Region 81 mscoCmd->SetGuidance(" regName : G4Region name");
82 mscoCmd->SetGuidance(" emType : G4EmStand 82 mscoCmd->SetGuidance(" emType : G4EmStandard, G4EmStandard_opt1, ...");
83 mscoCmd->AvailableForStates(G4State_PreInit) 83 mscoCmd->AvailableForStates(G4State_PreInit);
84 mscoCmd->SetToBeBroadcasted(false); 84 mscoCmd->SetToBeBroadcasted(false);
85 85
86 auto mregName = new G4UIparameter("regName", 86 auto mregName = new G4UIparameter("regName",'s',false);
87 mscoCmd->SetParameter(mregName); 87 mscoCmd->SetParameter(mregName);
88 88
89 auto mtype = new G4UIparameter("mscType",'s' 89 auto mtype = new G4UIparameter("mscType",'s',false);
90 mscoCmd->SetParameter(mtype); 90 mscoCmd->SetParameter(mtype);
91 mtype->SetParameterCandidates("G4EmStandard 91 mtype->SetParameterCandidates("G4EmStandard G4EmStandard_opt1 G4EmStandard_opt2 G4EmStandard_opt3 G4EmStandard_opt4 G4EmStandardGS G4EmStandardSS G4EmLivermore G4EmPenelope G4RadioactiveDecay");
92 92
93 SubSecCmd = new G4UIcmdWithAString("/process 93 SubSecCmd = new G4UIcmdWithAString("/process/eLoss/subsecRegion",this);
94 SubSecCmd->SetGuidance("Enable subcut genera 94 SubSecCmd->SetGuidance("Enable subcut generation per region.");
95 SubSecCmd->SetGuidance(" Region : region 95 SubSecCmd->SetGuidance(" Region : region name");
96 SubSecCmd->AvailableForStates(G4State_PreIni 96 SubSecCmd->AvailableForStates(G4State_PreInit);
97 SubSecCmd->SetToBeBroadcasted(false); 97 SubSecCmd->SetToBeBroadcasted(false);
98 98
99 StepFuncCmd = new G4UIcommand("/process/eLos 99 StepFuncCmd = new G4UIcommand("/process/eLoss/StepFunction",this);
100 StepFuncCmd->SetGuidance("Set the energy los 100 StepFuncCmd->SetGuidance("Set the energy loss step limitation parameters for e+-.");
101 StepFuncCmd->SetGuidance(" dRoverR : max 101 StepFuncCmd->SetGuidance(" dRoverR : max Range variation per step");
102 StepFuncCmd->SetGuidance(" finalRange: rang 102 StepFuncCmd->SetGuidance(" finalRange: range for final step");
103 StepFuncCmd->SetGuidance(" unit : unit 103 StepFuncCmd->SetGuidance(" unit : unit of finalRange");
104 StepFuncCmd->AvailableForStates(G4State_PreI 104 StepFuncCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
105 StepFuncCmd->SetToBeBroadcasted(false); 105 StepFuncCmd->SetToBeBroadcasted(false);
106 106
107 auto dRoverRPrm = new G4UIparameter("dRoverR 107 auto dRoverRPrm = new G4UIparameter("dRoverR",'d',false);
108 dRoverRPrm->SetParameterRange("dRoverR>0. && 108 dRoverRPrm->SetParameterRange("dRoverR>0. && dRoverR<=1.");
109 StepFuncCmd->SetParameter(dRoverRPrm); 109 StepFuncCmd->SetParameter(dRoverRPrm);
110 110
111 auto finalRangePrm = new G4UIparameter("fina 111 auto finalRangePrm = new G4UIparameter("finalRange",'d',false);
112 finalRangePrm->SetParameterRange("finalRange 112 finalRangePrm->SetParameterRange("finalRange>0.");
113 StepFuncCmd->SetParameter(finalRangePrm); 113 StepFuncCmd->SetParameter(finalRangePrm);
114 114
115 auto unitPrm = new G4UIparameter("unit",'s', 115 auto unitPrm = new G4UIparameter("unit",'s',true);
116 unitPrm->SetDefaultUnit("mm"); 116 unitPrm->SetDefaultUnit("mm");
117 StepFuncCmd->SetParameter(unitPrm); 117 StepFuncCmd->SetParameter(unitPrm);
118 118
119 StepFuncCmd1 = new G4UIcommand("/process/eLo 119 StepFuncCmd1 = new G4UIcommand("/process/eLoss/StepFunctionMuHad",this);
120 StepFuncCmd1->SetGuidance("Set the energy lo 120 StepFuncCmd1->SetGuidance("Set the energy loss step limitation parameters for muon/hadron.");
121 StepFuncCmd1->SetGuidance(" dRoverR : max 121 StepFuncCmd1->SetGuidance(" dRoverR : max Range variation per step");
122 StepFuncCmd1->SetGuidance(" finalRange: ran 122 StepFuncCmd1->SetGuidance(" finalRange: range for final step");
123 StepFuncCmd1->AvailableForStates(G4State_Pre 123 StepFuncCmd1->AvailableForStates(G4State_PreInit,G4State_Idle);
124 StepFuncCmd1->SetToBeBroadcasted(false); 124 StepFuncCmd1->SetToBeBroadcasted(false);
125 125
126 auto dRoverRPrm1 = new G4UIparameter("dRover 126 auto dRoverRPrm1 = new G4UIparameter("dRoverRMuHad",'d',false);
127 dRoverRPrm1->SetParameterRange("dRoverRMuHad 127 dRoverRPrm1->SetParameterRange("dRoverRMuHad>0. && dRoverRMuHad<=1.");
128 StepFuncCmd1->SetParameter(dRoverRPrm1); 128 StepFuncCmd1->SetParameter(dRoverRPrm1);
129 129
130 auto finalRangePrm1 = new G4UIparameter("fin 130 auto finalRangePrm1 = new G4UIparameter("finalRangeMuHad",'d',false);
131 finalRangePrm1->SetParameterRange("finalRang 131 finalRangePrm1->SetParameterRange("finalRangeMuHad>0.");
132 StepFuncCmd1->SetParameter(finalRangePrm1); 132 StepFuncCmd1->SetParameter(finalRangePrm1);
133 133
134 auto unitPrm1 = new G4UIparameter("unit",'s' 134 auto unitPrm1 = new G4UIparameter("unit",'s',true);
135 unitPrm1->SetDefaultValue("mm"); 135 unitPrm1->SetDefaultValue("mm");
136 StepFuncCmd1->SetParameter(unitPrm1); 136 StepFuncCmd1->SetParameter(unitPrm1);
137 137
138 StepFuncCmd2 = new G4UIcommand("/process/eLo 138 StepFuncCmd2 = new G4UIcommand("/process/eLoss/StepFunctionLightIons",this);
139 StepFuncCmd2->SetGuidance("Set the energy lo 139 StepFuncCmd2->SetGuidance("Set the energy loss step limitation parameters for light ions.");
140 StepFuncCmd2->SetGuidance(" dRoverR : max 140 StepFuncCmd2->SetGuidance(" dRoverR : max Range variation per step");
141 StepFuncCmd2->SetGuidance(" finalRange: ran 141 StepFuncCmd2->SetGuidance(" finalRange: range for final step");
142 StepFuncCmd2->AvailableForStates(G4State_Pre 142 StepFuncCmd2->AvailableForStates(G4State_PreInit,G4State_Idle);
143 StepFuncCmd2->SetToBeBroadcasted(false); 143 StepFuncCmd2->SetToBeBroadcasted(false);
144 144
145 auto dRoverRPrm2 = new G4UIparameter("dRover 145 auto dRoverRPrm2 = new G4UIparameter("dRoverRLIons",'d',false);
146 dRoverRPrm2->SetParameterRange("dRoverRLIons 146 dRoverRPrm2->SetParameterRange("dRoverRLIons>0. && dRoverRLIons<=1.");
147 StepFuncCmd2->SetParameter(dRoverRPrm2); 147 StepFuncCmd2->SetParameter(dRoverRPrm2);
148 148
149 auto finalRangePrm2 = new G4UIparameter("fin 149 auto finalRangePrm2 = new G4UIparameter("finalRangeLIons",'d',false);
150 finalRangePrm2->SetParameterRange("finalRang 150 finalRangePrm2->SetParameterRange("finalRangeLIons>0.");
151 StepFuncCmd2->SetParameter(finalRangePrm2); 151 StepFuncCmd2->SetParameter(finalRangePrm2);
152 152
153 auto unitPrm2 = new G4UIparameter("unit",'s' 153 auto unitPrm2 = new G4UIparameter("unit",'s',true);
154 unitPrm2->SetDefaultValue("mm"); 154 unitPrm2->SetDefaultValue("mm");
155 StepFuncCmd2->SetParameter(unitPrm2); 155 StepFuncCmd2->SetParameter(unitPrm2);
156 156
157 StepFuncCmd3 = new G4UIcommand("/process/eLo 157 StepFuncCmd3 = new G4UIcommand("/process/eLoss/StepFunctionIons",this);
158 StepFuncCmd3->SetGuidance("Set the energy lo 158 StepFuncCmd3->SetGuidance("Set the energy loss step limitation parameters for ions.");
159 StepFuncCmd3->SetGuidance(" dRoverR : max 159 StepFuncCmd3->SetGuidance(" dRoverR : max Range variation per step");
160 StepFuncCmd3->SetGuidance(" finalRange: ran 160 StepFuncCmd3->SetGuidance(" finalRange: range for final step");
161 StepFuncCmd3->AvailableForStates(G4State_Pre 161 StepFuncCmd3->AvailableForStates(G4State_PreInit,G4State_Idle);
162 StepFuncCmd3->SetToBeBroadcasted(false); 162 StepFuncCmd3->SetToBeBroadcasted(false);
163 163
164 auto dRoverRPrm3 = new G4UIparameter("dRover 164 auto dRoverRPrm3 = new G4UIparameter("dRoverRIons",'d',false);
165 dRoverRPrm3->SetParameterRange("dRoverRIons> 165 dRoverRPrm3->SetParameterRange("dRoverRIons>0. && dRoverRIons<=1.");
166 StepFuncCmd3->SetParameter(dRoverRPrm3); 166 StepFuncCmd3->SetParameter(dRoverRPrm3);
167 167
168 auto finalRangePrm3 = new G4UIparameter("fin 168 auto finalRangePrm3 = new G4UIparameter("finalRangeIons",'d',false);
169 finalRangePrm3->SetParameterRange("finalRang 169 finalRangePrm3->SetParameterRange("finalRangeIons>0.");
170 StepFuncCmd3->SetParameter(finalRangePrm3); 170 StepFuncCmd3->SetParameter(finalRangePrm3);
171 171
172 auto unitPrm3 = new G4UIparameter("unit",'s' 172 auto unitPrm3 = new G4UIparameter("unit",'s',true);
173 unitPrm3->SetDefaultValue("mm"); 173 unitPrm3->SetDefaultValue("mm");
174 StepFuncCmd3->SetParameter(unitPrm3); 174 StepFuncCmd3->SetParameter(unitPrm3);
175 175
176 bfCmd = new G4UIcommand("/process/em/setBias 176 bfCmd = new G4UIcommand("/process/em/setBiasingFactor",this);
177 bfCmd->SetGuidance("Set factor for the proce 177 bfCmd->SetGuidance("Set factor for the process cross section.");
178 bfCmd->SetGuidance(" procName : process n 178 bfCmd->SetGuidance(" procName : process name");
179 bfCmd->SetGuidance(" procFact : factor"); 179 bfCmd->SetGuidance(" procFact : factor");
180 bfCmd->SetGuidance(" flagFact : flag to c 180 bfCmd->SetGuidance(" flagFact : flag to change weight");
181 bfCmd->AvailableForStates(G4State_PreInit,G4 181 bfCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
182 bfCmd->SetToBeBroadcasted(false); 182 bfCmd->SetToBeBroadcasted(false);
183 183
184 auto procName = new G4UIparameter("procName" 184 auto procName = new G4UIparameter("procName",'s',false);
185 bfCmd->SetParameter(procName); 185 bfCmd->SetParameter(procName);
186 186
187 auto procFact = new G4UIparameter("procFact" 187 auto procFact = new G4UIparameter("procFact",'d',false);
188 bfCmd->SetParameter(procFact); 188 bfCmd->SetParameter(procFact);
189 189
190 auto flagFact = new G4UIparameter("flagFact" 190 auto flagFact = new G4UIparameter("flagFact",'s',false);
191 bfCmd->SetParameter(flagFact); 191 bfCmd->SetParameter(flagFact);
192 192
193 fiCmd = new G4UIcommand("/process/em/setForc 193 fiCmd = new G4UIcommand("/process/em/setForcedInteraction",this);
194 fiCmd->SetGuidance("Set factor for the proce 194 fiCmd->SetGuidance("Set factor for the process cross section.");
195 fiCmd->SetGuidance(" procNam : process n 195 fiCmd->SetGuidance(" procNam : process name");
196 fiCmd->SetGuidance(" regNam : region na 196 fiCmd->SetGuidance(" regNam : region name");
197 fiCmd->SetGuidance(" tlength : fixed tar 197 fiCmd->SetGuidance(" tlength : fixed target length");
198 fiCmd->SetGuidance(" unitT : length un 198 fiCmd->SetGuidance(" unitT : length unit");
199 fiCmd->SetGuidance(" tflag : flag to c 199 fiCmd->SetGuidance(" tflag : flag to change weight");
200 fiCmd->AvailableForStates(G4State_PreInit,G4 200 fiCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
201 fiCmd->SetToBeBroadcasted(false); 201 fiCmd->SetToBeBroadcasted(false);
202 202
203 auto procNam = new G4UIparameter("procNam",' 203 auto procNam = new G4UIparameter("procNam",'s',false);
204 fiCmd->SetParameter(procNam); 204 fiCmd->SetParameter(procNam);
205 205
206 auto regNam = new G4UIparameter("regNam",'s 206 auto regNam = new G4UIparameter("regNam",'s',false);
207 fiCmd->SetParameter(regNam); 207 fiCmd->SetParameter(regNam);
208 208
209 auto tlength = new G4UIparameter("tlength",' 209 auto tlength = new G4UIparameter("tlength",'d',false);
210 tlength->SetParameterRange("tlength>0"); 210 tlength->SetParameterRange("tlength>0");
211 fiCmd->SetParameter(tlength); 211 fiCmd->SetParameter(tlength);
212 212
213 auto unitT = new G4UIparameter("unitT",'s',t 213 auto unitT = new G4UIparameter("unitT",'s',true);
214 unitT->SetDefaultUnit("mm"); 214 unitT->SetDefaultUnit("mm");
215 fiCmd->SetParameter(unitT); 215 fiCmd->SetParameter(unitT);
216 216
217 auto flagT = new G4UIparameter("tflag",'b',t 217 auto flagT = new G4UIparameter("tflag",'b',true);
218 flagT->SetDefaultValue(true); 218 flagT->SetDefaultValue(true);
219 fiCmd->SetParameter(flagT); 219 fiCmd->SetParameter(flagT);
220 220
221 bsCmd = new G4UIcommand("/process/em/setSecB 221 bsCmd = new G4UIcommand("/process/em/setSecBiasing",this);
222 bsCmd->SetGuidance("Set bremsstrahlung or de 222 bsCmd->SetGuidance("Set bremsstrahlung or delta-e- splitting/Russian roulette per region.");
223 bsCmd->SetGuidance(" bProcNam : process nam 223 bsCmd->SetGuidance(" bProcNam : process name");
224 bsCmd->SetGuidance(" bRegNam : region name 224 bsCmd->SetGuidance(" bRegNam : region name");
225 bsCmd->SetGuidance(" bFactor : number of s 225 bsCmd->SetGuidance(" bFactor : number of split gamma or probability of Russian roulette");
226 bsCmd->SetGuidance(" bEnergy : max energy 226 bsCmd->SetGuidance(" bEnergy : max energy of a secondary for this biasing method");
227 bsCmd->SetGuidance(" bUnit : energy unit 227 bsCmd->SetGuidance(" bUnit : energy unit");
228 bsCmd->AvailableForStates(G4State_PreInit,G4 228 bsCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
229 bsCmd->SetToBeBroadcasted(false); 229 bsCmd->SetToBeBroadcasted(false);
230 230
231 auto bProcNam = new G4UIparameter("bProcNam" 231 auto bProcNam = new G4UIparameter("bProcNam",'s',false);
232 bsCmd->SetParameter(bProcNam); 232 bsCmd->SetParameter(bProcNam);
233 233
234 auto bRegNam = new G4UIparameter("bRegNam",' 234 auto bRegNam = new G4UIparameter("bRegNam",'s',false);
235 bsCmd->SetParameter(bRegNam); 235 bsCmd->SetParameter(bRegNam);
236 236
237 auto bFactor = new G4UIparameter("bFactor",' 237 auto bFactor = new G4UIparameter("bFactor",'d',false);
238 bsCmd->SetParameter(bFactor); 238 bsCmd->SetParameter(bFactor);
239 239
240 auto bEnergy = new G4UIparameter("bEnergy",' 240 auto bEnergy = new G4UIparameter("bEnergy",'d',false);
241 bsCmd->SetParameter(bEnergy); 241 bsCmd->SetParameter(bEnergy);
242 242
243 auto bUnit = new G4UIparameter("bUnit",'s',t 243 auto bUnit = new G4UIparameter("bUnit",'s',true);
244 bUnit->SetDefaultUnit("MeV"); 244 bUnit->SetDefaultUnit("MeV");
245 bsCmd->SetParameter(bUnit); 245 bsCmd->SetParameter(bUnit);
246 246
247 dirSplitCmd = new G4UIcmdWithABool("/process 247 dirSplitCmd = new G4UIcmdWithABool("/process/em/setDirectionalSplitting",this);
248 dirSplitCmd->SetGuidance("Enable directional 248 dirSplitCmd->SetGuidance("Enable directional brem splitting");
249 dirSplitCmd->AvailableForStates(G4State_PreI 249 dirSplitCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
250 dirSplitCmd->SetToBeBroadcasted(false); 250 dirSplitCmd->SetToBeBroadcasted(false);
251 251
252 qeCmd = new G4UIcmdWithABool("/process/em/Qu 252 qeCmd = new G4UIcmdWithABool("/process/em/QuantumEntanglement",this);
253 qeCmd->SetGuidance("Enable quantum entanglem 253 qeCmd->SetGuidance("Enable quantum entanglement");
254 qeCmd->AvailableForStates(G4State_PreInit,G4 254 qeCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
255 qeCmd->SetToBeBroadcasted(false); 255 qeCmd->SetToBeBroadcasted(false);
256 256
257 dirSplitTargetCmd = new G4UIcmdWith3VectorAn 257 dirSplitTargetCmd = new G4UIcmdWith3VectorAndUnit("/process/em/setDirectionalSplittingTarget",this);
258 dirSplitTargetCmd->SetGuidance("Position of 258 dirSplitTargetCmd->SetGuidance("Position of arget for directional splitting");
259 dirSplitTargetCmd->AvailableForStates(G4Stat 259 dirSplitTargetCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
260 260
261 dirSplitRadiusCmd = new G4UIcmdWithADoubleAn 261 dirSplitRadiusCmd = new G4UIcmdWithADoubleAndUnit("/process/em/setDirectionalSplittingRadius",this);
262 dirSplitRadiusCmd->SetGuidance("Radius of ta 262 dirSplitRadiusCmd->SetGuidance("Radius of target for directional splitting");
263 dirSplitRadiusCmd->AvailableForStates(G4Stat 263 dirSplitRadiusCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
264 dirSplitRadiusCmd->SetToBeBroadcasted(false) 264 dirSplitRadiusCmd->SetToBeBroadcasted(false);
265 } 265 }
266 266
267 //....oooOO0OOooo........oooOO0OOooo........oo 267 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
268 268
269 G4EmExtraParametersMessenger::~G4EmExtraParame 269 G4EmExtraParametersMessenger::~G4EmExtraParametersMessenger()
270 { 270 {
271 delete paiCmd; 271 delete paiCmd;
272 delete mscoCmd; 272 delete mscoCmd;
273 delete SubSecCmd; 273 delete SubSecCmd;
274 delete bfCmd; 274 delete bfCmd;
275 delete fiCmd; 275 delete fiCmd;
276 delete bsCmd; 276 delete bsCmd;
277 delete qeCmd; 277 delete qeCmd;
278 delete StepFuncCmd; 278 delete StepFuncCmd;
279 delete StepFuncCmd1; 279 delete StepFuncCmd1;
280 delete StepFuncCmd2; 280 delete StepFuncCmd2;
281 delete StepFuncCmd3; 281 delete StepFuncCmd3;
282 delete dirSplitCmd; 282 delete dirSplitCmd;
283 delete dirSplitTargetCmd; 283 delete dirSplitTargetCmd;
284 delete dirSplitRadiusCmd; 284 delete dirSplitRadiusCmd;
285 } 285 }
286 286
287 //....oooOO0OOooo........oooOO0OOooo........oo 287 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
288 288
289 void G4EmExtraParametersMessenger::SetNewValue 289 void G4EmExtraParametersMessenger::SetNewValue(G4UIcommand* command,
290 290 G4String newValue)
291 { 291 {
292 G4bool physicsModified = false; 292 G4bool physicsModified = false;
293 293
294 if (command == paiCmd) { 294 if (command == paiCmd) {
295 G4String s1(""),s2(""),s3(""); 295 G4String s1(""),s2(""),s3("");
296 std::istringstream is(newValue); 296 std::istringstream is(newValue);
297 is >> s1 >> s2 >> s3; 297 is >> s1 >> s2 >> s3;
298 theParameters->AddPAIModel(s1, s2, s3); 298 theParameters->AddPAIModel(s1, s2, s3);
299 } else if (command == mscoCmd) { 299 } else if (command == mscoCmd) {
300 G4String s1(""),s2(""); 300 G4String s1(""),s2("");
301 std::istringstream is(newValue); 301 std::istringstream is(newValue);
302 is >> s1 >> s2; 302 is >> s1 >> s2;
303 theParameters->AddPhysics(s1, s2); 303 theParameters->AddPhysics(s1, s2);
304 } else if (command == StepFuncCmd || command 304 } else if (command == StepFuncCmd || command == StepFuncCmd1 || command == StepFuncCmd2 || command == StepFuncCmd3) {
305 G4double v1,v2; 305 G4double v1,v2;
306 G4String unt; 306 G4String unt;
307 std::istringstream is(newValue); 307 std::istringstream is(newValue);
308 is >> v1 >> v2 >> unt; 308 is >> v1 >> v2 >> unt;
309 v2 *= G4UIcommand::ValueOf(unt); 309 v2 *= G4UIcommand::ValueOf(unt);
310 if(command == StepFuncCmd) { 310 if(command == StepFuncCmd) {
311 theParameters->SetStepFunction(v1,v2); 311 theParameters->SetStepFunction(v1,v2);
312 } else if(command == StepFuncCmd1) { 312 } else if(command == StepFuncCmd1) {
313 theParameters->SetStepFunctionMuHad(v1,v 313 theParameters->SetStepFunctionMuHad(v1,v2);
314 } else if(command == StepFuncCmd2) { 314 } else if(command == StepFuncCmd2) {
315 theParameters->SetStepFunctionLightIons( 315 theParameters->SetStepFunctionLightIons(v1,v2);
316 } else { 316 } else {
317 theParameters->SetStepFunctionIons(v1,v2 317 theParameters->SetStepFunctionIons(v1,v2);
318 } 318 }
319 physicsModified = true; 319 physicsModified = true;
320 } else if (command == SubSecCmd) { 320 } else if (command == SubSecCmd) {
321 theParameters->SetSubCutRegion(newValue); 321 theParameters->SetSubCutRegion(newValue);
322 } else if (command == bfCmd) { 322 } else if (command == bfCmd) {
323 G4double v1(1.0); 323 G4double v1(1.0);
324 G4String s0(""),s1(""); 324 G4String s0(""),s1("");
325 std::istringstream is(newValue); 325 std::istringstream is(newValue);
326 is >> s0 >> v1 >> s1; 326 is >> s0 >> v1 >> s1;
327 G4bool yes = false; 327 G4bool yes = false;
328 if(s1 == "true") { yes = true; } 328 if(s1 == "true") { yes = true; }
329 theParameters->SetProcessBiasingFactor(s0, 329 theParameters->SetProcessBiasingFactor(s0,v1,yes);
330 physicsModified = true; 330 physicsModified = true;
331 } else if (command == fiCmd) { 331 } else if (command == fiCmd) {
332 G4double v1(0.0); 332 G4double v1(0.0);
333 G4String s1(""),s2(""),s3(""),unt("mm"); 333 G4String s1(""),s2(""),s3(""),unt("mm");
334 std::istringstream is(newValue); 334 std::istringstream is(newValue);
335 is >> s1 >> s2 >> v1 >> unt >> s3; 335 is >> s1 >> s2 >> v1 >> unt >> s3;
336 G4bool yes = false; 336 G4bool yes = false;
337 if(s3 == "true") { yes = true; } 337 if(s3 == "true") { yes = true; }
338 v1 *= G4UIcommand::ValueOf(unt); 338 v1 *= G4UIcommand::ValueOf(unt);
339 theParameters->ActivateForcedInteraction(s 339 theParameters->ActivateForcedInteraction(s1,s2,v1,yes);
340 physicsModified = true; 340 physicsModified = true;
341 } else if (command == bsCmd) { 341 } else if (command == bsCmd) {
342 G4double fb(1.0),en(1.e+30); 342 G4double fb(1.0),en(1.e+30);
343 G4String s1(""),s2(""),unt("MeV"); 343 G4String s1(""),s2(""),unt("MeV");
344 std::istringstream is(newValue); 344 std::istringstream is(newValue);
345 is >> s1 >> s2 >> fb >> en >> unt; 345 is >> s1 >> s2 >> fb >> en >> unt;
346 en *= G4UIcommand::ValueOf(unt); 346 en *= G4UIcommand::ValueOf(unt);
347 theParameters->ActivateSecondaryBiasing(s1 347 theParameters->ActivateSecondaryBiasing(s1,s2,fb,en);
348 physicsModified = true; 348 physicsModified = true;
349 } else if (command == qeCmd) { 349 } else if (command == qeCmd) {
350 theParameters->SetQuantumEntanglement(qeCm 350 theParameters->SetQuantumEntanglement(qeCmd->GetNewBoolValue(newValue));
351 } else if (command == dirSplitCmd) { 351 } else if (command == dirSplitCmd) {
352 theParameters->SetDirectionalSplitting( 352 theParameters->SetDirectionalSplitting(
353 dirSplitCmd->GetNewBoolValue(newValue)); 353 dirSplitCmd->GetNewBoolValue(newValue));
354 physicsModified = true; 354 physicsModified = true;
355 } else if (command == dirSplitTargetCmd) { 355 } else if (command == dirSplitTargetCmd) {
356 G4ThreeVector t = dirSplitTargetCmd->GetNe 356 G4ThreeVector t = dirSplitTargetCmd->GetNew3VectorValue(newValue);
357 theParameters->SetDirectionalSplittingTarg 357 theParameters->SetDirectionalSplittingTarget(t);
358 physicsModified = true; 358 physicsModified = true;
359 } else if (command == dirSplitRadiusCmd) { 359 } else if (command == dirSplitRadiusCmd) {
360 G4double r = dirSplitRadiusCmd->GetNewDoub 360 G4double r = dirSplitRadiusCmd->GetNewDoubleValue(newValue);
361 theParameters->SetDirectionalSplittingRadi 361 theParameters->SetDirectionalSplittingRadius(r);
362 physicsModified = true; 362 physicsModified = true;
363 } 363 }
364 364
365 if(physicsModified) { 365 if(physicsModified) {
366 G4UImanager::GetUIpointer()->ApplyCommand( 366 G4UImanager::GetUIpointer()->ApplyCommand("/run/physicsModified");
367 } 367 }
368 } 368 }
369 369
370 //....oooOO0OOooo........oooOO0OOooo........oo 370 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
371 371