| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25
26 // File name: G4NistMessenger
27 //
28 // Author: Vladimir Ivanchenko
29 //
30 // Creation date: 23.12.2004
31 //
32 // Modifications:
33 //
34 //
35 // -------------------------------------------------------------------
36
37 #include "G4NistMessenger.hh"
38
39 #include "G4ApplicationState.hh"
40 #include "G4DensityEffectData.hh"
41 #include "G4IonisParamMat.hh"
42 #include "G4NistManager.hh"
43 #include "G4UIcmdWithAString.hh"
44 #include "G4UIcmdWithAnInteger.hh"
45 #include "G4UIdirectory.hh"
46 #include "G4UIcommand.hh"
47 #include "G4UIparameter.hh"
48
49 #include <sstream>
50
51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
52
53 G4NistMessenger::G4NistMessenger(G4NistManager* man) : manager(man)
54 {
55 matDir = new G4UIdirectory("/material/");
56 matDir->SetGuidance("Commands for materials");
57
58 verCmd = new G4UIcmdWithAnInteger("/material/verbose", this);
59 verCmd->SetGuidance("Set verbose level.");
60 verCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
61 verCmd->SetToBeBroadcasted(false);
62
63 nistDir = new G4UIdirectory("/material/nist/");
64 nistDir->SetGuidance("Commands for the nist dataBase");
65
66 prtElmCmd = new G4UIcmdWithAString("/material/nist/printElement", this);
67 prtElmCmd->SetGuidance("print element(s) in dataBase.");
68 prtElmCmd->SetGuidance("symbol = element.");
69 prtElmCmd->SetGuidance("all = all elements.");
70 prtElmCmd->SetParameterName("symbol", true);
71 prtElmCmd->SetDefaultValue("all");
72 prtElmCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
73 prtElmCmd->SetToBeBroadcasted(false);
74
75 przElmCmd = new G4UIcmdWithAnInteger("/material/nist/printElementZ", this);
76 przElmCmd->SetGuidance("print element Z in dataBase.");
77 przElmCmd->SetGuidance("0 = all elements.");
78 przElmCmd->SetParameterName("Z", true);
79 przElmCmd->SetDefaultValue(0);
80 przElmCmd->SetRange("0<=Z && Z<108");
81 przElmCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
82 przElmCmd->SetToBeBroadcasted(false);
83
84 lisMatCmd = new G4UIcmdWithAString("/material/nist/listMaterials", this);
85 lisMatCmd->SetGuidance("Materials in Geant4 dataBase.");
86 lisMatCmd->SetGuidance("simple - simple NIST materials.");
87 lisMatCmd->SetGuidance("compound - compound NIST materials.");
88 lisMatCmd->SetGuidance("hep - HEP materials.");
89 lisMatCmd->SetGuidance("bio - biomedical materials.");
90 lisMatCmd->SetGuidance("all - list of all Geant4 materials.");
91 lisMatCmd->SetParameterName("matlist", true);
92 // lisMatCmd->SetCandidates("simple compound hep bio all");
93 lisMatCmd->SetDefaultValue("all");
94 lisMatCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
95 lisMatCmd->SetToBeBroadcasted(false);
96
97 g4Dir = new G4UIdirectory("/material/g4/");
98 g4Dir->SetGuidance("Commands for G4MaterialTable");
99
100 g4ElmCmd = new G4UIcmdWithAString("/material/g4/printElement", this);
101 g4ElmCmd->SetGuidance("print Element from G4ElementTable.");
102 g4ElmCmd->SetGuidance("all - all elements.");
103 g4ElmCmd->SetParameterName("elm", true);
104 g4ElmCmd->SetDefaultValue("all");
105
106 g4MatCmd = new G4UIcmdWithAString("/material/g4/printMaterial", this);
107 g4MatCmd->SetGuidance("print Material from G4MaterialTable.");
108 g4MatCmd->SetGuidance("all - all materials");
109 g4MatCmd->SetParameterName("pmat", true);
110 g4MatCmd->SetDefaultValue("all");
111 g4MatCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
112 g4MatCmd->SetToBeBroadcasted(false);
113
114 g4DensCmd = new G4UIcmdWithAString("/material/g4/printDensityEffParam", this);
115 g4DensCmd->SetGuidance("print Material from G4DensityEffectData.");
116 g4DensCmd->SetGuidance("all - all materials");
117 g4DensCmd->SetParameterName("dmat", true);
118 g4DensCmd->SetDefaultValue("all");
119 g4DensCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
120 g4DensCmd->SetToBeBroadcasted(false);
121
122 densCmd = new G4UIcmdWithAString("/material/g4/enableDensityEffOnFly", this);
123 densCmd->SetGuidance("enable accurate computation of density effect.");
124 densCmd->SetGuidance("all - all materials.");
125 densCmd->SetParameterName("dens", true);
126 densCmd->SetDefaultValue("all");
127 densCmd->AvailableForStates(G4State_PreInit);
128 lisMatCmd->SetToBeBroadcasted(false);
129
130 adensCmd = new G4UIcmdWithAString("/material/g4/disableDensityEffOnFly", this);
131 adensCmd->SetGuidance("disable accurate computation of density effect.");
132 adensCmd->SetGuidance("all - all materials.");
133 adensCmd->SetParameterName("dens", true);
134 adensCmd->SetDefaultValue("all");
135 adensCmd->AvailableForStates(G4State_PreInit);
136 adensCmd->SetToBeBroadcasted(false);
137
138 fPosiCmd = new G4UIcommand("/material/g4/ortoPositroniumFraction", this);
139 fPosiCmd->SetGuidance("defined orto-positronium fraction for positron annihilation AtRest.");
140 fPosiCmd->SetGuidance("via material name, all - all materials.");
141 fPosiCmd->AvailableForStates(G4State_PreInit);
142 fPosiCmd->SetToBeBroadcasted(false);
143
144 auto p1 = new G4UIparameter("matname",'s', false);
145 fPosiCmd->SetParameter(p1);
146
147 auto p2 = new G4UIparameter("fraction", 'd', false);
148 p2->SetParameterRange("fraction>=0. && fraction <=1.");
149 fPosiCmd->SetParameter(p2);
150 }
151
152 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
153
154 G4NistMessenger::~G4NistMessenger()
155 {
156 delete verCmd;
157 delete prtElmCmd;
158 delete przElmCmd;
159 delete lisMatCmd;
160 delete nistDir;
161
162 delete g4ElmCmd;
163 delete g4MatCmd;
164 delete g4DensCmd;
165 delete densCmd;
166 delete adensCmd;
167 delete fPosiCmd;
168
169 delete g4Dir;
170 delete matDir;
171 }
172
173 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
174
175 void G4NistMessenger::SetNewValue(G4UIcommand* command, G4String newValue)
176 {
177 // G4cout << "G4NistMessenger::SetNewValue <" << newValue << ">" << G4endl;
178 if (command == verCmd) {
179 manager->SetVerbose(verCmd->GetNewIntValue(newValue));
180 }
181 else if (command == prtElmCmd) {
182 manager->PrintElement(newValue);
183 }
184 else if (command == przElmCmd) {
185 G4int Z = przElmCmd->GetNewIntValue(newValue);
186 if (Z >= 0 && Z < 108) {
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()->PrintData(newValue);
201 }
202 else if (command == densCmd) {
203 manager->SetDensityEffectCalculatorFlag(newValue, true);
204 }
205 else if (command == adensCmd) {
206 manager->SetDensityEffectCalculatorFlag(newValue, false);
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::GetMaterialTable();
217 for ( auto const & mat : *mtable ) {
218 mat->GetIonisation()->SetOrtoPositroniumFraction(f);
219 }
220 } else {
221 // set fraction for one material
222 auto mat = manager->FindOrBuildMaterial(mnam, true);
223 if (nullptr != mat) {
224 mat->GetIonisation()->SetOrtoPositroniumFraction(f);
225 }
226 }
227 }
228 }
229