| 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 // File name: G4NistMessenger 27 // File name: G4NistMessenger
27 // 28 //
28 // Author: Vladimir Ivanchenko 29 // Author: Vladimir Ivanchenko
29 // 30 //
30 // Creation date: 23.12.2004 31 // Creation date: 23.12.2004
31 // 32 //
32 // Modifications: 33 // Modifications:
33 // 34 //
34 // 35 //
35 // ------------------------------------------- 36 // -------------------------------------------------------------------
>> 37 //
>> 38
>> 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
36 41
37 #include "G4NistMessenger.hh" 42 #include "G4NistMessenger.hh"
38 43
39 #include "G4ApplicationState.hh" <<
40 #include "G4DensityEffectData.hh" <<
41 #include "G4IonisParamMat.hh" <<
42 #include "G4NistManager.hh" 44 #include "G4NistManager.hh"
43 #include "G4UIcmdWithAString.hh" <<
44 #include "G4UIcmdWithAnInteger.hh" <<
45 #include "G4UIdirectory.hh" <<
46 #include "G4UIcommand.hh" <<
47 #include "G4UIparameter.hh" <<
48 45
49 #include <sstream> << 46 #include "G4UIdirectory.hh"
>> 47 #include "G4UIcmdWithAnInteger.hh"
>> 48 #include "G4UIcmdWithAString.hh"
>> 49 #include "G4IonisParamMat.hh"
>> 50 #include "G4DensityEffectData.hh"
50 51
51 //....oooOO0OOooo........oooOO0OOooo........oo 52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
52 53
53 G4NistMessenger::G4NistMessenger(G4NistManager << 54 G4NistMessenger::G4NistMessenger(G4NistManager* man)
>> 55 : manager(man)
54 { 56 {
55 matDir = new G4UIdirectory("/material/"); 57 matDir = new G4UIdirectory("/material/");
56 matDir->SetGuidance("Commands for materials" 58 matDir->SetGuidance("Commands for materials");
57 59
58 verCmd = new G4UIcmdWithAnInteger("/material << 60 verCmd = new G4UIcmdWithAnInteger("/material/verbose",this);
59 verCmd->SetGuidance("Set verbose level."); 61 verCmd->SetGuidance("Set verbose level.");
60 verCmd->AvailableForStates(G4State_PreInit, << 62
61 verCmd->SetToBeBroadcasted(false); <<
62 <<
63 nistDir = new G4UIdirectory("/material/nist/ 63 nistDir = new G4UIdirectory("/material/nist/");
64 nistDir->SetGuidance("Commands for the nist 64 nistDir->SetGuidance("Commands for the nist dataBase");
65 << 65
66 prtElmCmd = new G4UIcmdWithAString("/materia << 66 prtElmCmd = new G4UIcmdWithAString("/material/nist/printElement",this);
67 prtElmCmd->SetGuidance("print element(s) in 67 prtElmCmd->SetGuidance("print element(s) in dataBase.");
68 prtElmCmd->SetGuidance("symbol = element."); 68 prtElmCmd->SetGuidance("symbol = element.");
69 prtElmCmd->SetGuidance("all = all element 69 prtElmCmd->SetGuidance("all = all elements.");
70 prtElmCmd->SetParameterName("symbol", true); 70 prtElmCmd->SetParameterName("symbol", true);
71 prtElmCmd->SetDefaultValue("all"); 71 prtElmCmd->SetDefaultValue("all");
72 prtElmCmd->AvailableForStates(G4State_PreIni << 72
73 prtElmCmd->SetToBeBroadcasted(false); << 73 przElmCmd = new G4UIcmdWithAnInteger("/material/nist/printElementZ",this);
74 <<
75 przElmCmd = new G4UIcmdWithAnInteger("/mater <<
76 przElmCmd->SetGuidance("print element Z in d 74 przElmCmd->SetGuidance("print element Z in dataBase.");
77 przElmCmd->SetGuidance("0 = all elements."); 75 przElmCmd->SetGuidance("0 = all elements.");
78 przElmCmd->SetParameterName("Z", true); 76 przElmCmd->SetParameterName("Z", true);
79 przElmCmd->SetDefaultValue(0); 77 przElmCmd->SetDefaultValue(0);
80 przElmCmd->SetRange("0<=Z && Z<108"); 78 przElmCmd->SetRange("0<=Z && Z<108");
81 przElmCmd->AvailableForStates(G4State_PreIni << 79
82 przElmCmd->SetToBeBroadcasted(false); << 80 lisMatCmd = new G4UIcmdWithAString("/material/nist/listMaterials",this);
83 <<
84 lisMatCmd = new G4UIcmdWithAString("/materia <<
85 lisMatCmd->SetGuidance("Materials in Geant4 81 lisMatCmd->SetGuidance("Materials in Geant4 dataBase.");
86 lisMatCmd->SetGuidance("simple - simple NIST 82 lisMatCmd->SetGuidance("simple - simple NIST materials.");
87 lisMatCmd->SetGuidance("compound - compound 83 lisMatCmd->SetGuidance("compound - compound NIST materials.");
88 lisMatCmd->SetGuidance("hep - HEP materials. 84 lisMatCmd->SetGuidance("hep - HEP materials.");
89 lisMatCmd->SetGuidance("bio - biomedical mat 85 lisMatCmd->SetGuidance("bio - biomedical materials.");
90 lisMatCmd->SetGuidance("all - list of all Ge 86 lisMatCmd->SetGuidance("all - list of all Geant4 materials.");
91 lisMatCmd->SetParameterName("matlist", true) 87 lisMatCmd->SetParameterName("matlist", true);
92 // lisMatCmd->SetCandidates("simple compound 88 // lisMatCmd->SetCandidates("simple compound hep bio all");
93 lisMatCmd->SetDefaultValue("all"); 89 lisMatCmd->SetDefaultValue("all");
94 lisMatCmd->AvailableForStates(G4State_PreIni << 90
95 lisMatCmd->SetToBeBroadcasted(false); <<
96 <<
97 g4Dir = new G4UIdirectory("/material/g4/"); 91 g4Dir = new G4UIdirectory("/material/g4/");
98 g4Dir->SetGuidance("Commands for G4MaterialT 92 g4Dir->SetGuidance("Commands for G4MaterialTable");
99 << 93
100 g4ElmCmd = new G4UIcmdWithAString("/material << 94 g4ElmCmd = new G4UIcmdWithAString("/material/g4/printElement",this);
101 g4ElmCmd->SetGuidance("print Element from G4 95 g4ElmCmd->SetGuidance("print Element from G4ElementTable.");
102 g4ElmCmd->SetGuidance("all - all elements.") 96 g4ElmCmd->SetGuidance("all - all elements.");
103 g4ElmCmd->SetParameterName("elm", true); 97 g4ElmCmd->SetParameterName("elm", true);
104 g4ElmCmd->SetDefaultValue("all"); 98 g4ElmCmd->SetDefaultValue("all");
105 << 99
106 g4MatCmd = new G4UIcmdWithAString("/material << 100 g4MatCmd = new G4UIcmdWithAString("/material/g4/printMaterial",this);
107 g4MatCmd->SetGuidance("print Material from G 101 g4MatCmd->SetGuidance("print Material from G4MaterialTable.");
108 g4MatCmd->SetGuidance("all - all materials") 102 g4MatCmd->SetGuidance("all - all materials");
109 g4MatCmd->SetParameterName("pmat", true); 103 g4MatCmd->SetParameterName("pmat", true);
110 g4MatCmd->SetDefaultValue("all"); 104 g4MatCmd->SetDefaultValue("all");
111 g4MatCmd->AvailableForStates(G4State_PreInit <<
112 g4MatCmd->SetToBeBroadcasted(false); <<
113 105
114 g4DensCmd = new G4UIcmdWithAString("/materia << 106 g4DensCmd = new G4UIcmdWithAString("/material/g4/printDensityEffParam",this);
115 g4DensCmd->SetGuidance("print Material from 107 g4DensCmd->SetGuidance("print Material from G4DensityEffectData.");
116 g4DensCmd->SetGuidance("all - all materials" 108 g4DensCmd->SetGuidance("all - all materials");
117 g4DensCmd->SetParameterName("dmat", true); 109 g4DensCmd->SetParameterName("dmat", true);
118 g4DensCmd->SetDefaultValue("all"); 110 g4DensCmd->SetDefaultValue("all");
119 g4DensCmd->AvailableForStates(G4State_PreIni <<
120 g4DensCmd->SetToBeBroadcasted(false); <<
121 111
122 densCmd = new G4UIcmdWithAString("/material/ << 112 densCmd = new G4UIcmdWithAString("/material/g4/enableDensityEffOnFly",this);
123 densCmd->SetGuidance("enable accurate comput 113 densCmd->SetGuidance("enable accurate computation of density effect.");
124 densCmd->SetGuidance("all - all materials.") 114 densCmd->SetGuidance("all - all materials.");
125 densCmd->SetParameterName("dens", true); 115 densCmd->SetParameterName("dens", true);
126 densCmd->SetDefaultValue("all"); 116 densCmd->SetDefaultValue("all");
127 densCmd->AvailableForStates(G4State_PreInit) <<
128 lisMatCmd->SetToBeBroadcasted(false); <<
129 117
130 adensCmd = new G4UIcmdWithAString("/material << 118 adensCmd = new G4UIcmdWithAString("/material/g4/disableDensityEffOnFly",this);
131 adensCmd->SetGuidance("disable accurate comp 119 adensCmd->SetGuidance("disable accurate computation of density effect.");
132 adensCmd->SetGuidance("all - all materials." 120 adensCmd->SetGuidance("all - all materials.");
133 adensCmd->SetParameterName("dens", true); 121 adensCmd->SetParameterName("dens", true);
134 adensCmd->SetDefaultValue("all"); 122 adensCmd->SetDefaultValue("all");
135 adensCmd->AvailableForStates(G4State_PreInit <<
136 adensCmd->SetToBeBroadcasted(false); <<
137 <<
138 fPosiCmd = new G4UIcommand("/material/g4/ort <<
139 fPosiCmd->SetGuidance("defined orto-positron <<
140 fPosiCmd->SetGuidance("via material name, al <<
141 fPosiCmd->AvailableForStates(G4State_PreInit <<
142 fPosiCmd->SetToBeBroadcasted(false); <<
143 <<
144 auto p1 = new G4UIparameter("matname",'s', f <<
145 fPosiCmd->SetParameter(p1); <<
146 <<
147 auto p2 = new G4UIparameter("fraction", 'd', <<
148 p2->SetParameterRange("fraction>=0. && fract <<
149 fPosiCmd->SetParameter(p2); <<
150 } 123 }
151 124
152 //....oooOO0OOooo........oooOO0OOooo........oo 125 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
153 126
154 G4NistMessenger::~G4NistMessenger() 127 G4NistMessenger::~G4NistMessenger()
155 { 128 {
156 delete verCmd; << 129 delete verCmd;
157 delete prtElmCmd; 130 delete prtElmCmd;
158 delete przElmCmd; 131 delete przElmCmd;
159 delete lisMatCmd; 132 delete lisMatCmd;
160 delete nistDir; 133 delete nistDir;
161 << 134
162 delete g4ElmCmd; << 135 delete g4ElmCmd;
163 delete g4MatCmd; 136 delete g4MatCmd;
164 delete g4DensCmd; 137 delete g4DensCmd;
165 delete densCmd; 138 delete densCmd;
166 delete adensCmd; 139 delete adensCmd;
167 delete fPosiCmd; <<
168 <<
169 delete g4Dir; 140 delete g4Dir;
170 delete matDir; << 141 delete matDir;
171 } 142 }
172 143
173 //....oooOO0OOooo........oooOO0OOooo........oo 144 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
174 145
175 void G4NistMessenger::SetNewValue(G4UIcommand* 146 void G4NistMessenger::SetNewValue(G4UIcommand* command, G4String newValue)
176 { 147 {
177 // G4cout << "G4NistMessenger::SetNewValue < << 148 //G4cout << "G4NistMessenger::SetNewValue <" << newValue << ">" << G4endl;
178 if (command == verCmd) { << 149 if (command == verCmd)
179 manager->SetVerbose(verCmd->GetNewIntValue << 150 { manager->SetVerbose(verCmd->GetNewIntValue(newValue));}
180 } << 151 else if (command == prtElmCmd)
181 else if (command == prtElmCmd) { << 152 { manager->PrintElement(newValue); }
182 manager->PrintElement(newValue); <<
183 } <<
184 else if (command == przElmCmd) { 153 else if (command == przElmCmd) {
185 G4int Z = przElmCmd->GetNewIntValue(newVal 154 G4int Z = przElmCmd->GetNewIntValue(newValue);
186 if (Z >= 0 && Z < 108) { << 155 if(Z >= 0 && Z < 108) { manager->PrintElement(Z); }
187 manager->PrintElement(Z); <<
188 } <<
189 } <<
190 else if (command == lisMatCmd) { <<
191 manager->ListMaterials(newValue); <<
192 } <<
193 else if (command == g4ElmCmd) { <<
194 manager->PrintG4Element(newValue); <<
195 } <<
196 else if (command == g4MatCmd) { <<
197 manager->PrintG4Material(newValue); <<
198 } <<
199 else if (command == g4DensCmd) { <<
200 G4IonisParamMat::GetDensityEffectData()->P <<
201 } <<
202 else if (command == densCmd) { <<
203 manager->SetDensityEffectCalculatorFlag(ne <<
204 } <<
205 else if (command == adensCmd) { <<
206 manager->SetDensityEffectCalculatorFlag(ne <<
207 } <<
208 else if (command == fPosiCmd) { <<
209 G4String mnam{""}; <<
210 G4double f{0.0}; <<
211 std::istringstream ss(newValue); <<
212 ss >> mnam >> f; <<
213 // set fraction for all materials <<
214 if (mnam == "all" || mnam == "none") { <<
215 if (mnam == "none" || f < 0.0) { f = 0.0 <<
216 auto mtable = G4Material::GetMaterialTab <<
217 for ( auto const & mat : *mtable ) { <<
218 mat->GetIonisation()->SetOrtoPositroniumFrac <<
219 } <<
220 } else { <<
221 // set fraction for one material <<
222 auto mat = manager->FindOrBuildMaterial( <<
223 if (nullptr != mat) { <<
224 mat->GetIonisation()->SetOrtoPositroniumFrac <<
225 } <<
226 } <<
227 } 156 }
>> 157 else if (command == lisMatCmd)
>> 158 { manager->ListMaterials(newValue); }
>> 159 else if (command == g4ElmCmd)
>> 160 { manager->PrintG4Element(newValue); }
>> 161 else if (command == g4MatCmd)
>> 162 { manager->PrintG4Material(newValue); }
>> 163 else if (command == g4DensCmd)
>> 164 { G4IonisParamMat::GetDensityEffectData()->PrintData(newValue); }
>> 165 else if (command == densCmd)
>> 166 { manager->SetDensityEffectCalculatorFlag(newValue,true); }
>> 167 else if (command == adensCmd)
>> 168 { manager->SetDensityEffectCalculatorFlag(newValue,false); }
228 } 169 }
>> 170
>> 171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
229 172