| 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 32
33 #include "G4UIcommand.hh" 33 #include "G4UIcommand.hh"
34 #include "G4UIparameter.hh" 34 #include "G4UIparameter.hh"
>> 35 #include "G4UIcmdWithAnInteger.hh"
35 #include "G4UIcmdWithABool.hh" 36 #include "G4UIcmdWithABool.hh"
36 #include "G4UIcmdWithAString.hh" <<
37 37
38 using namespace G4Analysis; << 38 #include <iostream>
39 39
40 //____________________________________________ 40 //_____________________________________________________________________________
41 G4HnMessenger::G4HnMessenger(G4HnManager& mana << 41 G4HnMessenger::G4HnMessenger(G4HnManager* manager)
42 : fManager(manager), << 42 : G4UImessenger(),
43 fHnType(manager.GetHnType()), << 43 fManager(manager),
44 fHnDimension(std::stoi(manager.GetHnType() << 44 fSetHnAsciiCmd(0),
45 { << 45 fSetHnActivationCmd(0),
>> 46 fSetHnActivationAllCmd(0)
>> 47 {
>> 48
46 SetHnAsciiCmd(); 49 SetHnAsciiCmd();
47 SetHnActivationCmd(); 50 SetHnActivationCmd();
48 SetHnActivationToAllCmd(); 51 SetHnActivationToAllCmd();
49 SetHnPlottingCmd(); <<
50 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 } 52 }
59 53
60 //____________________________________________ 54 //_____________________________________________________________________________
61 G4HnMessenger::~G4HnMessenger() = default; << 55 G4HnMessenger::~G4HnMessenger()
>> 56 {
>> 57 delete fSetHnAsciiCmd;
>> 58 delete fSetHnActivationCmd;
>> 59 delete fSetHnActivationAllCmd;
>> 60 }
62 61
63 // 62 //
64 // private functions 63 // private functions
65 // 64 //
66 65
67 //____________________________________________ 66 //_____________________________________________________________________________
68 G4String G4HnMessenger::GetObjectType() const << 67 G4String G4HnMessenger::GetCmdDirectoryName() 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 { 68 {
78 auto htId = new G4UIparameter("id", 'i', fal << 69 G4String hnType = fManager->GetHnType();
79 htId->SetGuidance("Histogram id"); << 70 hnType.toLower();
80 htId->SetParameterRange("id>=0"); << 71 G4String name("/analysis/");
81 command.SetParameter(htId); << 72 name.append(hnType);
>> 73 name.append("/");
>> 74 return name;
82 } 75 }
83 76
84 //____________________________________________ 77 //_____________________________________________________________________________
85 void G4HnMessenger::AddOptionParameter(G4UIcom << 78 G4String G4HnMessenger::GetHnDescription() const
86 { 79 {
87 auto param = new G4UIparameter(optionName, ' << 80 G4String hnDescription = fManager->GetHnType();
88 auto guidance = GetObjectType() + " " + opti << 81 hnDescription.erase(0);
89 param->SetGuidance(guidance.c_str()); << 82 hnDescription.append("D");
90 param->SetDefaultValue("true"); << 83 return hnDescription;
91 command.SetParameter(param); <<
92 } 84 }
93 85
94 //____________________________________________ 86 //_____________________________________________________________________________
95 void G4HnMessenger::SetHnAsciiCmd() 87 void G4HnMessenger::SetHnAsciiCmd()
96 { 88 {
97 fSetAsciiCmd = << 89 G4String name = GetCmdDirectoryName();
98 CreateCommand<G4UIcommand>("setAscii", "Pr << 90 name.append("setAscii");
99 << 91 fSetHnAsciiCmd = new G4UIcmdWithAnInteger(name, this);
100 AddIdParameter(*fSetAsciiCmd); << 92
101 AddOptionParameter(*fSetAsciiCmd, "hnAscii") << 93 G4String guidance("Print ");
102 << 94 guidance.append(GetHnDescription());
>> 95 guidance.append(" histogram of #Id on ascii file.");
>> 96
>> 97 fSetHnAsciiCmd->SetGuidance(guidance);
>> 98 fSetHnAsciiCmd->SetParameterName("Id",false);
>> 99 fSetHnAsciiCmd->SetRange("Id>=0");
>> 100 fSetHnAsciiCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
103 } 101 }
104 102
105 //____________________________________________ 103 //_____________________________________________________________________________
106 void G4HnMessenger::SetHnActivationCmd() 104 void G4HnMessenger::SetHnActivationCmd()
107 { 105 {
108 fSetActivationCmd = << 106 G4String name = GetCmdDirectoryName();
109 CreateCommand<G4UIcommand>("setActivation" << 107 name.append("setActivationToAll");
110 << 108 fSetHnAsciiCmd = new G4UIcmdWithAnInteger(name, this);
111 AddIdParameter(*fSetActivationCmd); << 109
112 AddOptionParameter(*fSetActivationCmd, "hnAc << 110 G4String guidance("Set activation to all ");
113 } << 111 guidance.append(GetHnDescription());
114 << 112 guidance.append(" histograms.");
>> 113
>> 114 fSetHnActivationAllCmd = new G4UIcmdWithABool(name, this);
>> 115 fSetHnActivationAllCmd->SetGuidance(guidance);
>> 116 fSetHnActivationAllCmd->SetParameterName("Activation",false);
>> 117 }
>> 118
115 //____________________________________________ 119 //_____________________________________________________________________________
116 void G4HnMessenger::SetHnActivationToAllCmd() 120 void G4HnMessenger::SetHnActivationToAllCmd()
117 { 121 {
118 fSetActivationAllCmd = << 122 G4UIparameter* hnId = new G4UIparameter("idActivation", 'i', false);
119 CreateCommand<G4UIcmdWithABool>( << 123 hnId->SetGuidance("Histogram id");
120 "setActivationToAll", "Set activation to << 124 hnId->SetParameterRange("idActivation>=0");
121 fSetActivationAllCmd->SetParameterName("Acti << 125
122 } << 126 G4UIparameter* hnActivation = new G4UIparameter("hnActivation", 's', true);
123 << 127 hnActivation->SetGuidance("Histogram activation");
124 //____________________________________________ << 128 hnActivation->SetDefaultValue("none");
125 void G4HnMessenger::SetHnPlottingCmd() << 129
126 { << 130 G4String name = GetCmdDirectoryName();
127 fSetPlottingCmd = << 131 name.append("setActivation");
128 CreateCommand<G4UIcommand>("setPlotting", << 132 fSetHnActivationCmd = new G4UIcommand(name, this);
129 << 133
130 AddIdParameter(*fSetPlottingCmd); << 134 G4String guidance("Set activation for the ");
131 AddOptionParameter(*fSetPlottingCmd, "hnPlot << 135 guidance.append(GetHnDescription());
132 } << 136 guidance.append(" histogram of #Id");
133 << 137
134 //____________________________________________ << 138 fSetHnActivationCmd->SetGuidance(guidance);
135 void G4HnMessenger::SetHnPlottingToAllCmd() << 139 fSetHnActivationCmd->SetParameter(hnId);
136 { << 140 fSetHnActivationCmd->SetParameter(hnActivation);
137 fSetPlottingAllCmd = << 141 fSetHnActivationCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
138 CreateCommand<G4UIcmdWithABool>( << 142 }
139 "setPlottingToAll", "(In)Activate batch <<
140 fSetPlottingAllCmd->SetParameterName("Plotti <<
141 } <<
142 <<
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 143
190 // 144 //
191 // public methods 145 // public methods
192 // 146 //
193 147
194 //____________________________________________ 148 //_____________________________________________________________________________
195 void G4HnMessenger::SetNewValue(G4UIcommand* c 149 void G4HnMessenger::SetNewValue(G4UIcommand* command, G4String newValues)
196 { 150 {
197 // process "All" commands first << 151 if ( command == fSetHnAsciiCmd ) {
198 if ( command == fSetActivationAllCmd.get() ) << 152 G4int id = fSetHnAsciiCmd->GetNewIntValue(newValues);
199 fManager.SetActivation(fSetActivationAllCm << 153 fManager->SetAscii(id, true);
200 return; << 154 }
201 } << 155 else if ( command == fSetHnActivationCmd ) {
202 << 156 G4int id;
203 if ( command == fSetPlottingAllCmd.get() ) { << 157 G4String sactivation;
204 fManager.SetPlotting(fSetPlottingAllCmd->G << 158 std::istringstream is(newValues.data());
205 return; << 159 is >> id >> sactivation;
206 } << 160 G4bool activation = G4UIcommand::ConvertToBool(sactivation);
207 << 161 fManager->SetActivation(id, activation);
208 if ( command == fSetFileNameAllCmd.get() ) { << 162 }
209 fManager.SetFileName(newValues); << 163 else if ( command == fSetHnActivationAllCmd ) {
210 return; << 164 G4bool activation = fSetHnActivationAllCmd->GetNewBoolValue(newValues);
211 } << 165 fManager->SetActivation(activation);
212 << 166 }
213 // Tokenize parameters in a vector << 167 }
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 } <<
257 } <<
258 } <<
259 168