Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/event/src/G4GeneralParticleSourceMessenger.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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 // G4GeneralParticleSourceMessenger class implementation
 27 //
 28 // Author: Fan Lei, QinetiQ ltd.
 29 // Customer: ESA/ESTEC
 30 // Revisions: Andrew Green, Andrea Dotti
 31 // --------------------------------------------------------------------
 32 
 33 #include "G4GeneralParticleSourceMessenger.hh"
 34 
 35 #include "G4PhysicalConstants.hh"
 36 #include "G4SystemOfUnits.hh"
 37 #include "G4Geantino.hh"
 38 #include "G4ThreeVector.hh"
 39 #include "G4ParticleTable.hh"
 40 #include "G4IonTable.hh"
 41 #include "G4UIdirectory.hh"
 42 #include "G4UIcmdWithoutParameter.hh"
 43 #include "G4UIcmdWithAString.hh"
 44 #include "G4UIcmdWithADoubleAndUnit.hh"
 45 #include "G4UIcmdWith3Vector.hh"
 46 #include "G4UIcmdWith3VectorAndUnit.hh"
 47 #include "G4UIcmdWithAnInteger.hh"
 48 #include "G4UIcmdWithADouble.hh"
 49 #include "G4UIcmdWithABool.hh"
 50 #include "G4ios.hh"
 51 
 52 #include "G4Tokenizer.hh"
 53 #include "G4SingleParticleSource.hh"
 54 #include "G4GeneralParticleSource.hh"
 55 
 56 #include "G4AutoLock.hh"
 57 
 58 namespace
 59 {
 60   G4Mutex creationM = G4MUTEX_INITIALIZER;
 61   G4GeneralParticleSourceMessenger* theInstance = nullptr;
 62 }
 63 
 64 G4GeneralParticleSourceMessenger*
 65 G4GeneralParticleSourceMessenger::GetInstance(G4GeneralParticleSource* psc)
 66 {
 67   G4AutoLock l(&creationM);
 68   if ( theInstance == nullptr )
 69     theInstance = new G4GeneralParticleSourceMessenger(psc);
 70   return theInstance;
 71 }
 72 
 73 void G4GeneralParticleSourceMessenger::Destroy()
 74 {
 75   G4AutoLock l(&creationM);
 76   if ( theInstance != nullptr )
 77   {
 78     delete theInstance;
 79     theInstance = nullptr;
 80   }
 81 }
 82 
 83 // --------------------------------------------------------------------
 84 //
 85 G4GeneralParticleSourceMessenger::
 86 G4GeneralParticleSourceMessenger(G4GeneralParticleSource* fPtclGun) 
 87   : fGPS(fPtclGun)
 88 {
 89   // A.Dotti - 10th October 2014
 90   // This messenger is special: it is instantiated in a user action but (e.g.
 91   // in a thread).
 92   // The UI commands it defines should be executed by the *master* thread
 93   // because they operate on shared resources and we want the UI commands to
 94   // take effect BEFORE the threads do some work (so all data are properly
 95   // initialized).
 96   // To achieve this behavior we set to true a base class protected
 97   // data member. Since it makes no sense to have more than one instance
 98   // of the messenger, we check that we actually have only one.
 99   // Note that the logic of implementing, in a given worker thread only one
100   // messenger is deleted/fated to the creator
101 
102   commandsShouldBeInMaster = true;
103 
104   particleTable = G4ParticleTable::GetParticleTable();
105   histtype = "biasx";
106 
107   // UI Commands only for master
108   //
109   G4bool broadcast = false;
110   gpsDirectory = new G4UIdirectory("/gps/",broadcast);
111     
112   gpsDirectory->SetGuidance("General Particle Source control commands.");
113 
114   // Now the commands for multiple sources
115   //
116   sourceDirectory = new G4UIdirectory("/gps/source/");
117   sourceDirectory->SetGuidance("Multiple source control sub-directory");
118 
119   addsourceCmd = new G4UIcmdWithADouble("/gps/source/add",this);
120   addsourceCmd->SetGuidance("Add a new source definition to the particle gun");
121   addsourceCmd->SetGuidance(" with the specified intensity");
122   addsourceCmd->SetParameterName("addsource",false,false);
123   addsourceCmd->SetRange("addsource > 0.");
124 
125   listsourceCmd = new G4UIcmdWithoutParameter("/gps/source/list",this);
126   listsourceCmd->SetGuidance("List the defined particle sources");
127 
128   clearsourceCmd = new G4UIcmdWithoutParameter("/gps/source/clear",this);
129   clearsourceCmd->SetGuidance("Remove all the defined particle sources");
130 
131   getsourceCmd = new G4UIcmdWithoutParameter("/gps/source/show",this);
132   getsourceCmd->SetGuidance("Show the current source index and intensity");
133 
134   setsourceCmd = new G4UIcmdWithAnInteger("/gps/source/set",this);
135   setsourceCmd->SetGuidance("Set the indexed source as the current one");
136   setsourceCmd->SetGuidance(" so one can change its source definition");
137   setsourceCmd->SetParameterName("setsource",false,false);
138   setsourceCmd->SetRange("setsource >= 0");
139 
140   deletesourceCmd = new G4UIcmdWithAnInteger("/gps/source/delete",this);
141   deletesourceCmd->SetGuidance("Delete the indexed source from the list");
142   deletesourceCmd->SetParameterName("deletesource",false,false);
143   deletesourceCmd->SetRange("deletesource > 0");
144 
145   setintensityCmd = new G4UIcmdWithADouble("/gps/source/intensity",this);
146   setintensityCmd->SetGuidance("Reset the current source to the specified intensity");
147   setintensityCmd->SetParameterName("setintensity",false,false);
148   setintensityCmd->SetRange("setintensity > 0."); 
149 
150   multiplevertexCmd = new G4UIcmdWithABool("/gps/source/multiplevertex",this);
151   multiplevertexCmd->SetGuidance("True for simultaneous generation multiple vertex");
152   multiplevertexCmd->SetGuidance(" Default is false");
153   multiplevertexCmd->SetParameterName("multiplevertex",true);
154   multiplevertexCmd->SetDefaultValue(false);
155 
156   flatsamplingCmd = new G4UIcmdWithABool("/gps/source/flatsampling",this);
157   flatsamplingCmd->SetGuidance("True for applying flat (biased) sampling among the sources");
158   flatsamplingCmd->SetGuidance("Default is false");
159   flatsamplingCmd->SetParameterName("flatsampling",true);
160   flatsamplingCmd->SetDefaultValue(false);
161 
162   // Below we reproduce commands awailable in G4Particle Gun
163   //
164   listCmd = new G4UIcmdWithoutParameter("/gps/List",this);
165   listCmd->SetGuidance("List available particles.");
166   listCmd->SetGuidance(" Invoke G4ParticleTable.");
167 
168   particleCmd = new G4UIcmdWithAString("/gps/particle",this);
169   particleCmd->SetGuidance("Set particle to be generated.");
170   particleCmd->SetGuidance(" (geantino is default)");
171   particleCmd->SetGuidance(" (ion can be specified for shooting ions)");
172   particleCmd->SetParameterName("particleName",true);
173   particleCmd->SetDefaultValue("geantino");
174   G4String candidateList; 
175   G4int nPtcl = particleTable->entries();
176   for(G4int i=0; i<nPtcl; ++i)
177   {
178     candidateList += particleTable->GetParticleName(i);
179     candidateList += " ";
180   }
181   candidateList += "ion ";
182   particleCmd->SetCandidates(candidateList);
183 
184   directionCmd = new G4UIcmdWith3Vector("/gps/direction",this);
185   directionCmd->SetGuidance("Set momentum direction.");
186   directionCmd->SetGuidance(" Direction needs not to be a unit vector.");
187   directionCmd->SetGuidance(" Angular distribution type is set to planar.");
188   directionCmd->SetParameterName("Px","Py","Pz",false,false); 
189   directionCmd->SetRange("Px != 0 || Py != 0 || Pz != 0");
190   
191   energyCmd = new G4UIcmdWithADoubleAndUnit("/gps/energy",this);
192   energyCmd->SetGuidance("Set kinetic energy.");
193   energyCmd->SetParameterName("Energy",false,false);
194   energyCmd->SetDefaultUnit("GeV");
195   //energyCmd->SetUnitCategory("Energy");
196   //energyCmd->SetUnitCandidates("eV keV MeV GeV TeV");
197 
198   positionCmd = new G4UIcmdWith3VectorAndUnit("/gps/position",this);
199   positionCmd->SetGuidance("Set starting position of the particle for a Point like source.");
200   positionCmd->SetGuidance(" Same effect as the two /gps/pos/type Point /gps/pos/centre commands.");
201   positionCmd->SetParameterName("X","Y","Z",false,false);
202   positionCmd->SetDefaultUnit("cm");
203   //positionCmd->SetUnitCategory("Length");
204   //positionCmd->SetUnitCandidates("microm mm cm m km");
205 
206   ionCmd = new G4UIcommand("/gps/ion",this);
207   ionCmd->SetGuidance("Set properties of ion to be generated.");
208   ionCmd->SetGuidance("[usage] /gps/ion Z A Q E");
209   ionCmd->SetGuidance("        Z:(int) AtomicNumber");
210   ionCmd->SetGuidance("        A:(int) AtomicMass");
211   ionCmd->SetGuidance("        Q:(int) Charge of Ion (in unit of e)");
212   ionCmd->SetGuidance("        E:(double) Excitation energy (in keV)");
213 
214   G4UIparameter* param;
215   param = new G4UIparameter("Z",'i',false);
216   param->SetDefaultValue("1");
217   ionCmd->SetParameter(param);
218   param = new G4UIparameter("A",'i',false);
219   param->SetDefaultValue("1");
220   ionCmd->SetParameter(param);
221   param = new G4UIparameter("Q",'i',true);
222   param->SetDefaultValue("0");
223   ionCmd->SetParameter(param);
224   param = new G4UIparameter("E",'d',true);
225   param->SetDefaultValue("0.0");
226   ionCmd->SetParameter(param);
227 
228   ionLvlCmd = new G4UIcommand("/gps/ionLvl",this);
229   ionLvlCmd->SetGuidance("Set properties of ion to be generated.");
230   ionLvlCmd->SetGuidance("[usage] /gps/ion Z A Q Lvl");
231   ionLvlCmd->SetGuidance("        Z:(int) AtomicNumber");
232   ionLvlCmd->SetGuidance("        A:(int) AtomicMass");
233   ionLvlCmd->SetGuidance("        Q:(int) Charge of Ion (in unit of e)");
234   ionLvlCmd->SetGuidance("        Lvl:(int) Number of metastable state excitation level (0-9)");
235 
236   G4UIparameter* paramL;
237   paramL = new G4UIparameter("Z",'i',false);
238   paramL->SetDefaultValue("1");
239   ionLvlCmd->SetParameter(paramL);
240   paramL = new G4UIparameter("A",'i',false);
241   paramL->SetDefaultValue("1");
242   ionLvlCmd->SetParameter(paramL);
243   paramL = new G4UIparameter("Q",'i',true);
244   paramL->SetDefaultValue("0");
245   ionLvlCmd->SetParameter(paramL);
246   paramL = new G4UIparameter("Lvl",'i',true);
247   paramL->SetDefaultValue("0.0");
248   ionLvlCmd->SetParameter(paramL);
249 
250   timeCmd = new G4UIcmdWithADoubleAndUnit("/gps/time",this);
251   timeCmd->SetGuidance("Set initial time of the particle.");
252   timeCmd->SetParameterName("t0",false,false);
253   timeCmd->SetDefaultUnit("ns");
254   //timeCmd->SetUnitCategory("Time");
255   //timeCmd->SetUnitCandidates("ns ms s");
256   
257   polCmd = new G4UIcmdWith3Vector("/gps/polarization",this);
258   polCmd->SetGuidance("Set polarization.");
259   polCmd->SetParameterName("Px","Py","Pz",false,false); 
260   polCmd->SetRange("Px>=-1.&&Px<=1.&&Py>=-1.&&Py<=1.&&Pz>=-1.&&Pz<=1.");
261 
262   numberCmd = new G4UIcmdWithAnInteger("/gps/number",this);
263   numberCmd->SetGuidance("Set number of particles to be generated per vertex.");
264   numberCmd->SetParameterName("N",false,false);
265   numberCmd->SetRange("N>0");
266 
267   // Verbosity
268   //
269   verbosityCmd = new G4UIcmdWithAnInteger("/gps/verbose",this);
270   verbosityCmd->SetGuidance("Set Verbose level for GPS");
271   verbosityCmd->SetGuidance(" 0 : Silent");
272   verbosityCmd->SetGuidance(" 1 : Limited information");
273   verbosityCmd->SetGuidance(" 2 : Detailed information");
274   verbosityCmd->SetParameterName("level",false);
275   verbosityCmd->SetRange("level>=0 && level <=2");
276 
277   volChkCmd = new G4UIcmdWithABool("/gps/checkVolume",this);
278   volChkCmd->SetGuidance("Switch on/off the check if the vertex position is inside the world volume.");
279   volChkCmd->SetGuidance("By default the check is on. There is a small performance gain if this check is off,");
280   volChkCmd->SetGuidance("but the user has to make sure setting the vertex position inside the world.");
281   volChkCmd->SetParameterName("switch",true,true);
282 
283   // Now extended commands
284   // Positional ones:
285   //
286   positionDirectory = new G4UIdirectory("/gps/pos/");
287   positionDirectory->SetGuidance("Positional commands sub-directory");
288 
289   typeCmd1 = new G4UIcmdWithAString("/gps/pos/type",this);
290   typeCmd1->SetGuidance("Sets source distribution type.");
291   typeCmd1->SetGuidance("Either Point, Beam, Plane, Surface or Volume");
292   typeCmd1->SetParameterName("DisType",false,false);
293   typeCmd1->SetDefaultValue("Point");
294   typeCmd1->SetCandidates("Point Beam Plane Surface Volume");
295 
296   shapeCmd1 = new G4UIcmdWithAString("/gps/pos/shape",this);
297   shapeCmd1->SetGuidance("Sets source shape for Plan, Surface or Volume type source.");
298   shapeCmd1->SetParameterName("Shape",false,false);
299   shapeCmd1->SetDefaultValue("NULL");
300   shapeCmd1->SetCandidates("Circle Annulus Ellipse Square Rectangle Sphere Ellipsoid Cylinder EllipticCylinder Para");
301 
302   centreCmd1 = new G4UIcmdWith3VectorAndUnit("/gps/pos/centre",this);
303   centreCmd1->SetGuidance("Set centre coordinates of source.");
304   centreCmd1->SetParameterName("X","Y","Z",false,false);
305   centreCmd1->SetDefaultUnit("cm");
306   // centreCmd1->SetUnitCandidates("micron mm cm m km");
307 
308   posrot1Cmd1 = new G4UIcmdWith3Vector("/gps/pos/rot1",this);
309   posrot1Cmd1->SetGuidance("Set the 1st vector defining the rotation matrix'.");
310   posrot1Cmd1->SetGuidance("It does not need to be a unit vector.");
311   posrot1Cmd1->SetParameterName("R1x","R1y","R1z",false,false); 
312   posrot1Cmd1->SetRange("R1x != 0 || R1y != 0 || R1z != 0");
313 
314   posrot2Cmd1 = new G4UIcmdWith3Vector("/gps/pos/rot2",this);
315   posrot2Cmd1->SetGuidance("Set the 2nd vector defining the rotation matrix'.");
316   posrot2Cmd1->SetGuidance("It does not need to be a unit vector.");
317   posrot2Cmd1->SetParameterName("R2x","R2y","R2z",false,false); 
318   posrot2Cmd1->SetRange("R2x != 0 || R2y != 0 || R2z != 0");
319 
320   halfxCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/halfx",this);
321   halfxCmd1->SetGuidance("Set x half length of source.");
322   halfxCmd1->SetParameterName("Halfx",false,false);
323   halfxCmd1->SetDefaultUnit("cm");
324   // halfxCmd1->SetUnitCandidates("micron mm cm m km");
325 
326   halfyCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/halfy",this);
327   halfyCmd1->SetGuidance("Set y half length of source.");
328   halfyCmd1->SetParameterName("Halfy",false,false);
329   halfyCmd1->SetDefaultUnit("cm");
330   // halfyCmd1->SetUnitCandidates("micron mm cm m km");
331 
332   halfzCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/halfz",this);
333   halfzCmd1->SetGuidance("Set z half length of source.");
334   halfzCmd1->SetParameterName("Halfz",false,false);
335   halfzCmd1->SetDefaultUnit("cm");
336   // halfzCmd1->SetUnitCandidates("micron mm cm m km");
337 
338   radiusCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/radius",this);
339   radiusCmd1->SetGuidance("Set radius of source.");
340   radiusCmd1->SetParameterName("Radius",false,false);
341   radiusCmd1->SetDefaultUnit("cm");
342   // radiusCmd1->SetUnitCandidates("micron mm cm m km");
343 
344   radius0Cmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/inner_radius",this);
345   radius0Cmd1->SetGuidance("Set inner radius of source when required.");
346   radius0Cmd1->SetParameterName("Radius0",false,false);
347   radius0Cmd1->SetDefaultUnit("cm");
348   // radius0Cmd1->SetUnitCandidates("micron mm cm m km");
349 
350   possigmarCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/sigma_r",this);
351   possigmarCmd1->SetGuidance("Set standard deviation in radial of the beam positional profile");
352   possigmarCmd1->SetGuidance(" applicable to Beam type source only");
353   possigmarCmd1->SetParameterName("Sigmar",false,false);
354   possigmarCmd1->SetDefaultUnit("cm");
355   // possigmarCmd1->SetUnitCandidates("micron mm cm m km");
356 
357   possigmaxCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/sigma_x",this);
358   possigmaxCmd1->SetGuidance("Set standard deviation of beam positional profile in x-dir");
359   possigmaxCmd1->SetGuidance(" applicable to Beam type source only");
360   possigmaxCmd1->SetParameterName("Sigmax",false,false);
361   possigmaxCmd1->SetDefaultUnit("cm");
362   // possigmaxCmd1->SetUnitCandidates("micron mm cm m km");
363 
364   possigmayCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/sigma_y",this);
365   possigmayCmd1->SetGuidance("Set standard deviation of beam positional profile in y-dir");
366   possigmayCmd1->SetGuidance(" applicable to Beam type source only");
367   possigmayCmd1->SetParameterName("Sigmay",false,false);
368   possigmayCmd1->SetDefaultUnit("cm");
369   // possigmayCmd1->SetUnitCandidates("micron mm cm m km");
370 
371   paralpCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/paralp",this);
372   paralpCmd1->SetGuidance("Angle from y-axis of y' in Para");
373   paralpCmd1->SetParameterName("paralp",false,false);
374   paralpCmd1->SetDefaultUnit("rad");
375   // paralpCmd1->SetUnitCandidates("rad deg");
376 
377   partheCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/parthe",this);
378   partheCmd1->SetGuidance("Polar angle through centres of z faces");
379   partheCmd1->SetParameterName("parthe",false,false);
380   partheCmd1->SetDefaultUnit("rad");
381   // partheCmd1->SetUnitCandidates("rad deg");
382 
383   parphiCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/pos/parphi",this);
384   parphiCmd1->SetGuidance("Azimuth angle through centres of z faces");
385   parphiCmd1->SetParameterName("parphi",false,false);
386   parphiCmd1->SetDefaultUnit("rad");
387   // parphiCmd1->SetUnitCandidates("rad deg");
388 
389   confineCmd1 = new G4UIcmdWithAString("/gps/pos/confine",this);
390   confineCmd1->SetGuidance("Confine source to volume (NULL to unset).");
391   confineCmd1->SetGuidance(" Usage: confine VolName");
392   confineCmd1->SetParameterName("VolName",false,false);
393   confineCmd1->SetDefaultValue("NULL");
394 
395   // Angular distribution commands
396   //
397   angularDirectory = new G4UIdirectory("/gps/ang/");
398   angularDirectory->SetGuidance("Angular commands sub-directory");
399 
400   angtypeCmd1 = new G4UIcmdWithAString("/gps/ang/type",this);
401   angtypeCmd1->SetGuidance("Sets angular source distribution type");
402   angtypeCmd1->SetGuidance(" Possible variables are: iso, cos, planar, beam1d, beam2d, focused or user");
403   angtypeCmd1->SetParameterName("AngDis",false,false);
404   angtypeCmd1->SetDefaultValue("iso");
405   angtypeCmd1->SetCandidates("iso cos planar beam1d beam2d focused user");
406 
407   angrot1Cmd1 = new G4UIcmdWith3Vector("/gps/ang/rot1",this);
408   angrot1Cmd1->SetGuidance("Sets the 1st vector for angular distribution rotation matrix");
409   angrot1Cmd1->SetGuidance(" Need not be a unit vector");
410   angrot1Cmd1->SetParameterName("AR1x","AR1y","AR1z",false,false);
411   angrot1Cmd1->SetRange("AR1x != 0 || AR1y != 0 || AR1z != 0");
412 
413   angrot2Cmd1 = new G4UIcmdWith3Vector("/gps/ang/rot2",this);
414   angrot2Cmd1->SetGuidance("Sets the 2nd vector for angular distribution rotation matrix");
415   angrot2Cmd1->SetGuidance(" Need not be a unit vector");
416   angrot2Cmd1->SetParameterName("AR2x","AR2y","AR2z",false,false);
417   angrot2Cmd1->SetRange("AR2x != 0 || AR2y != 0 || AR2z != 0");
418 
419   minthetaCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ang/mintheta",this);
420   minthetaCmd1->SetGuidance("Set minimum theta");
421   minthetaCmd1->SetParameterName("MinTheta",true,false);
422   minthetaCmd1->SetDefaultValue(0.);
423   minthetaCmd1->SetDefaultUnit("rad");
424   // minthetaCmd1->SetUnitCandidates("rad deg");
425 
426   maxthetaCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ang/maxtheta",this);
427   maxthetaCmd1->SetGuidance("Set maximum theta");
428   maxthetaCmd1->SetParameterName("MaxTheta",true,false);
429   maxthetaCmd1->SetDefaultValue(pi);
430   maxthetaCmd1->SetDefaultUnit("rad");
431   // maxthetaCmd1->SetUnitCandidates("rad deg");
432 
433   minphiCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ang/minphi",this);
434   minphiCmd1->SetGuidance("Set minimum phi");
435   minphiCmd1->SetParameterName("MinPhi",true,false);
436   minphiCmd1->SetDefaultValue(0.);
437   minphiCmd1->SetDefaultUnit("rad");
438   // minphiCmd1->SetUnitCandidates("rad deg");
439 
440   maxphiCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ang/maxphi",this);
441   maxphiCmd1->SetGuidance("Set maximum phi");
442   maxphiCmd1->SetParameterName("MaxPhi",true,false);
443   maxphiCmd1->SetDefaultValue(2.*pi);
444   maxphiCmd1->SetDefaultUnit("rad");
445   // maxphiCmd1->SetUnitCandidates("rad deg");
446 
447   angsigmarCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ang/sigma_r",this);
448   angsigmarCmd1->SetGuidance("Set standard deviation in direction for 1D beam.");
449   angsigmarCmd1->SetParameterName("Sigmara",false,false);
450   angsigmarCmd1->SetDefaultUnit("rad");
451   // angsigmarCmd1->SetUnitCandidates("rad deg");
452   
453   angsigmaxCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ang/sigma_x",this);
454   angsigmaxCmd1->SetGuidance("Set standard deviation in direction in x-direc. for 2D beam");
455   angsigmaxCmd1->SetParameterName("Sigmaxa",false,false);
456   angsigmaxCmd1->SetDefaultUnit("rad");
457   // angsigmaxCmd1->SetUnitCandidates("rad deg");
458 
459   angsigmayCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ang/sigma_y",this);
460   angsigmayCmd1->SetGuidance("Set standard deviation in direction in y-direc. for 2D beam");
461   angsigmayCmd1->SetParameterName("Sigmaya",false,false);
462   angsigmayCmd1->SetDefaultUnit("rad");
463   // angsigmayCmd1->SetUnitCandidates("rad deg");
464 
465   angfocusCmd = new G4UIcmdWith3VectorAndUnit("/gps/ang/focuspoint",this);
466   angfocusCmd->SetGuidance("Set the focusing point for the beam");
467   angfocusCmd->SetParameterName("x","y","z",false,false);
468   angfocusCmd->SetDefaultUnit("cm");
469   // angfocusCmd->SetUnitCandidates("micron mm cm m km");
470 
471   useuserangaxisCmd1 = new G4UIcmdWithABool("/gps/ang/user_coor",this);
472   useuserangaxisCmd1->SetGuidance("True for using user defined angular co-ordinates");
473   useuserangaxisCmd1->SetGuidance(" Default is false");
474   useuserangaxisCmd1->SetParameterName("useuserangaxis",true);
475   useuserangaxisCmd1->SetDefaultValue(false);
476 
477   surfnormCmd1 = new G4UIcmdWithABool("/gps/ang/surfnorm",this);
478   surfnormCmd1->SetGuidance("Makes a user-defined distribution with respect to surface normals rather than x,y,z axes.");
479   surfnormCmd1->SetGuidance(" Default is false");
480   surfnormCmd1->SetParameterName("surfnorm",true);
481   surfnormCmd1->SetDefaultValue(false);
482 
483   // Energy commands
484   //
485   energyDirectory = new G4UIdirectory("/gps/ene/");
486   energyDirectory->SetGuidance("Spectral commands sub-directory");
487 
488   energytypeCmd1 = new G4UIcmdWithAString("/gps/ene/type",this);
489   energytypeCmd1->SetGuidance("Sets energy distribution type");
490   energytypeCmd1->SetParameterName("EnergyDis",false,false);
491   energytypeCmd1->SetDefaultValue("Mono");
492   energytypeCmd1->SetCandidates("Mono Lin Pow Exp CPow Gauss Brem Bbody Cdg User Arb Epn LW");
493 
494   eminCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ene/min",this);
495   eminCmd1->SetGuidance("Sets minimum energy");
496   eminCmd1->SetParameterName("emin",false,false);
497   eminCmd1->SetDefaultUnit("keV");
498   // eminCmd1->SetUnitCandidates("eV keV MeV GeV TeV PeV");
499 
500   emaxCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ene/max",this);
501   emaxCmd1->SetGuidance("Sets maximum energy");
502   emaxCmd1->SetParameterName("emax",false,false);
503   emaxCmd1->SetDefaultUnit("keV");
504   // emaxCmd1->SetUnitCandidates("eV keV MeV GeV TeV PeV");
505 
506   monoenergyCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ene/mono",this);
507   monoenergyCmd1->SetGuidance("Sets a monocromatic energy (same as  gps/energy)");
508   monoenergyCmd1->SetParameterName("monoenergy",false,false);
509   monoenergyCmd1->SetDefaultUnit("keV");
510   // monoenergyCmd1->SetUnitCandidates("eV keV MeV GeV TeV PeV");
511 
512   engsigmaCmd1 = new G4UIcmdWithADoubleAndUnit("/gps/ene/sigma",this);
513   engsigmaCmd1->SetGuidance("Sets the standard deviation for Gaussian energy dist.");
514   engsigmaCmd1->SetParameterName("Sigmae",false,false);
515   engsigmaCmd1->SetDefaultUnit("keV");
516   // engsigmaCmd1->SetUnitCandidates("eV keV MeV GeV TeV PeV");
517 
518   alphaCmd1 = new G4UIcmdWithADouble("/gps/ene/alpha",this);
519   alphaCmd1->SetGuidance("Sets Alpha (index) for power-law energy dist.");
520   alphaCmd1->SetParameterName("alpha",false,false);
521   
522   tempCmd1 = new G4UIcmdWithADouble("/gps/ene/temp",this);
523   tempCmd1->SetGuidance("Sets the temperature for Brem and BBody distributions (in Kelvin)");
524   tempCmd1->SetParameterName("temp",false,false);
525 
526   ezeroCmd1 = new G4UIcmdWithADouble("/gps/ene/ezero",this);
527   ezeroCmd1->SetGuidance("Sets E_0 for exponential distribution (in MeV)");
528   ezeroCmd1->SetParameterName("ezero",false,false);
529 
530   gradientCmd1 = new G4UIcmdWithADouble("/gps/ene/gradient",this);
531   gradientCmd1->SetGuidance("Sets the gradient for Lin distribution (in 1/MeV)");
532   gradientCmd1->SetParameterName("gradient",false,false);
533 
534   interceptCmd1 = new G4UIcmdWithADouble("/gps/ene/intercept",this);
535   interceptCmd1->SetGuidance("Sets the intercept for Lin distributions (in MeV)");
536   interceptCmd1->SetParameterName("intercept",false,false);
537 
538   arbeintCmd1 = new G4UIcmdWithADouble("/gps/ene/biasAlpha",this);
539   arbeintCmd1->SetGuidance("Sets the power-law index for the energy sampling distri. )");
540   arbeintCmd1->SetParameterName("arbeint",false,false);
541 
542   calculateCmd1 = new G4UIcmdWithoutParameter("/gps/ene/calculate",this);
543   calculateCmd1->SetGuidance("Calculates the distributions for Cdg and BBody");
544 
545   energyspecCmd1 = new G4UIcmdWithABool("/gps/ene/emspec",this);
546   energyspecCmd1->SetGuidance("True for energy and false for momentum spectra");
547   energyspecCmd1->SetParameterName("energyspec",true);
548   energyspecCmd1->SetDefaultValue(true);
549 
550   diffspecCmd1 = new G4UIcmdWithABool("/gps/ene/diffspec",this);
551   diffspecCmd1->SetGuidance("True for differential and flase for integral spectra");
552   diffspecCmd1->SetParameterName("diffspec",true);
553   diffspecCmd1->SetDefaultValue(true);
554 
555   applyEnergyWeightCmd1 = new G4UIcmdWithABool("/gps/ene/applyEneWeight",this);
556   applyEnergyWeightCmd1->SetGuidance("Apply energy weight.");
557   applyEnergyWeightCmd1->SetGuidance("- Instead of using the Arb type histogram for sampling the energy spectrum,");
558   applyEnergyWeightCmd1->SetGuidance(" energy is sampled by Linear distribution, and the Arb type histogram is");
559   applyEnergyWeightCmd1->SetGuidance(" used for the weight of the generated particle.");
560   applyEnergyWeightCmd1->SetGuidance("- \"/gps/ene/type LW\" automatically applies this command.");
561   applyEnergyWeightCmd1->SetGuidance("- If this command has to be explicitly used, \"/gps/ene/type\" distribution mush be Lin.");
562   applyEnergyWeightCmd1->SetParameterName("flag",true);
563   applyEnergyWeightCmd1->SetDefaultValue(true);
564 
565   // Biasing + histograms in general
566   //
567   histDirectory = new G4UIdirectory("/gps/hist/");
568   histDirectory->SetGuidance("Histogram, biasing commands sub-directory");
569 
570   histnameCmd1 = new G4UIcmdWithAString("/gps/hist/type",this);
571   histnameCmd1->SetGuidance("Sets histogram type");
572   histnameCmd1->SetParameterName("HistType",false,false);
573   histnameCmd1->SetDefaultValue("biasx");
574   histnameCmd1->SetCandidates("biasx biasy biasz biast biasp biase biaspt biaspp theta phi energy arb epn");
575 
576   resethistCmd1 = new G4UIcmdWithAString("/gps/hist/reset",this);
577   resethistCmd1->SetGuidance("Reset (clean) the histogram ");
578   resethistCmd1->SetParameterName("HistType",false,false);
579   resethistCmd1->SetDefaultValue("energy");
580   resethistCmd1->SetCandidates("biasx biasy biasz biast biasp biase biaspt biaspp theta phi energy arb epn");
581 
582   histpointCmd1 = new G4UIcmdWith3Vector("/gps/hist/point",this);
583   histpointCmd1->SetGuidance("Allows user to define a histogram");
584   histpointCmd1->SetGuidance(" Enter: Ehi Weight");
585   histpointCmd1->SetParameterName("Ehi","Weight","Junk",true,true);
586   histpointCmd1->SetRange("Ehi >= 0. && Weight >= 0.");
587 
588   histfileCmd1 = new G4UIcmdWithAString("/gps/hist/file",this);
589   histfileCmd1->SetGuidance("Imports the arb energy hist in an ASCII file");
590   histfileCmd1->SetParameterName("HistFile",false,false);
591 
592   arbintCmd1 = new G4UIcmdWithAString("/gps/hist/inter",this);
593   arbintCmd1->SetGuidance("Sets the interpolation method for arbitrary distribution.");
594   arbintCmd1->SetGuidance("Spline interpolation may not be applicable for some distributions.");
595   arbintCmd1->SetParameterName("int",false,false);
596   arbintCmd1->SetDefaultValue("Lin");
597   arbintCmd1->SetCandidates("Lin Log Exp Spline");
598 }
599 
600 G4GeneralParticleSourceMessenger::~G4GeneralParticleSourceMessenger()
601 {
602   delete positionDirectory;
603   delete typeCmd1;
604   delete shapeCmd1;
605   delete centreCmd1;
606   delete posrot1Cmd1;
607   delete posrot2Cmd1;
608   delete halfxCmd1;
609   delete halfyCmd1;
610   delete halfzCmd1;
611   delete radiusCmd1;
612   delete radius0Cmd1;
613   delete possigmarCmd1;
614   delete possigmaxCmd1;
615   delete possigmayCmd1;
616   delete paralpCmd1;
617   delete partheCmd1;
618   delete parphiCmd1;
619   delete confineCmd1;
620 
621   delete angularDirectory;
622   delete angtypeCmd1;
623   delete angrot1Cmd1;
624   delete angrot2Cmd1;
625   delete minthetaCmd1;
626   delete maxthetaCmd1;
627   delete minphiCmd1;
628   delete maxphiCmd1;
629   delete angsigmarCmd1;
630   delete angsigmaxCmd1;
631   delete angsigmayCmd1;
632   delete angfocusCmd;
633   delete useuserangaxisCmd1;
634   delete surfnormCmd1;
635 
636   delete energyDirectory;
637   delete energytypeCmd1;
638   delete eminCmd1;
639   delete emaxCmd1;
640   delete monoenergyCmd1;
641   delete engsigmaCmd1;
642   delete alphaCmd1;
643   delete tempCmd1;
644   delete ezeroCmd1;
645   delete gradientCmd1;
646   delete interceptCmd1;
647   delete arbeintCmd1;
648   delete calculateCmd1;
649   delete energyspecCmd1;
650   delete diffspecCmd1;
651   delete applyEnergyWeightCmd1;
652 
653   delete histDirectory;
654   delete histnameCmd1;
655   delete resethistCmd1;
656   delete histpointCmd1;
657   delete histfileCmd1;
658   delete arbintCmd1;
659 
660   delete verbosityCmd;
661   delete volChkCmd;
662 
663   delete ionCmd;
664   delete ionLvlCmd;
665   delete particleCmd;
666   delete timeCmd;
667   delete polCmd;
668   delete numberCmd;
669   delete positionCmd;
670   delete directionCmd;
671   delete energyCmd;
672   delete listCmd;
673 
674   delete sourceDirectory;
675   delete addsourceCmd;
676   delete listsourceCmd;
677   delete clearsourceCmd;
678   delete getsourceCmd;
679   delete setsourceCmd;
680   delete setintensityCmd;
681   delete deletesourceCmd;
682   delete multiplevertexCmd;
683   delete flatsamplingCmd;
684 
685   delete gpsDirectory;
686   theInstance = nullptr;
687 }
688 
689 #define CHECKPG() { if (fParticleGun==nullptr) { \
690                       G4ExceptionDescription msg; \
691                       msg << "Command "<< command->GetCommandPath()<<"/";\
692                       msg << command->GetCommandName(); \
693                       msg << " used but no particle sources are set.";\
694                       msg <<" Add at least a source with: /gps/source/add.";\
695                       G4Exception("G4GeneralParticleSourceMessenger::SetNewValue","G4GPS003",\
696                                   FatalException,msg); return;\
697                   } }
698 
699 void G4GeneralParticleSourceMessenger::SetNewValue(G4UIcommand *command, G4String newValues)
700 {
701 //  if(command == typeCmd)
702 //    {
703 //      CHECKPG(); fParticleGun->GetPosDist()->SetPosDisType(newValues);
704 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
705 //             << " The command is obsolete and will be removed soon." << G4endl
706 //             << " Please try to use the new structured commands!" << G4endl;
707 //    }
708 //  else if(command == shapeCmd)
709 //    {
710 //      CHECKPG(); fParticleGun->GetPosDist()->SetPosDisShape(newValues);
711 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
712 //             << " The command is obsolete and will be removed soon." << G4endl
713 //             << " Please try to use the new structured commands!" << G4endl;
714 //    }
715 //  else if(command == centreCmd)
716 //    {
717 //      CHECKPG(); fParticleGun->GetPosDist()->SetCentreCoords(centreCmd->GetNew3VectorValue(newValues));
718 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
719 //             << " The command is obsolete and will be removed soon." << G4endl
720 //             << " Please try to use the new structured commands!" << G4endl;
721 //    }
722 //  else if(command == posrot1Cmd)
723 //    {
724 //      CHECKPG(); fParticleGun->GetPosDist()->SetPosRot1(posrot1Cmd->GetNew3VectorValue(newValues));
725 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
726 //             << " The command is obsolete and will be removed soon." << G4endl
727 //             << " Please try to use the new structured commands!" << G4endl;
728 //    }
729 //  else if(command == posrot2Cmd)
730 //    {
731 //      CHECKPG(); fParticleGun->GetPosDist()->SetPosRot2(posrot2Cmd->GetNew3VectorValue(newValues));
732 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
733 //             << " The command is obsolete and will be removed soon." << G4endl
734 //             << " Please try to use the new structured commands!" << G4endl;
735 //    }
736 //  else if(command == halfxCmd)
737 //    {
738 //      CHECKPG(); fParticleGun->GetPosDist()->SetHalfX(halfxCmd->GetNewDoubleValue(newValues));
739 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
740 //             << " The command is obsolete and will be removed soon." << G4endl
741 //             << " Please try to use the new structured commands!" << G4endl;
742 //    }
743 //  else if(command == halfyCmd)
744 //    {
745 //      CHECKPG(); fParticleGun->GetPosDist()->SetHalfY(halfyCmd->GetNewDoubleValue(newValues));
746 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
747 //             << " The command is obsolete and will be removed soon." << G4endl
748 //             << " Please try to use the new structured commands!" << G4endl;
749 //    }
750 //  else if(command == halfzCmd)
751 //    {
752 //      CHECKPG(); fParticleGun->GetPosDist()->SetHalfZ(halfzCmd->GetNewDoubleValue(newValues));
753 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
754 //             << " The command is obsolete and will be removed soon." << G4endl
755 //             << " Please try to use the new structured commands!" << G4endl;
756 //    }
757 //  else if(command == radiusCmd)
758 //    {
759 //      CHECKPG(); fParticleGun->GetPosDist()->SetRadius(radiusCmd->GetNewDoubleValue(newValues));
760 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
761 //             << " The command is obsolete and will be removed soon." << G4endl
762 //             << " Please try to use the new structured commands!" << G4endl;
763 //    }
764 //  else if(command == radius0Cmd)
765 //    {
766 //      CHECKPG(); fParticleGun->GetPosDist()->SetRadius0(radius0Cmd->GetNewDoubleValue(newValues));
767 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
768 //             << " The command is obsolete and will be removed soon." << G4endl
769 //             << " Please try to use the new structured commands!" << G4endl;
770 //    }
771 //  else if(command == possigmarCmd)
772 //    {
773 //      CHECKPG(); fParticleGun->GetPosDist()->SetBeamSigmaInR(possigmarCmd->GetNewDoubleValue(newValues));
774 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
775 //             << " The command is obsolete and will be removed soon." << G4endl
776 //             << " Please try to use the new structured commands!" << G4endl;
777 //    }
778 //  else if(command == possigmaxCmd)
779 //    {
780 //      CHECKPG(); fParticleGun->GetPosDist()->SetBeamSigmaInX(possigmaxCmd->GetNewDoubleValue(newValues));
781 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
782 //             << " The command is obsolete and will be removed soon." << G4endl
783 //             << " Please try to use the new structured commands!" << G4endl;
784 //    }
785 //  else if(command == possigmayCmd)
786 //    {
787 //      CHECKPG(); fParticleGun->GetPosDist()->SetBeamSigmaInY(possigmayCmd->GetNewDoubleValue(newValues));
788 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
789 //             << " The command is obsolete and will be removed soon." << G4endl
790 //             << " Please try to use the new structured commands!" << G4endl;
791 //    }
792 //  else if(command == paralpCmd)
793 //    {
794 //      CHECKPG(); fParticleGun->GetPosDist()->SetParAlpha(paralpCmd->GetNewDoubleValue(newValues));
795 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
796 //             << " The command is obsolete and will be removed soon." << G4endl
797 //             << " Please try to use the new structured commands!" << G4endl;
798 //    }
799 //  else if(command == partheCmd)
800 //    {
801 //      CHECKPG(); fParticleGun->GetPosDist()->SetParTheta(partheCmd->GetNewDoubleValue(newValues));
802 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
803 //             << " The command is obsolete and will be removed soon." << G4endl
804 //             << " Please try to use the new structured commands!" << G4endl;
805 //    }
806 //  else if(command == parphiCmd)
807 //    {
808 //      CHECKPG(); fParticleGun->GetPosDist()->SetParPhi(parphiCmd->GetNewDoubleValue(newValues));
809 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
810 //             << " The command is obsolete and will be removed soon." << G4endl
811 //             << " Please try to use the new structured commands!" << G4endl;
812 //    }
813 //  else if(command == confineCmd)
814 //    {
815 //      CHECKPG(); fParticleGun->GetPosDist()->ConfineSourceToVolume(newValues);
816 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
817 //             << " The command is obsolete and will be removed soon." << G4endl
818 //             << " Please try to use the new structured commands!" << G4endl;
819 //    }
820 //  else if(command == angtypeCmd)
821 //    {
822 //      CHECKPG(); fParticleGun->GetAngDist()->SetAngDistType(newValues);
823 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
824 //             << " The command is obsolete and will be removed soon." << G4endl
825 //             << " Please try to use the new structured commands!" << G4endl;
826 //    }
827 //  else if(command == angrot1Cmd)
828 //    {
829 //      CHECKPG();
830 //      G4String a = "angref1";
831 //      fParticleGun->GetAngDist()->DefineAngRefAxes(a,angrot1Cmd->GetNew3VectorValue(newValues));
832 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
833 //             << " The command is obsolete and will be removed soon." << G4endl
834 //             << " Please try to use the new structured commands!" << G4endl;
835 //    }
836 //  else if(command == angrot2Cmd)
837 //    {
838 //      CHECKPG();
839 //      G4String a = "angref2";
840 //      fParticleGun->GetAngDist()->DefineAngRefAxes(a,angrot2Cmd->GetNew3VectorValue(newValues));
841 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
842 //             << " The command is obsolete and will be removed soon." << G4endl
843 //             << " Please try to use the new structured commands!" << G4endl;
844 //    }
845 //  else if(command == minthetaCmd)
846 //    {
847 //      CHECKPG(); fParticleGun->GetAngDist()->SetMinTheta(minthetaCmd->GetNewDoubleValue(newValues));
848 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
849 //             << " The command is obsolete and will be removed soon." << G4endl
850 //             << " Please try to use the new structured commands!" << G4endl;
851 //    }
852 //  else if(command == minphiCmd)
853 //    {
854 //      CHECKPG(); fParticleGun->GetAngDist()->SetMinPhi(minphiCmd->GetNewDoubleValue(newValues));
855 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
856 //             << " The command is obsolete and will be removed soon." << G4endl
857 //             << " Please try to use the new structured commands!" << G4endl;
858 //    }
859 //  else if(command == maxthetaCmd)
860 //    {
861 //      CHECKPG(); fParticleGun->GetAngDist()->SetMaxTheta(maxthetaCmd->GetNewDoubleValue(newValues));
862 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
863 //             << " The command is obsolete and will be removed soon." << G4endl
864 //             << " Please try to use the new structured commands!" << G4endl;
865 //    }
866 //  else if(command == maxphiCmd)
867 //    {
868 //      CHECKPG(); fParticleGun->GetAngDist()->SetMaxPhi(maxphiCmd->GetNewDoubleValue(newValues));
869 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
870 //             << " The command is obsolete and will be removed soon." << G4endl
871 //             << " Please try to use the new structured commands!" << G4endl;
872 //    }
873 //  else if(command == angsigmarCmd)
874 //    {
875 //      CHECKPG(); fParticleGun->GetAngDist()->SetBeamSigmaInAngR(angsigmarCmd->GetNewDoubleValue(newValues));
876 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
877 //             << " The command is obsolete and will be removed soon." << G4endl
878 //             << " Please try to use the new structured commands!" << G4endl;
879 //    }
880 //  else if(command == angsigmaxCmd)
881 //    {
882 //      CHECKPG(); fParticleGun->GetAngDist()->SetBeamSigmaInAngX(angsigmaxCmd->GetNewDoubleValue(newValues));
883 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
884 //             << " The command is obsolete and will be removed soon." << G4endl
885 //             << " Please try to use the new structured commands!" << G4endl;
886 //    }
887 //  else if(command == angsigmayCmd)
888 //    {
889 //      CHECKPG(); fParticleGun->GetAngDist()->SetBeamSigmaInAngY(angsigmayCmd->GetNewDoubleValue(newValues));
890 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
891 //             << " The command is obsolete and will be removed soon." << G4endl
892 //             << " Please try to use the new structured commands!" << G4endl;
893 //    }
894 //  else if(command == useuserangaxisCmd)
895 //    {
896 //      CHECKPG(); fParticleGun->GetAngDist()->SetUseUserAngAxis(useuserangaxisCmd->GetNewBoolValue(newValues));
897 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
898 //             << " The command is obsolete and will be removed soon." << G4endl
899 //             << " Please try to use the new structured commands!" << G4endl;
900 //    }
901 //  else if(command == surfnormCmd)
902 //    {
903 //      CHECKPG(); fParticleGun->GetAngDist()->SetUserWRTSurface(surfnormCmd->GetNewBoolValue(newValues));
904 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
905 //             << " The command is obsolete and will be removed soon." << G4endl
906 //             << " Please try to use the new structured commands!" << G4endl;
907 //    }
908 //  else if(command == energytypeCmd)
909 //    {
910 //      CHECKPG(); fParticleGun->GetEneDist()->SetEnergyDisType(newValues);
911 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
912 //             << " The command is obsolete and will be removed soon." << G4endl
913 //             << " Please try to use the new structured commands!" << G4endl;
914 //    }
915 //  else if(command == eminCmd)
916 //    {
917 //      CHECKPG(); fParticleGun->GetEneDist()->SetEmin(eminCmd->GetNewDoubleValue(newValues));
918 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
919 //             << " The command is obsolete and will be removed soon." << G4endl
920 //             << " Please try to use the new structured commands!" << G4endl;
921 //    }
922 //  else if(command == emaxCmd)
923 //    {
924 //      CHECKPG(); fParticleGun->GetEneDist()->SetEmax(emaxCmd->GetNewDoubleValue(newValues));
925 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
926 //             << " The command is obsolete and will be removed soon." << G4endl
927 //             << " Please try to use the new structured commands!" << G4endl;
928 //    }
929 //  else if(command == monoenergyCmd)
930 //    {
931 //      CHECKPG(); fParticleGun->GetEneDist()->SetMonoEnergy(monoenergyCmd->GetNewDoubleValue(newValues));
932 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
933 //             << " The command is obsolete and will be removed soon." << G4endl
934 //             << " Please try to use the new structured commands!" << G4endl;
935 //    }
936 //  else if(command == engsigmaCmd)
937 //    {
938 //      CHECKPG(); fParticleGun->GetEneDist()->SetBeamSigmaInE(engsigmaCmd->GetNewDoubleValue(newValues));
939 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
940 //             << " The command is obsolete and will be removed soon." << G4endl
941 //             << " Please try to use the new structured commands!" << G4endl;
942 //    }
943 //  else if(command == alphaCmd)
944 //    {
945 //      CHECKPG(); fParticleGun->GetEneDist()->SetAlpha(alphaCmd->GetNewDoubleValue(newValues));
946 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
947 //             << " The command is obsolete and will be removed soon." << G4endl
948 //             << " Please try to use the new structured commands!" << G4endl;
949 //    }
950 //  else if(command == tempCmd)
951 //    {
952 //      CHECKPG(); fParticleGun->GetEneDist()->SetTemp(tempCmd->GetNewDoubleValue(newValues));
953 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
954 //             << " The command is obsolete and will be removed soon." << G4endl
955 //             << " Please try to use the new structured commands!" << G4endl;
956 //    }
957 //  else if(command == ezeroCmd)
958 //    {
959 //      CHECKPG(); fParticleGun->GetEneDist()->SetEzero(ezeroCmd->GetNewDoubleValue(newValues));
960 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
961 //             << " The command is obsolete and will be removed soon." << G4endl
962 //             << " Please try to use the new structured commands!" << G4endl;
963 //    }
964 //  else if(command == gradientCmd)
965 //    {
966 //      CHECKPG(); fParticleGun->GetEneDist()->SetGradient(gradientCmd->GetNewDoubleValue(newValues));
967 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
968 //             << " The command is obsolete and will be removed soon." << G4endl
969 //             << " Please try to use the new structured commands!" << G4endl;
970 //    }
971 //  else if(command == interceptCmd)
972 //    {
973 //      CHECKPG(); fParticleGun->GetEneDist()->SetInterCept(interceptCmd->GetNewDoubleValue(newValues));
974 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
975 //             << " The command is obsolete and will be removed soon." << G4endl
976 //             << " Please try to use the new structured commands!" << G4endl;
977 //    }
978 //  else if(command == calculateCmd)
979 //    {
980 //      CHECKPG(); fParticleGun->GetEneDist()->Calculate();
981 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
982 //             << " The command is obsolete and will be removed soon." << G4endl
983 //             << " Please try to use the new structured commands!" << G4endl;
984 //    }
985 //  else if(command == energyspecCmd)
986 //    {
987 //      CHECKPG(); fParticleGun->GetEneDist()->InputEnergySpectra(energyspecCmd->GetNewBoolValue(newValues));
988 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
989 //             << " The command is obsolete and will be removed soon." << G4endl
990 //             << " Please try to use the new structured commands!" << G4endl;
991 //    }
992 //  else if(command == diffspecCmd)
993 //    {
994 //      CHECKPG(); fParticleGun->GetEneDist()->InputDifferentialSpectra(diffspecCmd->GetNewBoolValue(newValues));
995 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
996 //             << " The command is obsolete and will be removed soon." << G4endl
997 //             << " Please try to use the new structured commands!" << G4endl;
998 //    }
999 //  else if(command == histnameCmd)
1000 //    {
1001 //      histtype = newValues;
1002 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
1003 //             << " The command is obsolete and will be removed soon." << G4endl
1004 //             << " Please try to use the new structured commands!" << G4endl;
1005 //    }
1006 //  else
1007 //      if(command == histpointCmd)
1008 //    {
1009 //      CHECKPG();
1010 //      if(histtype == "biasx")
1011 //        fParticleGun->GetBiasRndm()->SetXBias(histpointCmd->GetNew3VectorValue(newValues));
1012 //      if(histtype == "biasy")
1013 //        fParticleGun->GetBiasRndm()->SetYBias(histpointCmd->GetNew3VectorValue(newValues));
1014 //      if(histtype == "biasz")
1015 //        fParticleGun->GetBiasRndm()->SetZBias(histpointCmd->GetNew3VectorValue(newValues));
1016 //      if(histtype == "biast")
1017 //        fParticleGun->GetBiasRndm()->SetThetaBias(histpointCmd->GetNew3VectorValue(newValues));
1018 //      if(histtype == "biasp")
1019 //        fParticleGun->GetBiasRndm()->SetPhiBias(histpointCmd->GetNew3VectorValue(newValues));
1020 //      if(histtype == "biase")
1021 //        fParticleGun->GetBiasRndm()->SetEnergyBias(histpointCmd->GetNew3VectorValue(newValues));
1022 //      if(histtype == "theta")
1023 //        fParticleGun->GetAngDist()->UserDefAngTheta(histpointCmd->GetNew3VectorValue(newValues));
1024 //      if(histtype == "phi")
1025 //        fParticleGun->GetAngDist()->UserDefAngPhi(histpointCmd->GetNew3VectorValue(newValues));
1026 //      if(histtype == "energy")
1027 //        fParticleGun->GetEneDist()->UserEnergyHisto(histpointCmd->GetNew3VectorValue(newValues));
1028 //      if(histtype == "arb")
1029 //        fParticleGun->GetEneDist()->ArbEnergyHisto(histpointCmd->GetNew3VectorValue(newValues));
1030 //      if(histtype == "epn")
1031 //        fParticleGun->GetEneDist()->EpnEnergyHisto(histpointCmd->GetNew3VectorValue(newValues));
1032 //      G4cout << " G4GeneralParticleSourceMessenger - Warning: The command is obsolete and will be removed soon. Please try to use the new structured commands!" << G4endl;
1033 //    }
1034 //  else if(command == resethistCmd)
1035 //    {
1036 //      CHECKPG();
1037 //      if(newValues == "theta" || newValues == "phi") {
1038 //        fParticleGun->GetAngDist()->ReSetHist(newValues);
1039 //      } else if (newValues == "energy" || newValues == "arb" || newValues == "epn") {
1040 //        fParticleGun->GetEneDist()->ReSetHist(newValues);
1041 //      } else {
1042 //        fParticleGun->GetBiasRndm()->ReSetHist(newValues);
1043 //      }
1044 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
1045 //             << " The command is obsolete and will be removed soon." << G4endl
1046 //             << " Please try to use the new structured commands!" << G4endl;
1047 //    }
1048 //  else if(command == arbintCmd)
1049 //    {
1050 //      CHECKPG();
1051 //      fParticleGun->GetEneDist()->ArbInterpolate(newValues);
1052 //      G4cout << " G4GeneralParticleSourceMessenger - Warning:" << G4endl
1053 //             << " The command is obsolete and will be removed soon." << G4endl
1054 //             << " Please try to use the new structured commands!" << G4endl;
1055 //    }
1056 // else
1057        if( command==directionCmd )
1058     {
1059       CHECKPG();
1060       fParticleGun->GetAngDist()->SetAngDistType("planar");
1061       fParticleGun->GetAngDist()->SetParticleMomentumDirection(directionCmd->GetNew3VectorValue(newValues));
1062     }
1063   else if( command==energyCmd )
1064     {
1065       CHECKPG();
1066       fParticleGun->GetEneDist()->SetEnergyDisType("Mono");
1067       fParticleGun->GetEneDist()->SetMonoEnergy(energyCmd->GetNewDoubleValue(newValues));
1068     }
1069   else if( command==positionCmd )
1070     {
1071       CHECKPG();
1072       fParticleGun->GetPosDist()->SetPosDisType("Point");    
1073       fParticleGun->GetPosDist()->SetCentreCoords(positionCmd->GetNew3VectorValue(newValues));
1074     }
1075   else if(command == verbosityCmd)
1076     {
1077       fGPS->SetVerbosity(verbosityCmd->GetNewIntValue(newValues));
1078       // CHECKPG();
1079       // fParticleGun->SetVerbosity(verbosityCmd->GetNewIntValue(newValues));
1080     }
1081   else if( command==volChkCmd )
1082     {
1083       fGPS->CheckInside(volChkCmd->GetNewBoolValue(newValues));
1084     }
1085   else if( command==particleCmd )
1086     {
1087       if (newValues =="ion")
1088       {
1089         fShootIon = true;
1090       }
1091       else
1092       {
1093         fShootIon = false;
1094         G4ParticleDefinition* pd = particleTable->FindParticle(newValues);
1095         if(pd != nullptr)
1096         {
1097           CHECKPG();
1098           fParticleGun->SetParticleDefinition( pd );
1099         }
1100       }
1101     }
1102   else if( command==timeCmd )
1103     {
1104       CHECKPG();
1105       fParticleGun->SetParticleTime(timeCmd->GetNewDoubleValue(newValues));
1106     }
1107   else if( command==polCmd )
1108     {
1109       CHECKPG();
1110       fParticleGun->SetParticlePolarization(polCmd->GetNew3VectorValue(newValues));
1111     }
1112   else if( command==numberCmd )
1113     {
1114       CHECKPG();
1115       fParticleGun->SetNumberOfParticles(numberCmd->GetNewIntValue(newValues));
1116     }
1117   else if( command==ionCmd )
1118     {
1119       IonCommand(newValues);
1120     }
1121   else if( command==ionLvlCmd )
1122     {
1123       IonLvlCommand(newValues);
1124     }
1125   else if( command==listCmd )
1126     { 
1127       particleTable->DumpTable(); 
1128     }
1129   else if( command==addsourceCmd )
1130     {
1131       fGPS->AddaSource(addsourceCmd->GetNewDoubleValue(newValues));
1132     }
1133   else if( command==listsourceCmd )
1134     { 
1135       fGPS->ListSource();
1136     }
1137   else if( command==clearsourceCmd )
1138     { 
1139       fGPS->ClearAll();
1140       fParticleGun = nullptr;      
1141     }
1142   else if( command==getsourceCmd )
1143     { 
1144       G4cout << " Current source index:" << fGPS->GetCurrentSourceIndex() 
1145              << " ; Intensity:" << fGPS->GetCurrentSourceIntensity() << G4endl;
1146     }
1147   else if( command==setsourceCmd )
1148     { 
1149       // NOTE: This will also sets fParticleGun to the courrent source
1150       //       Not very clean, the GPS::SetCurrentSourceto will call:
1151       //       SetParticleSource( G4ParticleSource* )
1152       //       The point is that GPS has no public API to get a source by index
1153       // TODO: Can we add this API?
1154       const G4int sn = setsourceCmd->GetNewIntValue(newValues);
1155       if ( sn >= fGPS->GetNumberofSource() )
1156       {
1157         G4ExceptionDescription msg;
1158         msg << "Using command " << setsourceCmd->GetCommandPath() << "/"
1159             << setsourceCmd->GetCommandName() << " " << sn;
1160         msg << " is invalid " << fGPS->GetNumberofSource()
1161             << " source(s) are defined.";
1162         G4Exception("G4GeneralParticleSourceMessenger::SetNewValue",
1163                     "G4GPS005", FatalException, msg);
1164       }
1165       fGPS->SetCurrentSourceto(setsourceCmd->GetNewIntValue(newValues));
1166     }
1167   else if( command==setintensityCmd )
1168     { 
1169       fGPS->SetCurrentSourceIntensity(setintensityCmd->GetNewDoubleValue(newValues));
1170     }
1171   else if( command==deletesourceCmd )
1172     { 
1173       fGPS->DeleteaSource(deletesourceCmd->GetNewIntValue(newValues));
1174     }
1175   else if(command == multiplevertexCmd)
1176     {
1177       fGPS->SetMultipleVertex(multiplevertexCmd->GetNewBoolValue(newValues));
1178     }
1179   else if(command == flatsamplingCmd)
1180     {
1181       fGPS->SetFlatSampling(flatsamplingCmd->GetNewBoolValue(newValues));
1182     }
1183   //
1184   // new implementations
1185   //
1186   else if(command == typeCmd1)
1187     {
1188       CHECKPG();
1189       fParticleGun->GetPosDist()->SetPosDisType(newValues);
1190     }
1191   else if(command == shapeCmd1)
1192     {
1193       CHECKPG();
1194       fParticleGun->GetPosDist()->SetPosDisShape(newValues);
1195     }
1196   else if(command == centreCmd1)
1197     {
1198       CHECKPG();
1199       fParticleGun->GetPosDist()->SetCentreCoords(centreCmd1->GetNew3VectorValue(newValues));
1200     }
1201   else if(command == posrot1Cmd1)
1202     {
1203       CHECKPG();
1204       fParticleGun->GetPosDist()->SetPosRot1(posrot1Cmd1->GetNew3VectorValue(newValues));
1205     }
1206   else if(command == posrot2Cmd1)
1207     {
1208       CHECKPG();
1209       fParticleGun->GetPosDist()->SetPosRot2(posrot2Cmd1->GetNew3VectorValue(newValues));
1210     }
1211   else if(command == halfxCmd1)
1212     {
1213       CHECKPG();
1214       fParticleGun->GetPosDist()->SetHalfX(halfxCmd1->GetNewDoubleValue(newValues));
1215     }
1216   else if(command == halfyCmd1)
1217     {
1218       CHECKPG();
1219       fParticleGun->GetPosDist()->SetHalfY(halfyCmd1->GetNewDoubleValue(newValues));
1220     }
1221   else if(command == halfzCmd1)
1222     {
1223       CHECKPG();
1224       fParticleGun->GetPosDist()->SetHalfZ(halfzCmd1->GetNewDoubleValue(newValues));
1225     }
1226   else if(command == radiusCmd1)
1227     {
1228       CHECKPG();
1229       fParticleGun->GetPosDist()->SetRadius(radiusCmd1->GetNewDoubleValue(newValues));
1230     }
1231   else if(command == radius0Cmd1)
1232     {
1233       CHECKPG();
1234       fParticleGun->GetPosDist()->SetRadius0(radius0Cmd1->GetNewDoubleValue(newValues));
1235     }
1236   else if(command == possigmarCmd1)
1237     {
1238       CHECKPG();
1239       fParticleGun->GetPosDist()->SetBeamSigmaInR(possigmarCmd1->GetNewDoubleValue(newValues));
1240     }
1241   else if(command == possigmaxCmd1)
1242     {
1243       CHECKPG();
1244       fParticleGun->GetPosDist()->SetBeamSigmaInX(possigmaxCmd1->GetNewDoubleValue(newValues));
1245     }
1246   else if(command == possigmayCmd1)
1247     {
1248       CHECKPG();
1249       fParticleGun->GetPosDist()->SetBeamSigmaInY(possigmayCmd1->GetNewDoubleValue(newValues));
1250     }
1251   else if(command == paralpCmd1)
1252     {
1253       CHECKPG();
1254       fParticleGun->GetPosDist()->SetParAlpha(paralpCmd1->GetNewDoubleValue(newValues));
1255     }
1256   else if(command == partheCmd1)
1257     {
1258       CHECKPG();
1259       fParticleGun->GetPosDist()->SetParTheta(partheCmd1->GetNewDoubleValue(newValues));
1260     }
1261   else if(command == parphiCmd1)
1262     {
1263       CHECKPG();
1264       fParticleGun->GetPosDist()->SetParPhi(parphiCmd1->GetNewDoubleValue(newValues));
1265     }
1266   else if(command == confineCmd1)
1267     {
1268       CHECKPG();
1269       fParticleGun->GetPosDist()->ConfineSourceToVolume(newValues);
1270     }
1271   else if(command == angtypeCmd1)
1272     {
1273       CHECKPG();
1274       fParticleGun->GetAngDist()->SetAngDistType(newValues);
1275     }
1276   else if(command == angrot1Cmd1)
1277     {
1278       CHECKPG(); 
1279       G4String a = "angref1";
1280       fParticleGun->GetAngDist()->DefineAngRefAxes(a,angrot1Cmd1->GetNew3VectorValue(newValues));
1281     }
1282   else if(command == angrot2Cmd1)
1283     {
1284       CHECKPG(); 
1285       G4String a = "angref2";
1286       fParticleGun->GetAngDist()->DefineAngRefAxes(a,angrot2Cmd1->GetNew3VectorValue(newValues));
1287     }
1288   else if(command == minthetaCmd1)
1289     {
1290       CHECKPG(); 
1291       fParticleGun->GetAngDist()->SetMinTheta(minthetaCmd1->GetNewDoubleValue(newValues));
1292     }
1293   else if(command == minphiCmd1)
1294     {
1295       CHECKPG();
1296       fParticleGun->GetAngDist()->SetMinPhi(minphiCmd1->GetNewDoubleValue(newValues));
1297     }
1298   else if(command == maxthetaCmd1)
1299     {
1300       CHECKPG();
1301       fParticleGun->GetAngDist()->SetMaxTheta(maxthetaCmd1->GetNewDoubleValue(newValues));
1302     }
1303   else if(command == maxphiCmd1)
1304     {
1305       CHECKPG();
1306       fParticleGun->GetAngDist()->SetMaxPhi(maxphiCmd1->GetNewDoubleValue(newValues));
1307     }
1308   else if(command == angsigmarCmd1)
1309     {
1310       CHECKPG();
1311       fParticleGun->GetAngDist()->SetBeamSigmaInAngR(angsigmarCmd1->GetNewDoubleValue(newValues));
1312     }
1313   else if(command == angsigmaxCmd1)
1314     {
1315       CHECKPG();
1316       fParticleGun->GetAngDist()->SetBeamSigmaInAngX(angsigmaxCmd1->GetNewDoubleValue(newValues));
1317     }
1318   else if(command == angsigmayCmd1)
1319     {
1320       CHECKPG();
1321       fParticleGun->GetAngDist()->SetBeamSigmaInAngY(angsigmayCmd1->GetNewDoubleValue(newValues));
1322     }
1323   else if(command == angfocusCmd)
1324     {
1325       CHECKPG();
1326       fParticleGun->GetAngDist()->SetFocusPoint(angfocusCmd->GetNew3VectorValue(newValues));
1327     }
1328   else if(command == useuserangaxisCmd1)
1329     {
1330       CHECKPG();
1331       fParticleGun->GetAngDist()->SetUseUserAngAxis(useuserangaxisCmd1->GetNewBoolValue(newValues));
1332     }
1333   else if(command == surfnormCmd1)
1334     {
1335       CHECKPG();
1336       fParticleGun->GetAngDist()->SetUserWRTSurface(surfnormCmd1->GetNewBoolValue(newValues));
1337     }
1338   else if(command == energytypeCmd1)
1339     {
1340       CHECKPG();
1341       if(newValues=="LW")
1342       {
1343         fParticleGun->GetEneDist()->SetEnergyDisType("Lin");
1344         fParticleGun->GetEneDist()->SetGradient(0.);
1345         fParticleGun->GetEneDist()->SetInterCept(1.);
1346         fParticleGun->GetEneDist()->ApplyEnergyWeight(true);
1347       }
1348       else
1349       {
1350         fParticleGun->GetEneDist()->SetEnergyDisType(newValues);
1351         fParticleGun->GetEneDist()->ApplyEnergyWeight(false);
1352       }
1353     }
1354   else if(command == eminCmd1)
1355     {
1356       CHECKPG();
1357       fParticleGun->GetEneDist()->SetEmin(eminCmd1->GetNewDoubleValue(newValues));
1358     }
1359   else if(command == emaxCmd1)
1360     {
1361       CHECKPG();
1362       fParticleGun->GetEneDist()->SetEmax(emaxCmd1->GetNewDoubleValue(newValues));
1363     }
1364   else if(command == monoenergyCmd1)
1365     {
1366       CHECKPG();
1367       fParticleGun->GetEneDist()->SetMonoEnergy(monoenergyCmd1->GetNewDoubleValue(newValues));
1368     }
1369   else if(command == engsigmaCmd1)
1370     {
1371       CHECKPG();
1372       fParticleGun->GetEneDist()->SetBeamSigmaInE(engsigmaCmd1->GetNewDoubleValue(newValues));
1373     }
1374   else if(command == alphaCmd1)
1375     {
1376       CHECKPG();
1377       fParticleGun->GetEneDist()->SetAlpha(alphaCmd1->GetNewDoubleValue(newValues));
1378     }
1379   else if(command == tempCmd1)
1380     {
1381       CHECKPG();
1382       fParticleGun->GetEneDist()->SetTemp(tempCmd1->GetNewDoubleValue(newValues));
1383     }
1384   else if(command == ezeroCmd1)
1385     {
1386       CHECKPG();
1387       fParticleGun->GetEneDist()->SetEzero(ezeroCmd1->GetNewDoubleValue(newValues));
1388     }
1389   else if(command == gradientCmd1)
1390     {
1391       CHECKPG();
1392       fParticleGun->GetEneDist()->SetGradient(gradientCmd1->GetNewDoubleValue(newValues));
1393     }
1394   else if(command == interceptCmd1)
1395     {
1396       CHECKPG();
1397       fParticleGun->GetEneDist()->SetInterCept(interceptCmd1->GetNewDoubleValue(newValues));
1398     }
1399   else if(command == arbeintCmd1)
1400     {
1401       CHECKPG();
1402       fParticleGun->GetEneDist()->SetBiasAlpha(arbeintCmd1->GetNewDoubleValue(newValues));
1403     }
1404   else if(command == calculateCmd1)
1405     {
1406       CHECKPG();
1407       fParticleGun->GetEneDist()->Calculate();
1408     }
1409   else if(command == energyspecCmd1)
1410     {
1411       CHECKPG();
1412       fParticleGun->GetEneDist()->InputEnergySpectra(energyspecCmd1->GetNewBoolValue(newValues));
1413     }
1414   else if(command == diffspecCmd1)
1415     {
1416       CHECKPG();
1417       fParticleGun->GetEneDist()->InputDifferentialSpectra(diffspecCmd1->GetNewBoolValue(newValues));
1418     }
1419   else if(command == applyEnergyWeightCmd1)
1420     {
1421       CHECKPG();
1422       const auto& eDisType = fParticleGun->GetEneDist()->GetEnergyDisType();
1423       if(eDisType != "Lin")
1424       {
1425         G4ExceptionDescription ed;
1426         ed << "Energy distribution is defined as " << eDisType << ". /gps/ene/applyEneWeight is available only for Linear distribution.";
1427         command->CommandFailed(ed);
1428         return;
1429       }
1430       fParticleGun->GetEneDist()->ApplyEnergyWeight(applyEnergyWeightCmd1->GetNewBoolValue(newValues));
1431     } 
1432   else if(command == histnameCmd1)
1433     {
1434       histtype = newValues;
1435     }
1436   else if(command == histfileCmd1)
1437     {
1438       histtype = "arb";
1439       CHECKPG();
1440       fParticleGun->GetEneDist()->ArbEnergyHistoFile(newValues);
1441     }
1442   else if(command == histpointCmd1)
1443     {
1444       CHECKPG(); 
1445       if(histtype == "biasx")
1446         fParticleGun->GetBiasRndm()->SetXBias(histpointCmd1->GetNew3VectorValue(newValues));
1447       if(histtype == "biasy")
1448         fParticleGun->GetBiasRndm()->SetYBias(histpointCmd1->GetNew3VectorValue(newValues));
1449       if(histtype == "biasz")
1450         fParticleGun->GetBiasRndm()->SetZBias(histpointCmd1->GetNew3VectorValue(newValues));
1451       if(histtype == "biast")
1452         fParticleGun->GetBiasRndm()->SetThetaBias(histpointCmd1->GetNew3VectorValue(newValues));
1453       if(histtype == "biasp")
1454         fParticleGun->GetBiasRndm()->SetPhiBias(histpointCmd1->GetNew3VectorValue(newValues));
1455       if(histtype == "biaspt")
1456         fParticleGun->GetBiasRndm()->SetPosThetaBias(histpointCmd1->GetNew3VectorValue(newValues));
1457       if(histtype == "biaspp")
1458         fParticleGun->GetBiasRndm()->SetPosPhiBias(histpointCmd1->GetNew3VectorValue(newValues));
1459       if(histtype == "biase")
1460         fParticleGun->GetBiasRndm()->SetEnergyBias(histpointCmd1->GetNew3VectorValue(newValues));
1461       if(histtype == "theta")
1462         fParticleGun->GetAngDist()->UserDefAngTheta(histpointCmd1->GetNew3VectorValue(newValues));
1463       if(histtype == "phi")
1464         fParticleGun->GetAngDist()->UserDefAngPhi(histpointCmd1->GetNew3VectorValue(newValues));
1465       if(histtype == "energy")
1466         fParticleGun->GetEneDist()->UserEnergyHisto(histpointCmd1->GetNew3VectorValue(newValues));
1467       if(histtype == "arb")
1468         fParticleGun->GetEneDist()->ArbEnergyHisto(histpointCmd1->GetNew3VectorValue(newValues));
1469       if(histtype == "epn")
1470         fParticleGun->GetEneDist()->EpnEnergyHisto(histpointCmd1->GetNew3VectorValue(newValues));
1471     }
1472   else if(command == resethistCmd1)
1473     {
1474       CHECKPG(); 
1475       if(newValues == "theta" || newValues == "phi")
1476       {
1477         fParticleGun->GetAngDist()->ReSetHist(newValues);
1478       }
1479       else if (newValues == "energy" || newValues == "arb" || newValues == "epn")
1480       {
1481         fParticleGun->GetEneDist()->ReSetHist(newValues);
1482       }
1483       else
1484       {
1485         fParticleGun->GetBiasRndm()->ReSetHist(newValues);
1486       }
1487     }
1488   else if(command == arbintCmd1)
1489     {
1490       CHECKPG(); fParticleGun->GetEneDist()->ArbInterpolate(newValues);
1491     }
1492   else 
1493     {
1494       G4cout << "Error entering command" << G4endl;
1495     }
1496 }
1497 
1498 G4String G4GeneralParticleSourceMessenger::GetCurrentValue(G4UIcommand* command)
1499 {
1500   G4String cv;
1501   
1502   if( command==volChkCmd )
1503   { cv = volChkCmd->ConvertToString(fParticleGun->IfCheckInside()); }
1504   else
1505   {
1506   //  if( command==directionCmd )
1507   //  { cv = directionCmd->ConvertToString(fParticleGun->GetParticleMomentumDirection()); }
1508   //  else if( command==energyCmd )
1509   //  { cv = energyCmd->ConvertToString(fParticleGun->GetParticleEnergy(),"GeV"); }
1510   //  else if( command==positionCmd )
1511   //  { cv = positionCmd->ConvertToString(fParticleGun->GetParticlePosition(),"cm"); }
1512   //  else if( command==timeCmd )
1513   //  { cv = timeCmd->ConvertToString(fParticleGun->GetParticleTime(),"ns"); }
1514   //  else if( command==polCmd )
1515   //  { cv = polCmd->ConvertToString(fParticleGun->GetParticlePolarization()); }
1516   //  else if( command==numberCmd )
1517   //  { cv = numberCmd->ConvertToString(fParticleGun->GetNumberOfParticles()); }
1518   
1519   cv = "Not implemented yet";
1520   }
1521 
1522   return cv;
1523 }
1524 
1525 void G4GeneralParticleSourceMessenger::IonCommand(const G4String& newValues)
1526 {
1527   if (fShootIon)
1528   {
1529     G4Tokenizer next( newValues );
1530     // check argument
1531     fAtomicNumber = StoI(next());
1532     fAtomicMass = StoI(next());
1533     G4String sQ = next();
1534     if (sQ.empty())
1535     {
1536       fIonCharge = fAtomicNumber;
1537     }
1538     else
1539     {
1540       fIonCharge = StoI(sQ);
1541       sQ = next();
1542       if (sQ.empty())
1543       {
1544         fIonExciteEnergy = 0.0;
1545       }
1546       else
1547       {
1548         fIonExciteEnergy = StoD(sQ) * keV;
1549       }
1550     }
1551     G4ParticleDefinition* ion = G4IonTable::GetIonTable()
1552         ->GetIon(fAtomicNumber, fAtomicMass, fIonExciteEnergy);
1553     if (ion==nullptr)
1554     {
1555       G4ExceptionDescription ed;
1556       ed << "Ion with Z=" << fAtomicNumber;
1557       ed << " A=" << fAtomicMass << " is not defined";    
1558       ionCmd->CommandFailed(ed);
1559     }
1560     else
1561     {
1562       fParticleGun->SetParticleDefinition(ion);
1563       fParticleGun->SetParticleCharge(fIonCharge*eplus);
1564     }
1565   }
1566   else
1567   {
1568     G4ExceptionDescription ed;
1569     ed << "Set /gps/particle to ion before using /gps/ion command";
1570     ionCmd->CommandFailed(ed);
1571   }
1572 }
1573 
1574 void G4GeneralParticleSourceMessenger::IonLvlCommand(const G4String& newValues)
1575 {
1576   if (fShootIon)
1577   {
1578     G4Tokenizer next(newValues);
1579     // check argument
1580     fAtomicNumberL = StoI(next());
1581     fAtomicMassL = StoI(next());
1582     G4String sQ = next();
1583     if (sQ.empty())
1584     {
1585       fIonChargeL = fAtomicNumberL;
1586     }
1587     else
1588     {
1589       fIonChargeL = StoI(sQ);
1590       sQ = next();
1591       if (sQ.empty())
1592       {
1593         fIonEnergyLevel = 0;
1594       }
1595       else
1596       {
1597         fIonEnergyLevel = StoI(sQ);
1598       }
1599     }
1600 
1601     G4ParticleDefinition* ion = G4IonTable::GetIonTable()
1602         ->GetIon(fAtomicNumberL, fAtomicMassL, fIonEnergyLevel);
1603     if (ion == nullptr)
1604     {
1605       G4ExceptionDescription ed;
1606       ed << "Ion with Z=" << fAtomicNumberL;
1607       ed << " A=" << fAtomicMassL << " is not defined";
1608       ionLvlCmd->CommandFailed(ed);
1609     }
1610     else
1611     {
1612       fParticleGun->SetParticleDefinition(ion);
1613       fParticleGun->SetParticleCharge(fIonChargeL*eplus);
1614     }
1615 
1616   }
1617   else
1618   {
1619     G4ExceptionDescription ed;
1620     ed << "Set /gps/particle to ion before using /gps/ionLvl command";
1621     ionLvlCmd->CommandFailed(ed);
1622   }
1623 }
1624