| 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 // $Id: G4HnMessenger.cc 66310 2012-12-17 11:56:35Z ihrivnac $
26 27
27 // Author: Ivana Hrivnacova, 18/06/2013 (ivan 28 // Author: Ivana Hrivnacova, 18/06/2013 (ivana@ipno.in2p3.fr)
28 29
29 #include "G4HnMessenger.hh" 30 #include "G4HnMessenger.hh"
30 #include "G4HnManager.hh" 31 #include "G4HnManager.hh"
31 #include "G4AnalysisUtilities.hh" 32 #include "G4AnalysisUtilities.hh"
>> 33 #include "G4AnalysisMessengerHelper.hh"
32 34
33 #include "G4UIcommand.hh" 35 #include "G4UIcommand.hh"
34 #include "G4UIparameter.hh" 36 #include "G4UIparameter.hh"
>> 37 #include "G4UIcmdWithAnInteger.hh"
35 #include "G4UIcmdWithABool.hh" 38 #include "G4UIcmdWithABool.hh"
36 #include "G4UIcmdWithAString.hh" <<
37 39
38 using namespace G4Analysis; 40 using namespace G4Analysis;
39 41
40 //____________________________________________ << 42 #include <iostream>
>> 43
41 G4HnMessenger::G4HnMessenger(G4HnManager& mana 44 G4HnMessenger::G4HnMessenger(G4HnManager& manager)
42 : fManager(manager), << 45 : G4UImessenger(),
43 fHnType(manager.GetHnType()), << 46 fManager(manager),
44 fHnDimension(std::stoi(manager.GetHnType() << 47 fHelper(nullptr),
45 { << 48 fSetHnAsciiCmd(nullptr),
>> 49 fSetHnActivationCmd(nullptr),
>> 50 fSetHnActivationAllCmd(nullptr),
>> 51 fSetHnPlottingCmd(nullptr),
>> 52 fSetHnPlottingAllCmd(nullptr)
>> 53 {
>> 54 G4String hnType = fManager.GetHnType();
>> 55 hnType.toLower();
>> 56 fHelper = G4Analysis::make_unique<G4AnalysisMessengerHelper>(hnType);
>> 57
46 SetHnAsciiCmd(); 58 SetHnAsciiCmd();
47 SetHnActivationCmd(); 59 SetHnActivationCmd();
48 SetHnActivationToAllCmd(); 60 SetHnActivationToAllCmd();
49 SetHnPlottingCmd(); 61 SetHnPlottingCmd();
50 SetHnPlottingToAllCmd(); 62 SetHnPlottingToAllCmd();
51 SetHnFileNameCmd(); <<
52 SetHnFileNameToAllCmd(); <<
53 <<
54 auto maxDim = (fHnDimension < kMaxDim) ? fHn <<
55 for (unsigned int idim = 0; idim < maxDim; + <<
56 fSetAxisLogCmd.push_back(CreateSetAxisLogC <<
57 } <<
58 } 63 }
59 64
60 //____________________________________________ 65 //_____________________________________________________________________________
61 G4HnMessenger::~G4HnMessenger() = default; << 66 G4HnMessenger::~G4HnMessenger()
>> 67 {}
62 68
63 // 69 //
64 // private functions 70 // private functions
65 // 71 //
66 72
67 //____________________________________________ 73 //_____________________________________________________________________________
68 G4String G4HnMessenger::GetObjectType() const <<
69 { <<
70 return (fHnType[0] == 'h') ? <<
71 fHnType.substr(1,1) + "D histogram" : <<
72 fHnType.substr(1,1) + "D profile"; <<
73 } <<
74 <<
75 //____________________________________________ <<
76 void G4HnMessenger::AddIdParameter(G4UIcommand <<
77 { <<
78 auto htId = new G4UIparameter("id", 'i', fal <<
79 htId->SetGuidance("Histogram id"); <<
80 htId->SetParameterRange("id>=0"); <<
81 command.SetParameter(htId); <<
82 } <<
83 <<
84 //____________________________________________ <<
85 void G4HnMessenger::AddOptionParameter(G4UIcom <<
86 { <<
87 auto param = new G4UIparameter(optionName, ' <<
88 auto guidance = GetObjectType() + " " + opti <<
89 param->SetGuidance(guidance.c_str()); <<
90 param->SetDefaultValue("true"); <<
91 command.SetParameter(param); <<
92 } <<
93 <<
94 //____________________________________________ <<
95 void G4HnMessenger::SetHnAsciiCmd() 74 void G4HnMessenger::SetHnAsciiCmd()
96 { 75 {
97 fSetAsciiCmd = << 76 fSetHnAsciiCmd
98 CreateCommand<G4UIcommand>("setAscii", "Pr << 77 = G4Analysis::make_unique<G4UIcmdWithAnInteger>(fHelper->Update("/analysis/HNTYPE_/setAscii"), this);
99 << 78 fSetHnAsciiCmd->SetGuidance(
100 AddIdParameter(*fSetAsciiCmd); << 79 fHelper->Update("Print NDIM_D LOBJECT of given id on ascii file."));
101 AddOptionParameter(*fSetAsciiCmd, "hnAscii") << 80
102 << 81 fSetHnAsciiCmd->SetParameterName("id",false);
>> 82 fSetHnAsciiCmd->SetRange("id>=0");
>> 83 fSetHnAsciiCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
103 } 84 }
104 85
105 //____________________________________________ 86 //_____________________________________________________________________________
106 void G4HnMessenger::SetHnActivationCmd() 87 void G4HnMessenger::SetHnActivationCmd()
107 { 88 {
108 fSetActivationCmd = << 89 auto hnId = new G4UIparameter("id", 'i', false);
109 CreateCommand<G4UIcommand>("setActivation" << 90 hnId->SetGuidance(fHelper->Update("OBJECT id"));
110 << 91 hnId->SetParameterRange("id>=0");
111 AddIdParameter(*fSetActivationCmd); << 92
112 AddOptionParameter(*fSetActivationCmd, "hnAc << 93 auto hnActivation = new G4UIparameter("hnActivation", 's', true);
113 } << 94 hnActivation->SetGuidance(fHelper->Update("OBJECT activation"));
>> 95 hnActivation->SetDefaultValue("none");
>> 96
>> 97 fSetHnActivationCmd
>> 98 = G4Analysis::make_unique<G4UIcommand>(fHelper->Update("/analysis/HNTYPE_/setActivation"), this);
>> 99 fSetHnActivationCmd->SetGuidance(
>> 100 fHelper->Update("Set activation for the NDIM_D LOBJECT of given id"));
>> 101 fSetHnActivationCmd->SetParameter(hnId);
>> 102 fSetHnActivationCmd->SetParameter(hnActivation);
>> 103 fSetHnActivationCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
>> 104 }
114 105
115 //____________________________________________ 106 //_____________________________________________________________________________
116 void G4HnMessenger::SetHnActivationToAllCmd() 107 void G4HnMessenger::SetHnActivationToAllCmd()
117 { 108 {
118 fSetActivationAllCmd = << 109 fSetHnActivationAllCmd
119 CreateCommand<G4UIcmdWithABool>( << 110 = G4Analysis::make_unique<G4UIcmdWithABool>(fHelper->Update("/analysis/HNTYPE_/setActivationToAll"), this);
120 "setActivationToAll", "Set activation to << 111 fSetHnActivationAllCmd->SetGuidance(
121 fSetActivationAllCmd->SetParameterName("Acti << 112 fHelper->Update("Set activation to all NDIM_D LOBJECTs"));
122 } << 113 fSetHnActivationAllCmd->SetParameterName("Activation",false);
123 << 114 }
>> 115
124 //____________________________________________ 116 //_____________________________________________________________________________
125 void G4HnMessenger::SetHnPlottingCmd() 117 void G4HnMessenger::SetHnPlottingCmd()
126 { 118 {
127 fSetPlottingCmd = << 119 auto hnId = new G4UIparameter("id", 'i', false);
128 CreateCommand<G4UIcommand>("setPlotting", << 120 hnId->SetGuidance(fHelper->Update("OBJECT id"));
129 << 121 hnId->SetParameterRange("id>=0");
130 AddIdParameter(*fSetPlottingCmd); << 122
131 AddOptionParameter(*fSetPlottingCmd, "hnPlot << 123 auto hnPlotting = new G4UIparameter("hnPlotting", 's', true);
132 } << 124 hnPlotting->SetGuidance(fHelper->Update("(In)Activate OBJECT plotting"));
>> 125 hnPlotting->SetDefaultValue("none");
>> 126
>> 127 fSetHnPlottingCmd
>> 128 = G4Analysis::make_unique<G4UIcommand>(fHelper->Update("/analysis/HNTYPE_/setPlotting"), this);
>> 129 fSetHnPlottingCmd->SetGuidance(
>> 130 fHelper->Update("(In)Activate batch plotting of the NDIM_D LOBJECT of given id"));
>> 131 fSetHnPlottingCmd->SetParameter(hnId);
>> 132 fSetHnPlottingCmd->SetParameter(hnPlotting);
>> 133 fSetHnPlottingCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
>> 134 }
133 135
134 //____________________________________________ 136 //_____________________________________________________________________________
135 void G4HnMessenger::SetHnPlottingToAllCmd() 137 void G4HnMessenger::SetHnPlottingToAllCmd()
136 { 138 {
137 fSetPlottingAllCmd = << 139 fSetHnPlottingAllCmd
138 CreateCommand<G4UIcmdWithABool>( << 140 = G4Analysis::make_unique<G4UIcmdWithABool>(fHelper->Update("/analysis/HNTYPE_/setPlottingToAll"), this);
139 "setPlottingToAll", "(In)Activate batch << 141 fSetHnPlottingAllCmd->SetGuidance(
140 fSetPlottingAllCmd->SetParameterName("Plotti << 142 fHelper->Update("(In)Activate batch plotting of all NDIM_D LOBJECTs"));
141 } << 143 fSetHnPlottingAllCmd->SetParameterName("Plotting",false);
142 << 144 }
143 //____________________________________________ <<
144 void G4HnMessenger::SetHnFileNameCmd() <<
145 { <<
146 fSetFileNameCmd = <<
147 CreateCommand<G4UIcommand>("setFileName", <<
148 <<
149 AddIdParameter(*fSetFileNameCmd); <<
150 <<
151 auto param = new G4UIparameter("hnFileName", <<
152 auto guidance = GetObjectType() + " output f <<
153 param->SetGuidance(guidance.c_str()); <<
154 fSetFileNameCmd->SetParameter(param); <<
155 } <<
156 <<
157 //____________________________________________ <<
158 void G4HnMessenger::SetHnFileNameToAllCmd() <<
159 { <<
160 fSetFileNameAllCmd = <<
161 CreateCommand<G4UIcmdWithAString>( <<
162 "setFileNameToAll", "Set output file nam <<
163 fSetFileNameAllCmd->SetParameterName("FileN <<
164 } <<
165 <<
166 //____________________________________________ <<
167 std::unique_ptr<G4UIcommand> <<
168 G4HnMessenger::CreateSetAxisLogCommand(unsigne <<
169 { <<
170 G4String xyz{"XYZ"}; <<
171 auto axis = xyz.substr(idim, 1); <<
172 <<
173 G4String commandName = "set" + axis + "axisL <<
174 G4String guidance = "Activate " + axis + "- <<
175 <<
176 auto command = CreateCommand<G4UIcommand>(st <<
177 command->AvailableForStates(G4State_PreInit, <<
178 <<
179 // Add Id parameter <<
180 AddIdParameter(*command); <<
181 <<
182 auto parAxisLog = new G4UIparameter("axis", <<
183 guidance = GetObjectType() + " " + std::move <<
184 parAxisLog->SetGuidance(guidance.c_str()); <<
185 command->SetParameter(parAxisLog); <<
186 <<
187 return command; <<
188 } <<
189 145
190 // 146 //
191 // public methods 147 // public methods
192 // 148 //
193 149
194 //____________________________________________ 150 //_____________________________________________________________________________
195 void G4HnMessenger::SetNewValue(G4UIcommand* c 151 void G4HnMessenger::SetNewValue(G4UIcommand* command, G4String newValues)
196 { 152 {
197 // process "All" commands first << 153 if ( command == fSetHnAsciiCmd.get() ) {
198 if ( command == fSetActivationAllCmd.get() ) << 154 auto id = fSetHnAsciiCmd->GetNewIntValue(newValues);
199 fManager.SetActivation(fSetActivationAllCm << 155 fManager.SetAscii(id, true);
200 return; << 156 }
201 } << 157 else if ( command == fSetHnActivationCmd.get() ) {
202 << 158 // tokenize parameters in a vector
203 if ( command == fSetPlottingAllCmd.get() ) { << 159 std::vector<G4String> parameters;
204 fManager.SetPlotting(fSetPlottingAllCmd->G << 160 G4Analysis::Tokenize(newValues, parameters);
205 return; << 161 // check consistency
206 } << 162 if ( G4int(parameters.size()) == command->GetParameterEntries() ) {
207 << 163 auto counter = 0;
208 if ( command == fSetFileNameAllCmd.get() ) { << 164 auto id = G4UIcommand::ConvertToInt(parameters[counter++]);
209 fManager.SetFileName(newValues); << 165 auto activation = G4UIcommand::ConvertToBool(parameters[counter++]);
210 return; << 166 fManager.SetActivation(id, activation);
211 } <<
212 <<
213 // Tokenize parameters in a vector <<
214 std::vector<G4String> parameters; <<
215 G4Analysis::Tokenize(newValues, parameters); <<
216 // check consistency <<
217 if ( parameters.size() != command->GetParame <<
218 // Should never happen but let's check any <<
219 G4Analysis::Warn( <<
220 "Got wrong number of \"" + command->GetC <<
221 "\" parameters: " + std::to_string(param <<
222 " instead of " + std::to_string(command- <<
223 fkClass, "WarnAboutParameters"); <<
224 return; <<
225 } <<
226 <<
227 auto counter = 0; <<
228 auto id = G4UIcommand::ConvertToInt(paramete <<
229 <<
230 if ( command == fSetAsciiCmd.get() ) { <<
231 fManager.SetAscii(id, G4UIcommand::Convert <<
232 return; <<
233 } <<
234 <<
235 if ( command == fSetActivationCmd.get() ) { <<
236 fManager.SetActivation(id, G4UIcommand::Co <<
237 return; <<
238 } <<
239 <<
240 if ( command == fSetPlottingCmd.get() ) { <<
241 fManager.SetPlotting(id, G4UIcommand::Conv <<
242 return; <<
243 } <<
244 <<
245 if ( command == fSetFileNameCmd.get() ) { <<
246 fManager.SetFileName(id, parameters[counte <<
247 return; <<
248 } <<
249 <<
250 auto maxDim = (fHnDimension < kMaxDim) ? fHn <<
251 for (unsigned int idim = 0; idim < maxDim; + <<
252 if ( command == fSetAxisLogCmd[idim].get() <<
253 auto axisLog = G4UIcommand::ConvertToBoo <<
254 fManager.SetAxisIsLog(idim, id, axisLog) <<
255 return; <<
256 } 167 }
257 } << 168 else {
258 } << 169 // Should never happen but let's check anyway for consistency
>> 170 fHelper->WarnAboutParameters(command, parameters.size());
>> 171 }
>> 172 }
>> 173 else if ( command == fSetHnActivationAllCmd.get() ) {
>> 174 auto activation = fSetHnActivationAllCmd->GetNewBoolValue(newValues);
>> 175 fManager.SetActivation(activation);
>> 176 }
>> 177 else if ( command == fSetHnPlottingCmd.get() ) {
>> 178 // tokenize parameters in a vector
>> 179 std::vector<G4String> parameters;
>> 180 G4Analysis::Tokenize(newValues, parameters);
>> 181 // check consistency
>> 182 if ( G4int(parameters.size()) == command->GetParameterEntries() ) {
>> 183 auto counter = 0;
>> 184 auto id = G4UIcommand::ConvertToInt(parameters[counter++]);
>> 185 auto activation = G4UIcommand::ConvertToBool(parameters[counter++]);
>> 186 fManager.SetPlotting(id, activation);
>> 187 }
>> 188 else {
>> 189 // Should never happen but let's check anyway for consistency
>> 190 fHelper->WarnAboutParameters(command, parameters.size());
>> 191 }
>> 192 }
>> 193 else if ( command == fSetHnPlottingAllCmd.get() ) {
>> 194 auto activation = fSetHnPlottingAllCmd->GetNewBoolValue(newValues);
>> 195 fManager.SetPlotting(activation);
>> 196 }
>> 197 }
259 198